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Dynamische Dunkle Energie und Variation funda- 
mentaler "Konstanten" 



Diese Arbeit beschaftigt sich mit den Auswirkungen moglicher Variationen funda- 

mentaler "Konstanten" auf den Prozcss dcr primordialen Nuklcosynthcsc (BBN). Die 
gewonnenen Ergebnisse zur Nukleosynthese werden mit Untersuchungen zu variieren- 
den Konstanten in anderen physikalischen Prozessen kombiniert, um Modelle der 
groBcn Vcrcinhcitlichung (GUT) und Quintessence zu iiberpriifen. Unsere Unter- 
suchungen ergeben, dass das ''Li-Problem der Nukleosynthese stark gemildert werden 
kann, sofern man Variationen von Konstanten zulasst, wobei insbesondere eine Vari- 
ation dcr Icichtcn Quarkmasscn cincn starkcn Einfluss hat. Wcitcrhin findcn wir, 
dass aktuelle Messungen zu variablen Konstanten im Rahmen von sechs exemplar- 
ischen GUT Modellen nicht miteinander und mit BBN in Einklang gebracht werden 
konncn, sofern eine monotone zcitlichc Variation angcnommcn wird. Wir folgern, 
dass aktuelle Messungen nichtverschwindender Variationen in starkem Widerspruch 
zueinander stehen und entweder selbst revidiert werden miissen, oder in der Natur 
crhcblich komplcxcrc GUT-Zusammcnhangc (und/odcr nicht-monotonc Variationen) 
vorliegen. Die im Rahmen dieser Dissertation vorgestellten Methoden erweisen sich 
hierbei als machtige Werkzeuge, um per Experiment unzugangliche Bereiche weit 
jcnscits dcs StandardmodcUs dcr Tcilchcnphysik bzw. dcs concordance ModcUs dcr 
Kosmologie auf ihre intrinsische Konsistenz sowie auch Vereinbarkeit miteinander zu 
iiberpriifen, sofern einmal erste unumstofiliche Beweise fiir Variationen von Naturkon- 
stanten vorliegen sollten. 



Dynamical dark energy and variation of fundamental 
"constants" 

In this thesis we study the influence of a possible variation of fundamental "constants" 
on the process of Big Bang Nucleosynthesis (BBN) . Our findings are combined with 
further studies on variations of constants in other physical processes to constrain mod- 
els of grand unification (GUT) and quintessence. We will find that the 'Li problem of 
BBN can be ameliorated if one allows for varying constants, where especially varying 
light quark masses show a strong influence. Furthermore, we show that recent studies 
of varying constants are in contradiction with each other and BBN in the frame- 
work of six exemplary GUT scenarios, if one assumes monotonic variation with time. 
We conclude that there is strong tension between recent claims of varying constants, 
hence either some claims have to be revised, or there are much more sophisticated 
GUT relations (and/or non- monotonic variations) realized in nature. The methods 
introduced in this thesis prove to be powerful tools to probe regimes well beyond the 
Standard Model of particle physics or the concordance model of cosmology, which 
are currently inaccessible by experiments. Once the first irrefutable proofs of varying 
constants are available, our method will allow for probing the consistency of models 
beyond the standard theories like GUT or quintessence and also the compatibility 
between these models. 
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Chapter 1 

Introduction 



The constants of nature 

Since the time of Newton, the constancy of the fundamental laws of nature has been 
undoubted. Comparing and reproducing experiments have been at the root of the 
scientific approach: A physical experiment which we perform today will have the same 
outcome as the same experiment performed tomorro Neglecting local gravitational 
effects, it should also not matter where we perform the experiment. Hence, it has been 
unquestionable for a long time that the laws of nature are constant over space and 
time. Moreover, Einstein formulated this space- and time independence of physics in 
his strong equivalence principle, making it an essential part of his theory of general 
relativity. 

Today's view of this question is somewhat different, at least from theoretical as- 
pects. Even though compelling evidence for changes in the laws of physics has up 
to now not been found, we have to admit that we are still lacking a profound test 
of this constancy. In the past, the laws of physics have only been thoroughly tested 
on time and length scales accessible by mankind, i.e. on timescales of years and on 
length scales that do not go beyond the size of our solar systerrH. Only recently as- 
trophysics and cosmology have opened a door to test physics on immensely broader 
scales, reaching out to unimaginable length scales of several gigaparsecs and going 
back in time to the very beginning of our Universe. 

This thesis will deal with probes of possible variations of constants throughout 
the whole accessible history of the Universe. In a first part, we will study one of 
the most distant (in time and space) events where physics can be applied and tested, 
primordial nucleosynthesis. It is the process during which the light elements of our 
Universe were formed and which happened when our Universe was only one minute 
old, extremely hot and dense. If physics was really subject to variations, primordial 
nucleosynthesis is a prime candidate for any studies of this kind. In a second step 
the obtained results will be combined together with further tests of varying constants 
at later times to derive a "history of variations" . Finally, we will show how these 
results can be used to test models beyond standard physics which currently cannot 
be accessed directly by experiments. 

^Neglecting experiments which incorporate probabilities, for instance quantum mechanical effects. 
^Note that general relativity has furthermore not been tested on length scales smaller than about 
1mm. 
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Outline 

In this thesis I wiU work out the influence of varying "constants" on the process 
of primordial nucleosynthesis and implied constraints to theories beyond standard 
physics. In the next chapter, I will give a short introduction to variations of physical 
constants, some historical remarks and theoretical motivations. Chapter [3] will lay the 
theoretical framework for our understanding of the Universe as a whole, explaining 
general relativity and the main laws of cosmology. In Chapter 0] I will introduce 
the concepts of supersymmetry and grand unified theories (GUTs) which are widely 
accepted as extensions of the Standard Model of particle physics. Chapter [5] will 
introduce quintessence models which can yield variations of constants. 

Part II will focus on the details of one particular process in the history of our Uni- 
verse, Big Bang Nucleosynthesis (BBN). Chapter [5] will explain the standard process 
of BBN and the physics behind. Chapter [7] will introduce the possibility of varying 
constants in the process of BBN, and Chapter [S] will demonstrate how one can relate 
the results to variations of the Standard Model parameters. Finally, in Chapter HI 
the observed element abundances will be used to derive constraints on variations of 
fundamental parameters. 

In Part III I will study relevant tests of varying constants from the Big Bang until 
today. Chapter 1101 gives an overview over tests of varying constants, and in Chapter 
mi l will combine these tests within six different GUT models, showing how variations 
of constants can in principle be used to probe models of grand unification. Using the 
six GUT models. Chapter [T2l shows how models of quintessence can be probed under 
the assumption of grand unification. 

Finally, in Chapter ll3l I will sum up the findings of this thesis and give some final 
conclusions and outlook. 

The work on this thesis has led to four main publications: 

• Michael Doran, Steffen Stern, Eduard Thommes, 
Baryon Acoustic Oscillations and Dynamical Dark Energy, 
JCAP 0704:015 (2007) |UST06j (not in focus of this thesis) 

• Thomas Dent, Steffen Stern, Christof Wettcrich, 

Primordial nucleosynthesis as a probe of fundamental physics parameters, 
Phys. Rev. D 76, 063513 (2007) jDSW07| 

• Thomas Dent, Steffen Stern, Christof Wetterich, 

Unifying cosmological and recent time variations of fundamental couplings, 
Phys. Rev. D 78, 103518 (2008) jDSWOS.l) 

• Thomas Dent, Steffen Stern, Christof Wetterich, 

Time variation of fundamental couplings and dynamical dark energy. 
Preprint .arXiv:0809.462 8, accepted by JCAP tDSW08.2] 



Chapter 2 

Variation of "constants 



2.1 The laws of physics and the constants of nature 

The fundamental laws of physics, represented by the Standard Model of particle 
physics and Einstein's theory of general relativity, consist of two parts. One part 
is the mathematical form of the laws {e.g. the l/r behavior of Newton's theory of 
gravity), the other part is the actual strength of the interactions relative to each other. 
Whilst the first part, the mathematical form of the laws of nature, can be derived 
from considerations of fundamental symmetries of natur^, the second part has to be 
put into the theories "by hand" in form of about 27 - from a theory standpoint a 
priori absolutely arbitrary - numerical values, the constants of nature. Tab. 12.11 gives 
a list of these fundamental constant^ for the Standard Model of particle physics and 
general relativit}|l. Note that the list of fundamental parameters gets much larger 
when going to theories beyond the Standard Model, e.g. supersymmetry (see Sec. 
14. 4p . Up to now it is unclear where these constants come from and if they are "real" 
constants in the sense that their numerical values are fixed once and for all. 

2.2 The question of constancy 

The question if the constants of nature are actually constant was probably first raised 
by Dirac [Dirac37l lDirac38, Dirac79| . In his "large numbers hypothesis", he argues 
that very large (or small) dimensionless constants must not enter in basic laws of 
physics. Based on his numerological principle, he suggests that very large numbers 
rather characterize the state of the Universe, specifically the time which has passed 
since the Big Bang. For instance, he finds that the age of the Universe in atomic 
time, Hoe'^/meC^ w 2 x 10~^^, is of the same order of magnitude as the ratio of 
electrostatic to gravitational force between proton and electron, Gf^mpirie/^^ « 

^For example, the Standard Model of particle physics is obtained when demanding a local 5(7(3) X 
SU{2) X [/(I) symmetry. See Sec. KT\ 

^As will be explained in Sec. 12.51 only ratios of masses are measurable fundamental parameters. 
Hence, in fact one can get rid of one the mass terms in Tab. 12.11 for instance by defining all masses 
with respect to the Planck mass. This would reduce the number of fundamental parameters by one. 

■^In cosmology some more free parameters turn up which have to be determined by observations, 
for instance those describing the composition of our Universe. However, it is assumed that these 
parameters can in principle be obtained from some fundamental laws of physics once the processes 
in the very early stage of the Universe are better understood. 
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Number of 


Type of constant 


parameters 


3 coupling constants 


3 


masses of 6 quarks 


6 


CKM matrix (3 angles + 1 complex phase) 


4 


masses of 3 leptons 


3 


Higgs mechanism 


2 


strong CP phase 


1 


masses of 3 neutrinos 


3 


PMNS mixing matrix for neutrinos 


4 


gravitational constant 


1 


in summa 


27 



Table 2.1: The fundamental constants of nature. 



4 X lO^'"'. Consequently, he suggests that also the latter quantity should vary with 
cosmic time. Attributing the variation to the gravitational sector, the intensity of 
all gravitational effects would then decrease with a rate of about 10~^° y~^. It 
was quickly found that this would lead to astrophysical effects |Chandrasekhar37] 
which could not be detected in the following time. Hence, Dirac's theory was finally 
abandoned, but the discussion on varying constants had startec0. 

In 1961, Brans and Dicke |BransDicke61| used Mach's principlelfl to derive what 
we now call a "scalar-tensor theory". In their model, the gravitational constant is 
replaced by a scalar field which can vary in space and time. Besides others, models of 
this kind are still being considered as theoretical arguments for variation of constants. 

2.3 Theoretical arguments for variation of constants 

In high-energy theories such as string theory, which unifies gravity with the Standard 
Model of particle physics, our low-energy laws of physics appear as an effective theory 
whose parameters are set dynamically by vacuum expectation values which break the 
"higher" symmetry. In particular string theory offers a plethora of possibilities to 
introduce variations of constants, for instance due to the fact that it is formulated 
with 10 (or 11) spacetime dimensions which need to be compactified in order to 
arrive at the 4 spacetime dimensions of the Standard Model (we will give some more 
details in Sec. 15. 5p . Similar considerations also apply to other theories with extra 
dimensions, for instance the possibility of varying constants in Kaluza-Klein theories 
has been studied in [Marciano83| . Hence, both temporal and spatial variations of 
constants are from a theoretical standpoint well founded, even though those high- 
energy theories mostly do not give any hint on the actual size of the variations. 

Also, "low-energy" theories, for instance theories which extend the concordance 
model of cosmology by introducing a cosmological scalar field, allow variations of 
constants. In this thesis we will concentrate on theories of coupled quintessence in 
which constants can depend on cosmic time and the environment. 

This thesis will examine the possibility of variations of constants from today back 

*See for instance |Uzan02] for a more complete review of the history of varying constant theories. 
^There are different formulations of Mach's principle. In Brans' and Dicke's argument it states 
|Brans05] that the gravitational constant should be a function of the mass distribution in the universe. 
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to the time of Big Bang Nucleosynthesis (BBN). During BBN, the composition of the 
Universe was quite different from today's composition (concerning temperature and 
pressure) . Hence, composition dependent effects which might cause spatial variations 
today might have caused variations at BBN time. However, since the Universe was 
almost homogeneous at BBN, these variations can effectively be treated as a time- 
dependent effect. This thesis will not evoke the question of space-dependence of 
constants but treat possible variations at BBN as purely temporal effects. 

2.4 Equivalence principles and possible violations 

Particle theory is based on Poincare covariance. In quantum field theory (QFT), we 
demand that each of the fundamental fields is a representation of the Poincare group. 
Hence, amongst others, invariance under spacetime translations is automatically built 
in. However, we can still implement spacetime variations by introducing additional 
dynamical fields, whose values are determined by the fields' own actions and their 
couplings to the rest of the theory. While the theory as a whole remains Poincare 
invariant, variations in measurable quantities can still arise if the solution for the 
additional fields has a nontrivial spacetime dependence. 

This discussion can be extended to general relativity (GR), which is also based 
on symmetry principles that are apparently violated by variations of constants. In 
particular, GR is based on the strong equivalence principle, which can be decomposed 
into the following symmetries 

• Weak equivalence principle: The trajectory of a freely falling test body only 
depends on its initial position and velocity and is independent of its composition. 

• Local Lorentz invariance: The outcomes of any experiments (whether gravita- 
tionally or not) in a laboratory moving in an inertial frame of reference are 
independent of the velocity of the laboratory. 

• Local position invariance: Outcomes of experiments (whether gravitationally or 
not) do not depend on their position in space and time. 

A space or time variation of fundamental constants obviously violates local position 
invariance. Also, as the gradient of any varying fundamental parameter defines a 
direction in spacetime, local Lorentz invariance is violated. Finally, it has been shown 
(see e.g. |NordtvedtO"2) ) that any space-time variation of fundamental constants will 
necessarily lead to an additional gravitational force, hence also the weak equivalence 
principle will be violatecH. As probes of general relativity so far do not find any 
violation of the theory, we can immediately conclude that variations of constants 
must be extremely tiny. Note, however, that GR has only been tested on relatively 
small time scales and also only on length scales from 1mm to the size of our solar 
system. 



°We will work out the relation between violation of the weak equivalence principle and variation 
of constants in Sec. 112.31 
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2.5 Variation of dimensionful parameters 

When measuring or estimating possible variations of constants, one always has to 
keep in mind that the variation of any dimensionful quantity is not physically well- 
defined, as one always has to specify how the dimension {e.g. the unit [Energy]) is 
defined. In general, a dimensionful quantity can only be measured by comparison with 
another dimensionful quantity, so in fact only dimensionless ratios are measurable. For 
example, measurements of variations of the electron mass rrie are only well-defined 
if one states how the mass unit is defined, for instance by choosing a system of 
units where the mass scale is kept fixed. Popular system of units are the "Einstein 
frame" where the Planck mass Mp is kept constant, or the "Jordan frame" where 
some particular particle physics scale is kept fixed. Considering the electron mass in 
the Einstein frame, the actually measured varying quantity (without system of units 
ambiguities) is then rather m^/AI-p. 

In the part of this thesis which deals with Big Bang Nucleosynthesis, we use 
a system of units where the QCD invariant scale Aqcd is kept constant. This is 
convenient for dealing with nuclear reactions, where the energy scales are mainly 
determined by the strong interaction. Thus the variations of dimensionful parameters 
include implicitly some power of Aqcd- For example, if we take the electron mass 
TOe as a varying parameter we are implicitly considering a variation of mg/ Aqcd- In 
the last part of this thesis we will work with theories of grand unification. There, 
the grand unified scale AIcut enters as natural scale which we choose to be constant. 
The appropriate conversion from a constant Aqcd to a constant AIcut system of 
units is explained in Sec. 14.6.31 

2.5.1 The chiral limit 

Many studies on varying constants work with the chiral limit, i.e. they assume that all 
quarks are massless |Epelbaum02[ |BeaneSavage02[ |Donoghue06| . Then all dimension- 
ful QCD parameters are simply proportional to a power of Aqc d , which ameliorates 
their treatment considerably. For instance, QCD masses like the proton mass simply 
scale like 

A In mp = A In Aqcd (2.1) 

and any other dimensionful QCD parameter according to its mass scale (for instance, 
cross sections with [a] ~ [Energy]"^ scale like A In cr = —2A In Aqcd)- Switching on 
the quark masses, one obtains a finite range for pion-mediated interactions, which 
may greatly affect static and dynamical properties of nuclei. Also, the masses of 
all hadrons are affected at some order in chiral perturbation theory [Gasser82] . In 
this thesis we will work with the full quark contributions, which are for most QCD 
parameters known at least in first order chiral perturbation theory, i. e. to terms linear 
in the quark masses. 

2.6 Probes of varying constants 

A possible variation of constants can be tested in various ways. Common tests are 
laboratory based measurements, for instance of atomic transitions. Also, a multi- 
tude of astrophysical and cosmological effects can be studied under the question of 
constancy, which allow to probe physics over a timescale unreachable with labora- 
tory measurements. In recent years probes of variations in the constants of nature 
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have been performed with increasingly high accuracy. Whilst direct laboratory mea- 
surements do not point towards any variation, some astrophysical tests yield slight 
variations. In part III (Sec. I10() we will list all recent relevant probes of varying con- 
stants, followed by detailed studies on how one can combine the different outcomes 
in unified scenarios. BBN as a probe of varying constants will be examined in detail 
in part II of this thesis. 

2.7 Fine-tuning of constants and the anthropic 
principle 

Connected to the question of constancy of fundamental constants is the question of 
fine-tuning of these constants. Even though this question can be seen as a rather 
philosophical one, we will shortly comment on it. 

As far as we know today, the value of most of the 27 fundamental constants is 
extremely fine-tuned in order to allow life to appear. It has been argued |Tegmark"97| 
that even small deviations (less than or order of 1%) will make the appearance of 
any life impossible. For example, if the strong force was slightly weaker, multi-proton 
nuclei would not be stable, and if it was slightly stronger, hydrogen could fuse into 
helium-2. Similar arguments can be found for the electromagnetic and weak force 
and for many other natural constants. 

This fine-tuning problem can be ameliorated, like all problems of this kind, by 
evoking the anthropic principle[l|. In short, this principle states that the Universe 
which we observe has to be capable to develop intelligent life like us. Otherwise we 
would not be here and could not ask the question why the Universe has exactly the 
laws of nature which it has. The final outcome is that the question why we are living in 
such a highly fine-tuned, i.e. extremely improbable, universe has simply disappeared, 
because the actual probability we have to discuss is rather the probability under the 
condition of our existence, which is no longer vanishingly small. 

In recent years scientists have come up with the idea of "multiverses" , stating 
that universes with many different kinds of physical properties are constantly formed 
[Linde86j . This is supported by candidate theories of everything (like string theory) 
which ab initio do not seem to have hard constraints which would exclusively se- 
lect our physics. Rather, they allow an extremely high number of different physical 
configurations. In the framework of those theories, universes with many different 
physical configurations bubble out constantly, and the anthropic principle states that 
our universe is the one of these many universes which allowed us to appear. 

These considerations are not directly connected to the investigations which are 
subject of this thesis, except the fact that varying constants would lead to an even 
more fine-tuned universe: Not only the values of the constants today, but also their 
whole time evolution would need to be tuned such that we could appear. We will not 
comment on the point of fine-tuning in the following, but it has become clear that the 
problems we are tackling have some deeper connection to philosophy and the question 
of why we are actually here. 



'^The concept of the anthropic principle was systematically introduced by Brandon Carter in a 
contribution to a symposium honoring Copernicus' 500th birthday in 1973 jCarte^l], even though 
the idea of the anthropic principle has already been used long before. 



Chapter 3 

Cosmology 



In this thesis we will consider probes for varying constants from today back to the 
first minute after the Big Bang. Hence it is essential to understand the evolution of 
our Universe from the Big Bang until today. This chapter gives a short review of 
our current picture of the Universe, its history and present status, and the important 
equations that govern its evolution. 



3.1 General relativity and the basics of cosmology 
3.1.1 General relativity 

General relativity is an extension of the theory of special relativity, which states that 
gravity is a purely geometric effect, generated by the curvature of spacetime. The 
relation between curvature and stress-energy is given by the Einstein field equations 

where -R^j^ is the Ricci tensor, R the Ricci scalar, g^i, the metric tensor and T^j, the 
stress-energy tensor. Equation (j3.ip is a complicated differential equation which can 
in general only be solved if one makes simplifying assumptions and/or uses numeric 
techniques. 



3.1.2 The basics of cosmology 

In cosmology one is interested in the evolution of the Universe as a whole. Thus, 
one usually confines oneself to physics on large scales which allows to make some 
simplifying assumptions that dramatically reduce the complexity of equation (|3.ip . 

Assumption 1. The main assumption o J cosmology is that the Universe is homoge- 
neous and isotropic on large scales. 

Of course, the existence of objects like the earth, sun etc. contradicts this as- 
sumption locally. However, if one averages over distances (> 1000 Mpc), it turns out 
that Assumption [T] is observationally well-justifiec(3. Demanding all quantities to be 

-"^Thc biggest known structure is the Sloan great wall which is 1.37 billion lightyears long. 
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Composition 


w 


non-relativistic matter 





ultra-relativistic matter (radiation) 


1/3 


curvature 


-1/3 


cosmological constant 


-1 



Table 3.1: Equation-of-state parameters for different types of cosmological compo- 
nents 



homogeneous and isotropic, one can show |WeinbergGRT| that the metric takes the 
simple fornix 

ds^ = dt^ - a{tf { Y^Tk^ + r'^dQ^ + sin^ ©^'/'^ j > (3-2) 

where k describes the curvature and a{t) is the scale parameter, related to the redshift 
z via 

^ (3.3) 



l + z 



The metric (I3.2p is called Friedmann- Robertson- Walker metric (FRW metric). The 
scale parameter fulfills the Friedmann equations 

(3.4) 

3- ^-IttGn fp+%V (3.5) 
a \ J 

where H is the Hubble constant and p and p denote the total energy and pressure 
density. These two densities are usually split up into the different components which 
are assumed to be present in today's Universe, baryonic and dark matter, dark energy 
(denoted with the symbol A), photons, neutrinos and curvatur^, 

P = Pb + PDA! + PK + P-i + pu + Pk ■ (3.6) 

The pressure is related to the energy via an equation of state. 

Pi = WiPi , (3.7) 

where the equation-of-state parameter Wi depends on the composition of the compo- 
nents as shown in Table 13.11 

With the critical density defined as 

37?'^ 

PC ■= , (3-8) 

07r{_rM 



^There are theories claiming that the averaged Einstein tensor G^i/ = R/iv — ^g^v which enters 
in Eq. I I3.1I I is not equivalent to the Einstein tensor derived from an averaged metric as given in Eq. 
1)3. 2| l. Since the actual outcome of these considerations is still unclear, we will not consider those 
theories in this thesis. See IBuchertOTI for a recent review. 

^In the very early Universe, also electrons will make a substantial contribution to the expansion 
rate. This applies to the epoch of BBN and will be explained in more detail in chapter \E\ 
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all densities are usually given as fractional densities 

r!, (3.9) 

PC 

Note that equation (j3.4|) yields 

QB + ^DM + ^A + ^-y + ^u + ^K^l (3.10) 

at all times. 

In the course of the evolution of the Universe, the energy densities scale like 

pcxa-3(i+-), (3.11) 

which means that the values of fli do not stay constant over time since we have 
different w for different kinds of energy densities fTab. [XT]) . As baryons and dark 
matter follow the same equation of state, one can combine these to the matter energy 
density 

■■= + ■ (3.12) 

Given today's values il?, one can combine Eqs. (j3.4p . (|3.8p . p.9p and (|3.1ip and use 
Tab. l3.1l to derive the time evolution of the Hubble constant, 

where we have neglected the neutrinos which have no substantial contribution to to- 
day's content of the Univers^ (see Tab. 13. 2p . Eq. (|3.13p shows that at early times 
(a ^ 1) non-relativistic and relativistic matter become dominant and any cosmolog- 
ical constant component irrelevant, whilst at late times (a ^ 1) J^a dominates. The 
flow of cosmological components in the ACDM concordance model (see Sec. 13. 2p is 
depicted in Fig. 13.11 where the time evolution of the fractional components is given 

by 

" r!5ja-4 + nia-^ + n^a-^ + nl ■ ^^-^^^ 

As can be seen in Fig. 13.11 today's Universe (z = 0) is dominated by dark energy 
(JIa) but did undergo 2 transitions, from radiation dominated to matter dominated 
and from matter to dark energy dominated: 



• In the early Universe, the expansion was almost completely due to relativistic 
particles radiation-dominated era. 

• At z w 5000, about 70,000 years after the Big Bang, we have matter-radiation- 
equality and the Universe becomes matter dominated. 

• At z « 0.4, about 4.3 Gyrs ago, the Universe becomes dominated by dark energy 
(in a ACDM model). 



^Further note that due to the tiny but non- vanishing mass of the neutrinos, the neutrino equation 
of state might change during the evolution of the universe. 
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z 



Figure 3.1: Evolution of the components in a ACDM Universe 

3.2 The concordance model: Our current picture of 
the Universe 

3.2.1 Historical development 

Presumably, the question of where we come from and where we will go is as old as 
mankind. As a first modern physical approach to questions of origin, evolution and 
fate of the Universe, one usually considers Einstein's paper "Cosmological Consider- 
ations in the General Theory of Relativity" from 1917 [EinsteinTT] . One might say 
that high-precision observational cosmology started with the Hubble space mission in 
1990. It was followed by further astrophysical and cosmological investigations, and 
basically all of these (mainly observational) tests point towards a coherent picture of 
our Universe, which is called the "concordance model". 

3.2.2 Our current picture of the Universe 

According to the concordance model, the Universe started in a Big Ban^ and has been 
expanding since then. All observational evidence points towards a so-called ACDM 
cosmology, stating that the Universe is geometrically flat {^Ik = 0) and consists 
besides known baryonic matter, leptons and photons of an unknown "dark matter" 
component which has the property of non-relativistic, only gravitationally interacting 
heavy particles, and a "dark energy" component, in the simplest version described by 

^Even though it is hoped that physics will once be able to explain the actual origin of this singular 
event, one is lacking an accepted theory of quantum gravitation which would allow to go beyond the 
time of the Planck epoch. 
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Figure 3.2: History of our Universe from Particle Data Group 2000 



a cosmological constant A. See for instance |Bartelmann06| for a compilation of the 
major observational evidences for the Big Bang and [WMAPS] for recent parameter 
determinations including all major probes of the concordance model. 

As can be seen in Fig. 13.11 we live in a dark-energy dominated universe just now 
but have undergone both radiation- and matter-dominated phases. The question why 
dark energy is taking over "just now" is unclear; this problem is called the "coincidence 
problem" or "why now problem" (see Sec. 15. ip . 

Looking back to the very beginning of our Universe, the concordance model sug- 
gests the following history of our Universe which is depicted in Fig. 13.21 Here, t is the 
age of the Universe, i.e. the time after the Big Bang, and T is the temperature of the 
Universe, defined by the photon temperaturqj. 

• t « 10-43 seconds (T « lO^^ GeV): 

Planck epoch, needs to be described by a quantum theory of gravity. 

• t « 10-35 seconds {T w lO^^ GeV): 

Inflation (exponentially fast expansion of the Universe), 
Baryogenesis (production of matter- antimatter asymmetry) 

®See Sec. 16.4.11 for details on the photon temperature. 
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• 10"^ seconds < t < 10^2 seconds (T « 0.1 GeV): 
Quark-hadron transition: protons and neutrons form 

• 1 second < t < 3 minutes, (T w 1 MeV, z w 10^°): 
Nucleosynthesis: light elements (D, He, Li) form 

• tfnTO, 000 years (T « 1 eV, z « 5000): 
Beginning of matter dominated era 

• t « 300,000 years (T w 0.25 eV, z w 1100): 
Recombination, the cosmic microwave background (CMB) forms 

• After that: 
Galaxy/star formation. 

Experimental tests of high-energy physics currently do not go beyond the TeV 
region. The processes which happened during the Planck epoch, inflation and also 
baryogenesis go beyond the physics of the Standard Model and are hence subject to 
a lot of speculation. Also the quark-hadron transition is hard to describe since it 
involves QCD at low energies, in the regime where QCD effects cannot be evaluated 
perturbatively. Hence, the oldest cosmological events which can be reasonably de- 
scribed quantitatively with known physics are primordial nucleosynthesis and CMB 
formation. 

3.3 Cosmological parameter values 

In recent years, important and quite extensive missions have been undertaken to 
deepen our understanding of cosmological relations. In particular, WMAP, SDSS and 
Supernovae la |WMAP5l IHinshawWMAPS] (see also references therein) have yielded 
a coherent set of cosmological parameters of a precision which had been inconceivable 
10 years ago. However, compared to the Standard Model of particle physics, the 
concordance model of cosmology is rather new and by far less tested. The set of 
cosmological parameters for a ACDM cosmology is given in Tab. 13.21 It turns out 
that only 4.6% of our Universe is made of "known" ordinary baryonic matter, the 
rest of the Universe is dark matter and dark energy. The energy composition of our 
Universe today is shown in Fig. 13.31 

Using the evolution equations of Sec. 13.11 the present-day values of cosmological 
parameters allow to deduce the content of our Universe in the past and also in the 
futurel!!. As observations also allow us to look back in time, the picture for the past 
is nowadays quite clear and observationally probed. The composition history of our 
Universe in a ACDM model is shown in Fig. 13.11 



^Cosmological models allow to extrapolate cosmology into the future, however models are not 
tested sufficiently to allow a definite prediction of what the ultimate fate of the Universe will be. 
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Quantity 


Symbol 


Value 


Hubble expansion rate 


Ho 


100/ikm/s Mpc"^ 


normalized Hubble expansion rate 


h 


0.701 ±0.013 


baryon density 


= Pb/Pc 


0.0462 ±0.0015 


dark matter density 


^drn — Pdrn / Pc 


0.233 ±0.013 


matter density 


— ^b ^" ^dm 


0.279 ±0.013 


dark energy density 


= Pa/Pc 


0.721 ±0.015 


radiation density 


^7 = Pl/Pc 


(5.0 ±0.2) • 10-5 


neutrino density 


= Piy/Pc 


< 0.013 (95% CL) 


baryon to photon ratio 




(6.21 ±0.16) • 10-1° 


CMB temperature 


T 


2.725 K 



Table 3.2: Parameters describing our Universe. WMAP "recommended parameter 
values" IWMAP51 [HinshawWMAPSj from WMAP5, BAO and SN for a ACDM cos- 
mology 




Figure 3.3: Today's energy content of our Universe. From NASA / WMAP Science 
Team 



Chapter 4 



The Standard Model and 
beyond 

4.1 The Standard Model of particle physics 

The Standard Model (SM) of particle physics describes the elementary particles and 
three of the four fundamental interactions, the strong, weak and electromagnetic in- 
teraction. This section will only give a rough overview over some particular aspects 
of the Standard Model which will be of relevance later. For a more comprehensive 
introduction, see for instance [HalzenMartinj or any other textbook on modern par- 
ticle physics. 

The list of the fundamental particles of the Standard Model comprises 

• six leptons {e, ^,T,Ve,i'^,Vr), 

• six quarks (u, d, s, c, b, t), 

• the gauge bosons as mediators of the fundamental interactions, 

• a Higgs boson. 

The matter particles enter as pointlike massless fermions, and the interactions are 
introduced by demanding a local SU{i) x SU{2) x C/(l) gauge symmetry. Here, the 
group SU{3) is responsible for quantum chromodynamics (QCD), with 8 massless 
gauge bosons called gluons as mediators of the strong force. At small momenta, 
the strong coupling constant becomes large (see Sec. 14. 2p . which is thought to be the 
explanation for confinement, i. e. the fact that only color-neutral particles are observed 
in nature. 

The theory of electroweak interaction goes back to the seminal work of Glashow, 
Salam and Weinberg |Weinberg67| ISalam69[ rGlashow70| (Nobel Prize 1979). It de- 
scribes the electroweak interaction by a SU{2) x U (!) gauge symmetry which is broken 
by the Higgs mechanism into the weak interaction (with massive gauge bosons 
and Z") and the electromagnetic sector with the photon as massless mediator of the 
electromagnetic force. In particular, the SU{2) x U{1) gauge symmetry implies four 
massless gauge bosons, written as W^, Wj^ and S^. Additionally, one introduces 
a scalar Higgs field cj) (as a weak doublet under SU{2) which has altogether 4 real 
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components) and gives it a potential V{(f>) which results in a vacuum that is not 
symmetric under the SU{2) x C/(l) gauge symmetry. This leads to a spontaneous 
symmetry breakdown of the electroweak symmetry, and the Higgs field obtains a 
nonzero vacuum expectation value (</>) « 246 GeV. One of the four Higgs compo- 
nents becomes a massive scalar particle, which is the only particle of the SM which 
has not yet been observed. The and a combination of Wj^ and i?^ obtain masses 
proportional to (0) and become the massive mediators of the weak force, and 
Z^, where the three remaining components of the Higgs form the longitudinal modes 
of the and Z^. The coupling constant ttem of the remaining electromagnetic 
symmetry group U{l)em can be obtained from the coupling constants of the original 
SU{2) X U{1) couphng constants ai, 012 (see £.17. [HalzenMartinj ). 

(^e^ - + ■ (4-1) 

Also, the SM fermions obtain masses via the Higgs mechanism, their mass is a product 
of Higgs v.e.v. and a Yukawa coupling /i^, for instance for the electron 

me = he{4>) . (4.2) 



4.2 Running of couplings 

The influence of fluctuations with different momenta leads to scale dependent coupling 
constants. See for instance [Wilson? ll |Wegner72l IWilson73| or any good textbook 
on quantum field theory for details of this process. Generally, physical systems at 
slightly different scales are described by the similar laws of physics, with slightly 
changed parameters. In quantum field theory, this behavior is described by the famous 
beta function, which describes the behavior of the coupling parameter g under slight 
changes of the energy scale /j,, 

= m • (4.3) 

Using the coupling constant <^ '•= ^ instead, one can also define a (3 function for a, 

= /?(«) ■ (4.4) 

The mathematical apparatus to investigate these changes of physical systems un- 
der scale transformations is called the renormalization group (RG). In quantum field 
theory, the renormalization group equation (j4.3p can only be computed perturba- 
tively as the exact RG equation would in principle include an infinite order of loop 
corrections. For our purpose the first-order (one-loop) RG equations are sufficient, 
and these are known for all three interactions of the Standard Model. In particular, 
the beta function for QED (with only photons and electrons present) at first order is 
given by 

2^2 

/3em(ae™) = ^ , (4.5) 



which is solved by 
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10"^ 10^ 10^2 1016 1020 104 108 1012 1016 1020 

Q (GeV) Q (GeV) 



Figure 4.1: Running of coupling constants of the three gauge groups SU (3) x SU{2) x 
U{1) in the Standard Model (left) and in SUSY (right), ai is scaled by a factor |, 
see Eq. From [Peskin97j . 



The fine structure constant a is defined in the limit of zero momentum transfer, i.e. 
for /i < nie. For QCD, the beta function is 

f3strongias) = -[n-^\^ , (4.7) 



asip-) = 7 — (4.8) 



solved by 



Here n/ is the number of quark flavors present, i.e. the number of quarks with mass 



l + ^(33-2n,)ln(-iL 



m 



< fj,. As rt/ < 6 in the Standard Model, the beta function (3strong is negative. 



The running of coupling constants is shown in Fig. 14.11 As the beta function for 
QCD is negative, the QCD coupling diverges when going to low energies. This effect, 
which was found by Wilczek, Politzer and Gross (Nobel price 2004), is thought to 
be the reason for confinement. As opposed to the electroweak theory, QCD thus has 
an intrinsic energy scale induced by the RG equation, the scale where as becomes 
formally infinite. Choosing rrig < /ip < rric such that n/ remains constant {uf = 3), 
this happens when the denominator in equation (|4.8p becomes zero, i.e. 



""(''°)(33-2n,.)lnfii)=-l, (4.9) 



which happens at the QCD invariant scale fi = Aqcd, defined by 

Aqcu Mocxp ( — - — — - — -) . (4.10) 
V(33 - 27i/)as(Aio)/ 

Hence Eq. (|4.8p can be rewritten as (/i < nic) 

6tt 

"^(^)=(33-2n,)WAQCZ.)- ^''^'^ 



Note that when the energy /i becomes of the order of Aqcd, perturbation theory 
breaks down, and a world of quarks and gluons becomes a world of pions, protons 
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Figure 4.2: Unification of tlic four forces in the string theory picture 

and so on. This is revealed in the fact that Aqcd is of the order of light meson masses 
|Berger06| , 

Aqcd ~ 200 MeV . (4.12) 

The beta functions given in Eqs. (|4.5p and (|4.7p are simplified versions. In the first 
case, Eq. (|4.5p only holds for one particle present, and in the second case one has to 
note that rif is not constant at all energies /i. In the full functions, particles only con- 
tribute when energies are above the particle's threshold energies, which are typically 
the corresponding particle masses. Hence every particle contributes one threshold 
term to the renormalization group equation. We will give the full renormalization 
group equations, including extra terms coming from additional supersymmetric par- 
ticles (see Sec. \4A^ in Sec. |4?5l 

4.3 The necessity of a "theory beyond" 

From the point of a theorist, the established Standard Model of particle physics cannot 
be the end of the story. One can definitely say that at latest at the Planck energy 
scale Mp « 10^^ GeV quantum gravity effects will become important, demanding a 
quantized description of gravity. A further hint that the Standard Model of particle 
physics might not be the end of the story is the running of the coupling constants. 
They seem to meet at a energy scale of Mqut ~ 10^^ GeV as depicted in Fig. 14.11 
Hence, it appears likely that electroweak and strong interaction can be unified in a 
grand unified theory (GUT). Within such grand unified theories, it is most likely that 
if any coupling constant of the Standard Model varies, all coupling constants will vary. 
In the later chapters, we will assume that some kind of GUT is realized, and hence the 
electroweak and strong coupling constants are related to each other. Further details 
on grand unified theories are given in Sec. 14.51 
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Actually, one can even go further with the idea of unification. String theory, for 
instance, implements an unification of the interactions of the Standard Model and 
also gravity. The unification of the four forces including gravity in a string theory 
picture is shown in Fig. 14.21 In such theories, parameters related to the mass of the 
SM particles {e.g. the Yukawa couplings) should also derive from some sort of unified 
greater theory, whereas at the level of the Standard Model and also at the level of 
simple GUTs, there is no direct relation to the gauge coupling sectoi0. In this thesis, 
we will assume constant Yukawa couplings. 

4.4 Supersymmetry and the MSSM 

Many high-energy theories {e.g. string theory) contain supersymmetry as an essen- 
tial part of the theory. Supersymmetry establishes a symmetry between bosons and 
fermions. Every boson gets a fermionic partner and vice versa with the same quantum 
numbers. For this thesis the motivation and theoretical framework of supersymmetry 
are not needed, so I will refrain from going into too much detail. Introductions to 
supersymmetry can be found in many textbooks and reviews, e.g. [NillesSUSYj . 

Within the Standard Model of particle physics, no supersymmetric partner can be 
found, hence in a supersymmetric extension of the Standard Model one has to intro- 
duce additional supersymmetric partners for every single particle of the SM. These 
supersymmetric partners are assumed to be heavier than the current experimentally 
tested mass region (« 100 GeV). 

The minimal supersymmetric extension of the Standard Model is called the MSSM 
(minimal supersymmetric standard model). It contains an additional supersymmetric 
partner for every SM particle. Furthermore, with supersymmetry, a single Higgs 
doublet would result in a gauge anomaly, so a second Higgs doublet is introduced. 
Hence the MSSM contains 2 additional (heavy) neutral Higgs scalars and two charged 
Higgs scalars, supplemented by the appropriate superpartners. 

The complete particle spectrum of the MSSM is given in Tab. 14.11 When working 
with supersymmetric theories in this thesis, we will always assume that the MSSM 
particle spectrum is realized. 

4.5 Grand unification 

In grand unified theories, the gauge group of the Standard Model SU{3)x SU{2)xU{l) 
with coupling constants gs, 92,91 is unified into a bigger Lie group {e.g. SU{5) or 
5*0(10)) with a single coupling constant gx at a certain energy scale, 

Mgut ~ 10^^ GeV , (4.13) 

which is assumed an independent parameter and can also vary with time. Its actual 
value depends on the specific form of the grand unified theory {e.g. SUSY/non-SUSY). 



^Due to renormalization group effects, afso the Yukawa couplings get contributions from coupling 
parameters. However, we will show in Sec. 111. T1 that these effects are small and can be neglected 
when studying variations of parameters. 
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SM particle 


El. charge 


Spin 


SUSY partner 


Spin 


quarks u,c,t 


2/3 


1/2 


squarks m, c, t 





quarks d,s,b 


-1/3 


1/2 


squarks d, s, h 





charged leptons e, /i, r 


-1 


1/2 


sleptons 





neutrinos Ve,v^i, Vt 





1/2 


sneutrinos (y) 





photon 7 





1 


photino 7 


1/2 







1 


Zino (Z) 


1/2 


neutral Higgs scalar 








Zinc Higgs 


1/2 




±1 


1 


Wine (T^) 


1/2 


charged Higgs scalar 


±1 





Wino Higgs 


1/2 


8 gluons 





1 


8 gluinos (g) 


1/2 


neutral Higgs H'^ 








higgsino 


1/2 


2 MSSM neutral Higgs 








2 neutral higgsinos 


1/2 


2 MSSM charged Higgs 


±1 





2 charged higgsinos 


1/2 



Table 4.1: The MSSM particle spectrum 



It can be shown (see e.g. |WeinbergQFT2| ) that for any unification of 5C/(3) x 
SU{2) X U{1) with couplings 33 = gs, (?2 and gi into a simple Lie group with coupling 
gx, one obtains the relation 

9% = 9 s = 9I = I?? (4.14) 

which holds a± E = Mqut- For the electromagnetic interaction, Eqs. (|4.ip and (|4.14p 
yield 

5 8 

a<7m = "2"^ + = otx^ + 3"3f^ = 3"^^ > (4-15) 

which actually only holds at E = Mqut but will be of relevance when studying the 
running of aem in a GUT framework. 

The value of the unified coupling ax can roughly be estimated from the RG 
running as displayed in Fig. 14.11 showing that is of the order 38 — 45 in the 
non-SUSY case and 23 — 29 in the SUSY case [Amaldi9T| . We take as representative 
values |Dent03| 

ax = 1/40 (non-SUSY) (4.16) 
ax = 1/24 (SUSY). (4.17) 

At lower energies, the GUT symmetry is broken and the relation (|4.14p does 
not hold any longer. The coupling constants of the SM evolve separately according 
to the renormalization group equations (see section r4.2|l . Generalizing the running 
of couplings as given in equations (|4.6p and (|4.8p to the full SM/MSSM particle 
spectrum, we obtain for QCD 

^ - ± ,n ^ - ± I„ {£-) ,4.18, 

where the first sum goes over all particles i with threshold mass m*^ < fi and the 
second sum goes over all particles with m*'' > /i. For the hi and 6*'* the values from 
Tab. 14.21 have to be applied. 
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Type of particle 




quarks 


2/3 


gluons 


-11/8 


squarks 


1/3 


gluinos 


1/4 



Table 4.2: Renormalization group coefficients for the strong interaction 



Type of particle 




chiral (or Majorana) fermion 


2/3 


complex scalar 


1/3 


vector boson 


-11/3 



Table 4.3: Renormalization group coefficients for the electromagnetic interaction 



The corresponding expression for the fine-structure constant a := aemifne) is 

-''A-x-lY.Q'etf'Ai^), (4.19) 



27r ^ \ Mgut 

th ^ 



where the factor | derives from Eq. (I4.15|) . Qem denotes the electric charge and for 



the values given in Tab. 14.31 are applied. An analogous equation also holds for 
the weak coupling, where the electric charge is replaced by the weak isospir(^. When 
dealing with weak interactions in this thesis, we will only be working with terms that 
contain the weak coupling in terms of the Fermi constanlH, 



V2 



J2 



As My/ = SjsMl with (0) the Higgs v.e.v., the weak coupling drops out and the 
Fermi constant can be expressed only in terms of the Higgs v.e.v.. 



4.6 Variations in a GUT framework 

The GUT relations which were introduced in the preceding section show that within 
a GUT framework, the coupling constants are usually related to further fundamental 
parameters, in particular the GUT coupling ax and threshold masses. In this section 
we will derive the equations which relate variations in the GUT coupling constant ax 
and particle masses to variations in the SM coupling constants. 



^For the SU(2) weak interaction, the weak isospin effectively acts like a multiplicative charge 
factor, hence it can be treated analogously to the electric charge. 
^See for instance the weak decay of the neutron, Sec. 18.1.31 
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4.6.1 Variation of the electromagnetic coupling 

For the MSSM particle spectrum, we obtain for the fine-structure constant from Eq. 

(TO 




+ --3-(l)^ln 

3 ^ ' Mgut 

4-2-(l)'l-ir^ + ^2.(lfln-^ + iha--i^ (4-22) 
3 Mgut 3 Mgut 3 Mgut _ 

where it has been used that 

• The Ught quarks u, d, s decouple at m*'' = A.qcd- 

• The quarks enter in 3 different colors 

• The charged leptons enter as both left- and righthanded particles (2 chiral 
fermions) 

• The massive gauge bosons have to be supplemented by a charged complex 
Higgs scalar (longitudinal DOFs) 

• m^± is the mass of the additional charged Higgs scalars which have to be 
introduced in MSSM. 

Taking the linear variation of Eq. (|4.22p . we obtain for the variation of the fine 
structure constant, including the MSSM particles, 

A In a 8 A In ax , 1/8 Aqcd 16., rricmt 4 rrib 
— H 1 -Am — 1 Am — ^ \ — Am 

a 3 ax 27r V3 Mgut 9 M^ut 9 Mgut 

4 . , mf,mi,m-r 21 . , Mw . , rh 1 . , mui \ 
-|--Aln Aln ,^^ -|-8Aln— + -Aln—- ^ . 4.23 

3 M3yj, 3 Mgut Mgut 3 Mgut J ^ ' 

Not knowing the actual mass or mass generating mechanism for the superpartners, we 
assume that the mechanism is the same for all superpartners and define rh as the av- 
erage superpartner mass. To obtain the corresponding relation in nonsupersymmetric 
models, one simply has to leave out the terms with rh and m//± When later dealing 
with variations in supersymmetric models, we will further assume A In m = A In m//± , 
so the last two terms can be combined into one. 



"'Note that the RG equations 114.221 1. 114.231 1 and also 114.241 1 and I I4.26I I only hold under the condition 
that all threshold masses are smaller than Mqut- Hence, the nonsupersymmetric case is not obtained 
in the limit msusy — ► oo. 
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4.6.2 Variation of the QCD scale 

As as diverges at ^ = Aqcd ~ 200 MeV, we are interested in the value of as(/^) 
in the regime Aqcd < fJ. < rric (in this regime, n/ = 3). For the MSSM particle 
spectrum, we obtain 

_i _i9 / n \ 1/2 mcmtmt 6 8 rrig 

[u] = o-v H in — In 5 — m — In ■ 



2-K \MgutJ 27r V3 M^[/t 3 Mqut 4 Mqut , 

(4.24) 

Inserting Eq. (|4.24p into Eq. (|4.10p {Aqcd < Ho = (J- < mc, nj = 3) yields 
A / \ 2/27 / X 2/9 



McuT V ^4c/T / \MgutJ 

and the linear variation gives 

Aln42££^^Alnax + AAln^^^^ + ?fAln^+Aln. ™^ ^ 



Mgc/t 9ax 27 M^^^j, 9 V Mgc/t Mgut J ' 

(4.26) 

When later dealing with variations in supersymmetric models, we will further assume 
Alnmg = AlnTOg, so the last two terms can be combined into one. 



4.6.3 Conversion of units 

As has been explained in Sec. 12.51 we will work with two different systems of units. 
During the discussion of BBN processes, we choose units with Aqcd = const, as the 
BBN energy scale is of roughly the same order of magnitude. When applying grand 
unified theories, AIqut is the more appropriate energy scale to keep constant. How- 
ever, we usually neglect the reference scale and write A In me instead for A In y^™^^ , 
for instance. The conversion to a different system is then performed by keeping track 
of all reference scales, 

A In = A In + A In . (4.27) 

AIgut ^qcd AIgut 

Obviously, the term A In ^^J^ naturally enters when converting from the Aqcd = 
const, to the Mgut = const, system of units, its explicit dependence on the unified 
coupling and particle masses is given in Eq. (|4.26p . Keep in mind, however, that 
the reference scale has to enter in the correct power, for instance the gravitational 
constant has units [Gn] — [Energy]^^, hence it enters as AIuGnAqj^^. 



Chapter 5 

Models of quintessence 



5.1 Problems of the cosmological constant 

Recent observations show that roughly 75% of the energy content of our Universe is 
made from dark energy (see Sec. 13. 2p . However, the nature of dark energy is still 
far from being clear. The assumption that the cosmological constant derives from a 
vacuum energy density sufferes from a severe fine-tuning problem. In particular, the 



oberseved dark energy density evaluates to Copeland06 



PA = ^ « 10-47 GeV\ (5.1) 
on 

while the vacuum energy density of particle physics which is evaluated by summing 
up the zero-point energies of the present quantum fields gives |Copeland06| 

Pvac^ ^^10'' GeY\ (5.2) 

Here we have chosen Mp as a natural cut-off scale where we assume that the known 
quantum field theory is no longer applicable. Obviously, there is a discrepancy of the 
order of 10^^^. Assuming that the dark energy comes from a particle physics origin, 
one would have to introduce counter terms which have to be extremely fine-tuned. 
Hence this problem is called the "finctuning" problem. 

A further problem which is related with dark energy can be seen in Fig. 13.11 In 
a ACDM model, the dark energy is only recently becoming important, and the time 
when the universe switched from a matter-dominated to a dark energy dominated 
epoch is only 4.3 billion years ago. There is no natural reason why these two pre- 
sumably completely independent constituents of our Universe are of about the same 
order of magnitude and / or why we live in a period of time where this is the case. 
This problem is called the "coincidence problem" or "why now problem" , and typi- 
cally cosmological models with a cosmological constant (like the ACDM model) fail 
to address this issue. 



5.2 Basics of quintessence 

Models of quintessence jWetterichSS.H IRatraPeebles88| can offer an explanation to 
the issues mentioned in the previous section. A good review on dark energy models 
can be found in I Copeland06| . 
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In quintessence theories, one introduces a scalar field (called the cosmon) which 
is coupled to gravity and, most times, also to matter and gauge fields. A typical 
Lagrangian for a quintessence theory including couplings to matter and the electro- 
magnetic gauge field looks like [Wetterich02.1[ |Copeland06| 

£ = M^R + ^(dip)^ + V{ip) - V^rn + C,^ (5.3) 
with gauge field coupling 

C,,^ = -\{l + K,^f)-^F^,F^'•' (5.4) 

and matter term 



V^,n = me(93)ee + mu{ip)uu + ... + mdark(v)^dark^'dark + 



(5.5) 



Here Afp denotes the reduced Planck mass 

1 



Mp 



/8n 



(5.6) 



If there are only slight changes in the cosmon field, the dependence of the mass, m((p), 
can be linearized, i.e. 

m{(p) = (1 + A(p)mo (5.7) 

with some coupling A. However, there are also models with significant changes in the 
masses, for instance in the models of growing neutrinos in Sec. 15.41 where we apply a 
more advanced nonlinear expression. 

In order to derive one of the main properties of quintessence, its capability of 
producing accelerated expansion, we can neglect couplings to matter and gauge fields 
and work instead with the action 



S 



(5.8) 



In the background of a flat FRW cosmology (Sec. I3.1.2p . and assuming that ip is 
homogeneous, i.e. it only depends on time, a variation of the action ()5.8p with respect 
to (fi yields the equation of motion 



(fi + 3Htp + — = . 
dip 

The corresponding energy momentum tensor 

^ _ ~2 6S 



(5.9) 



(5.10) 



yields the energy and pressure densities 



(5.11) 
(5.12) 
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Then Eqs. ([23) and yield the relations 

l^^ + V{^) (5.13) 

■ - V{^)] , (5.14) 

showing that we get an accelerating universe when ip'^ < V{(p). Introducing the kinetic 
energy T := </5^/2, we define the equation of state parameter of quintessence, 

In the next section we introduce two specific models of quintessence, crossover quint- 
essence which has been introduced 20 years ago |Wetterich88.1j and a very recent 
model, where quintessence is strongly coupled to neutrinos. 

As the cosmon also couples to other fields and matter, one question one might 
ask is whether the cosmon evolution decouples in a local cluster with high "cosmon 
charge density" from the cosmological evolution. It has been shown jWetterich02.H 
rMota03t ,Shaw05 that for a very light field weakly coupled to matter the local per- 
turbations are generally small relative to the cosmological evolution. In other words, 
the evolution of the scalar field in a cluster of galaxies or on Earth does not decouple 
from the cosmological evolution (in distinction to the gravitational field), such that 
its cosmological time evolution is reflected in a universal variation of couplings, both 
on Earth and in the distant Universe. 



5.3 Crossover quintessence models 

Our first class of models is "crossover quintessence" |HebcckerOO' 'Doran07' ' Wett erich02 .2] . 
Here the scalar field follows tracking solutions [Wetterich88 . 1 , , RatraPcebles88j at large 
redshift. In this early epoch the equation of state Wh is equal to that of the dominant 
energy component (matter or radiation). One particular difference to cosmological 
constant models is that this type of quintessence models yields a non vanishing amount 
of early dark energy, flh,e 0- Typically, such models have an exponential potential, 
for instance 

V{ip) = M^e^"^ . (5.16) 

In this specific case, a is related to the early dark energy fraction jWetterich88.H 
lAmendolaOS] . 

^Ke = (5.17) 

with n = 3(4) for the matter (radiation) epoch. Late-time acceleration can be 
achieved, for instance, by slight modifications in the cosmon potential or, equiva- 
lently, in the kinetic term [HebeckerOOj . 

At some intermediate redshift before the onset of acceleration, the time evolution 
of the cosmon slows down. In consequence, there is a crossover to a negative equation 
of state and the fraction of energy density due to the scalar begins to grow. In recent 
epochs the field has an effective equation of state Wh > —1. The aim of this thesis is 



^The effects of early dark energy on tfie measurements of baryon acoustic peaks have been studied 
by the author in |DST06I . However, these considerations are not subject of this thesis. 
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Figure 5.1: Equation of state of quintessence in the crossover quintessence model 
{who = -0.99, z+ = 7) 



not building and solving models of this type in detail, but rather estimating general 
properties of the scalar evolution. Hence we will not start with appropriate potentials 
for crossover quintessence and evolve the cosmon over time, but rather simulate the 
behavior of the field by defining the quintessence equation of state by hand. We set the 
dark energy equation of state to constant at late times with value Who- Above some 
given redshift the equation of state crosses over to the scaling condition Wh = 
in the matter dominated era; then for z > Zgq, before matter-radiation equality, we 
again have scaling through = 1/3, where be obtained via 



(5.18) 



Then the general relation (a = (1 + z) ^) 

rfln p 



din a 



-3(l + w(a)) (5.19) 



may be used to find the matter, radiation and dark energy densities over cosmological 
time. Combining Eq. ()5.15p with Eq. ()5.11|) , we can estimate the scalar kinetic energy 
via 

ipy2^ Phil + Wh)/2 (5.20) 

and thus integrate dip /da = (p/aH from the present back to any previous redshift. 
The initial conditions are set by specifying the present densities of matter, radiation 
and dark energy and the model parameters Who and z+. For illustration, we set 
WfiQ — —0.99 and z+ — 7. The resulting equation of state is displayed in Fig. 15.11 
the corresponding evolution of energy components in Fig. 15.21 and the dimensionless 
cosmon field ip/Mp in Fig. 15.31 As can be seen from Fig. 15.31 in this type of models 
the scalar field has a monotonic evolution. Assuming a constant coupling S to the 
fundamental varying parameter, usually ax, the variation is given by 

Alnaxiz) ^ 6iip{z) - ifiO)) . (5.21) 



Hence this ansatz implies a monotonic evolution of variations. 
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Figure 5.3: Dimcnsionless cosmon field f in the crossover quintessence model 
(wha = -0.99, 2+ = 7) 
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5.4 Growing neutrino mass models 

Growing neutrino models |Amendola08[ I WetterichOS] explain the value of today's dark 
energy density by the "principle of cosmological selection" . The present fraction of 
dark energy, fi^, is set by a dynamical mechanism. The essential ingredient of this 
class of models is a neutrino mass that depends on the cosmon field if and grows in the 
course of the cosmological evolution. As soon as the neutrinos become non-relativistic, 
their coupling to the cosmon triggers an effective stop (or substantial slowing) of the 
evolution of the cosmon. Before this event, the quintessence field follows the tracking 
behavior described in the preceding section. In the models which we will study, the 
cosmon is assumed to have the potential from Eq. (j5.16p . 

V{ip) = M^e^"^ . 

The present dark energy density, pho , can be expressed in terms of the average present 
neutrino mass, TOiy(io)i ^^nd a dimensionless parameter C of order unity [AmendolaOSj . 

(p.o)^/^ = 1.07 (^^^) lO-^eV . (5.22) 

We follow again our simple proportionality assumption, namely that the cosmon 
coupling to a typical fundamental parameter is given by Eq. (|5.2ip . 

Alnax{z) = 6iip{z)-ipiO)) , 

with a proportional variation for other couplings according to the unification scenario 
that we will study. This is the only contribution to the variation of the unified 
coupling ax and Mgut/Mp. However, a new ingredient is an additional variation 
of the Higgs v.e.v. (c^) with respect to Mgut, which only becomes relevant at late 
time |Wetterich08| . It is due to the effect of a changing weak triplet operator on the 
v.e.v. of the Higgs doublet. If the dominant contribution to the neutrino mass arises 
from the "cascade mechanism" (or "induced triplet mechanism" ) via the expectation 
value of this triplet, this changing triplet value is directly related to the growing 
neutrino mass [Wetterich08| . To understand this mechanism, we start with the most 
general mass matrix for the light neutrinos, 

= MdMj^^mI + Ml ■ (5.23) 

The first term is responsible for the seesaw mechanism [Minkowski??) with the mass 
matrix for heavy "right handed" neutrinos Mb. and a Dirac mass term Mjj. The 
second term accounts for the "induced triplet mechanism" |Magg80| 

Ml^^-^, (5.24) 

where a heavy SU (2)L-triplet field t with mass Mt enters the equation (see |Wetterich08j 
for details). It is assumed that the mass of the triplet depends on the cosmon field, 
Mt = Mtiifi). 

The (^-dependence of the Higgs v.e.v. (0) is introduced by assuming a general 
effective potential U{ip, cj), t). Solving the field equations for the Higgs doublet field (f) 
and the triplet field t, dU/d4> = 0, dU/dt — 0, the cosmon potential is then obtained 
as 

F(^) = [/(^,0(^),i(^)) . (5.25) 
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In |Wetterich08j the simple potential 

A 



(5.26) 



is assumed, with 7 and A some coupling parameters. Solving the field equations for 
the Higgs doublet, it is found |Wetterich08| 



^GUT 



^"GUT 



XM^iip) 



(5.27) 



where {(j))^ has to be chosen such that the measured Higgs v.e.v. is obtained today. 

In the following we consider two models, with slightly different functional depen- 
dence of the Higgs v.e.v. and neutrino mass on the scalar field. 



5.4.1 Stopping growing neutrino model 

In the first model, the cosmon asymptotically approaches a constant value ("stopping 
growing neutrino model" ) [WetterichOSj and the neutrino mass is given by 



With a triplet mass dependence 



Mi 



1 — exp 



1 — exp 



-if -ft) 



the additional Higgs variation is given according to 



Mgut 



{z)^H{l-R{z)) 



-0.5 



where 



R{z) 



Rq 

l-exp(--^((^(z)-^t))' 



(5.28) 



(5.29) 



(5.30) 



(5.31) 



Here, ipt ~ 27.6 is the asymptotic value (choosing the parameter a = 10 in the 
exponential potential [WetterichOS] ). For illustration we take the set of parameters 
given in ^Wet terichOS] , e = —0.05, and H is set by demanding the Higgs v.e.v. being 
consistent with measurements today, {(f)){z = 0) = 175 GeV. We set i?o ~ 10^^, 
however in general we only require R{z = 0) <C 1. The resulting variations are shown 
in Fig. [O] and Fig. [531 

The stopping growing neutrino model has an oscillation in (0) that grows both 
in frequency and amplitude at late times as ip approaches its asymptotic value. Such 
oscillations must not be too strong as measurements between z — 2 and today would 
measure a high rate of change. The oscillation may be made arbitrarily small by 
choosing small Rq. However, the linear variation (j5.2ip is independent of Rq. 
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Figure 5.4: Evolution of the dimensionless cosmon field in the stopping growing neu- 
trino model of [WetterichOS] . 




Figure 5.5: Additional variation of the Higgs v.e.v. according to Eq. (|5.30p in the 
stopping growing neutrino model of ^WetterichOS] . 
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Figure 5.6: Evolution of the dimensionless cosmon field in the scaling growing neutrino 
model of [AmendolaOS] . 
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Figure 5.7: Additional variation of the Higgs v.e.v. according to Eq. (|5.30p in the 
scaling growing neutrino model of |Amendola08j . 



5.4.2 Scaling growing neutrino model 



The second growing neutrino model [AmendolaOS] does not lead to an asymptotically 
constant ip. Now the coupling of the neutrino to the cosmon ip is given by a constant 
f3, according to 

TO^ = m^e"^"^. (5.32) 

This "scaling growing neutrino model" leads in the future to a scaling solution with a 
constant ratio between the neutrino and cosmon contributions to the energy density. 

With the choice of parameters (3 = —52, a = 10 and m^.o = 2.3 eV |Amendol a08) . 
and given the triplet mechanism of [WetterichOS] . the Higgs v.e.v. varies as Eq. (|5.30p . 
where now R is given by 

i?(z) = i?oe~'''^(^^ (5.33) 

Here the Higgs oscillations remain comparatively small in amplitude, while the abso- 
lute value of {(p) grows overall with time: see Figs. 15.61 and 15.71 with the parameter 
choice i?o ~ lO^*'. Compared to the "stopping growing neutrino model", this model 
has milder oscillations at late time. 
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5.5 A short note on string theory 

String theory is a candidate for a "theory of everything" . The underlying concept 
is that fundamental particles are one dimensional objects (strings) which live in a 
10 dimensional spacetime. As the world as we know it has only four spacetime di- 
mensions, the remaining 6 dimensions are usually compactified in a way that only 4 
"large" spacetime dimensions remain. During compactification, several scalar fields 
appear which can in principle have all sorts of couplings to the Standard Model fields. 
A further scalar field, the dilaton, is present in string theory from the very beginning. 
Its v.e.v. sets the string theory coupling constant. 

With a multitude of possible fields and couplings, one should in principle be able 
to model any quintessence scenario within the framework of string theory. However, 
there are strong arguments against it, in particular based on the small mass and 
potential of the quintessence field (see e.g. [BanksOl| V A deeper study of variations 
implied by string theory or possible tests of string theory is out of scope of this thesis, 
but in principle our methods will also apply to any string theory induced variations 
of constants, as long as they can be described with effective field theories. Then our 
methods also allow to constrain the allowed regions in the landscape of string theory. 
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Bang Nucleosynthesis 
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Chapter 6 

Big Bang Nucleosynthesis 



6.1 Why BBN? 

As we have mentioned in Sec. 12. 3[ changes in fundamental constants are most Hkely 
to appear over large time scales and / or different environments. Big Bang Nuc- 
leosynthesis is the earliest process in the history of the Universe which can be both 
reasonably described with standard physics and astrophysically probed. Also, the 
Universe was much denser at this timq^, even though locally the environment one 
finds for instance in heavy stars, supernovae or black holes is much more extreme 
than at BBN. Hence, it is very reasonable to carefully study the influence of varying 
parameters on the process of primordial element production. 



6.2 How will we study BBN? 

BBN is a complex process involving a lot of nuclear reactions and particles. There are 
two main approaches to a theoretical prediction of primordial element abundances. 
The first one is purely analytical and only gives rough estimates, whereas the sec- 
ond one numerically simulates the whole process of BBN and gives high-precision 
abundance predictions. 

For our purpose, analytical estimates will turn out to have insufficient accuracy, 
hence we will utilize a numerical procedure to obtain our findings on the influence of 
varying constants on BBN. Numerical codes for the simulation of primordial nucleo- 
synthesis go back to the late 60s and 70s |Peeble s66', 'Wagoner69' 'W agoner72] . Our 
BBN code is based on the Kawano 1992 code jKawano88, Kawano92 , with updated 
nuclear reaction cross sections as given in [NACRE99] and |NETGEN] . We have im- 
proved the code in terms of numerical accuracy in order to be able to derive high 
precision parameter dependences (see Sec. l6.7.T] for details). In the following sections 
we explain the process of BBN and the underlying physics in greater detail. As has 
been explained in Sec. 12.51 we will work in a system of units where Aqcd is kept 
constant. 



^The energy density at BBN was roughly the same as the the density of normal water. 
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Figure 6.1: Network of main reactions responsible for primordial nucleosynthesis. 
From [Bartelmann06| . 

6.3 The process of BBN 

In the standard BBN process (SBBN) , the element synthesis does not depend on any 
pre-BBN phase. At a temperature of, for instance, T = lO^^K, baryogenesis and 
quark condensation have ended and the Universe only contains protons and neutrons 
in the baryonic sector. These are held in equilibrium via the weak reactions 



n ■ 



n + Vf, 



(6.1) 



As the Universe expands and cools down, these reactions freeze out at T w 800 
keV, and the neutrons decay freely. The next step in nucleosynthesis is the fusion of 
deuterium, which happens when reaction 



P 



drops out of equilibrium. Due to the large photodissociation cross section, any pro- 
duced deuterium is immediately destroyed by photons from the background radiation 
until their temperature has dropped considerably below the binding energy of deu- 
terium, T Ki Bd — 2.2 MeV. In fact, since photons are much more abundant than 
baryons, the high-energy photons in the Maxwell tail keep the reaction in equilibrium 
until the temperature reaches T « 80 keV. Once deuterium fusion sets in, elements 
up to mass number A=7 are synthesized via a network of 11 main nuclear reactions 
displayed in Fig. 16.11 As the Universe continues expanding and cooling down, the 
temperature and density will fall below the level required for the nuclear reactions at 
some point. This happens a few minutes after the Big Bang, when the reaction rates 
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Figure 6.2: Element abundances F as a function of the decreasing temperature of the 
universe Tg^Tx IQ-^K^^ 



become slower than the expansion rate of the Universe, 

n{av) < H , (6.2) 

with n the particle number density and (av) the Maxwell-Boltzmann averaged cross 
section (see Sec. 16.4. 5| . At the end of the BBN process, the element composition of 
the Universe is roughly 75% hydrogen, 25% helium (numbers w.r.t. mass) and small 
amounts of deuterium and ''Li. Two elements produced during BBN, ''Be and tritium, 
decay via a slow /3-decay to ^Li and ^He respectively after BBN has ended. In fact, 
the primordial ''Li abundance primarily derives from ^Bc produced during BBN. 

The evolution of element abundances during the process of BBN is shown in 
Fig. 16. 2[ where we define the mass abundance Yi for a nucleus i as 

Y,:=A^ (6.3) 

with Ub the baryon number density, Ui the number density of nucleus i and Ai its 
mass number. 



Beyond A = 8 elements 

During primordial nucleosynthesis, essentially no elements are produced beyond the 
A = 7 mass limit. Heavier elements, including carbon and oxygen which are ubiqui- 
tous on earth, have been produced later in stars and were thrown back into space 
during supernova explosions. The process leading to those elements proceeds over the 
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unstable nucleus ®Be via the "triple- alpha process" 

^Hc + ^He ^Be + 7 

«Bc + 4He ^ '^C + j. 

Due to the low density and low temperature at BBN (compared to stars), the prob- 
ability that a ^Be nucleus meets a ^He nucleus during its lifetime of 67 x lO^^^s is 
extremely low. Hence, the absence of stable A = 8 nuclei prohibits any element fusion 
beyond ^Li during BBN. 



6.4 The physics of BBN 

In this section we will describe the underlying processes of BBN in greater detail and 
introduce the relevant equations. 



6.4.1 Cosmological background equations 

BBN is described in the framework of the FRW metric. The basic equations describing 
the expansion of the Universe are the two Friedmann equations (|3.4[) and (|3.5p . The 
process of BBN happened at the a redshift of z '--^ 10^°, a time when the Universe 
was radiation dominated (see Fig. 13. ip . i.e. the expansion rate of the Universe is 
completely controlled by relativistic particles. The temperature T of the Universe is 
defined by the photon temperature Tj, and the Stefan-Boltzmann law can be used to 
express energy and number density of the photons in terms of T-y, 



2C(3) f keT^ 
he 



P. - -TzhT,' (6.4) 



(6.5) 



BBN is taking place in the temperature regime 3 x lO^K < T < 0.1 x lO^K. 
Combining Eq. (|6.4p with the Friedmann equation (|3.4p yields a relation between 
temperature and time after the Big Bang, 

^ = C'bT) . (6.6) 



5^3 c5 

This relation can be applied in the early Universe until the number of photons in the 
Universe is increased by e+e^ annihilation. 

Before BBN, neutrinos, photons and electrons were in equilibrium, for instance 

via 

v + i? < y e+ + e~ < > 2j . (6.7) 

This reaction freezes out at T w 1.5 MeV |Fornengo97| , but as also the neutrinos are 
ultra relativistic, both Tj and depend on the scale factoi@ 

T-i (6.8) 

a 



^Equation l|6.8p also gives a relation between temperature and redshift, T ~ {1 + z). Note however 
that this relation does not hold during reheating. 
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and consequently should evolve equally over time. Due to e+e^ annihilation, the 
photon temperature increases whilst the neutrino temperature strictly follows the 
T ~ i law. With T^ and T^, evolving differently, the two temperatures have to be 
tracked separately during BBN. For the neutrino temperature we get an expression 
equivalent to the Stefan-Boltzmann law, 

where N^ is the number of neutrino generations (in SBBN N^, = 3). 

Note that temperatures can be given both in units of Kelvin (K) and energy (MeV) 
via E = kgT , with the conversion factor 

1 MeV ~ 11.6 • lO^K (6.10) 



At T « 10i°K {EKil MeV) the reaction e+ + «-> 27 drops out of thermodynamic 
equilibrium. Electrons and positrons annihilate to photons, the number of photons 
and hence also the temperature of the photon gas rises (this is called "reheating"). 
Entropy conservation yields that the number of photons, and hence the photon energy- 
and number density, rise by a factor 11/4. This is also the reason for today's difference 
in photon and neutrino temperature, r° = 2.73 K, = 1.95 K (r°/T° = (11/4)1/^). 

6.4.2 Initial conditions 

We start our simulation of the BBN process at a temperature of T = 10^^ K. As has 
been explained in Sec. 16. 31 protons and neutrons are in thermodynamic equilibrium via 
the reaction (|6.ip . hence the neutron-to-proton ratio is given by the thermodynamic 
relation 

^ _ g-(m„-mp)c^/fcBT ^ (6.11) 
P 

Knowing the starting temperature (which has to be well above the freezeout temper- 
ature of the reaction (|6.ip ). the initial ratio of neutrons to protons is well defined. 
Via Eq. (j6.6p the initial temperature can be translated into the initial time tinu- 

The only additional cosmological parameter necessary for BBN is the initial baryon 
number densitjH, 

initial initial i initial (a -i 

— ^n + ■ lO.lZj 

In principle, this parameter should drop out of a theory of baryogenesis. However, 
lacking this theory, we have to take ns as a cosmological parameter which has to be 
plugged into the BBN simulation by hand. Since ns scales like the photon number 
density during cosmological expansion, 

ub (X rij oc , (6.13) 

the parameter 

r] := UB/n^ = 2.74 • lO^^f^b/i^ (6.14) 



further cosmological parameter often considered in BBN studies is the 'number of neutrino 
species' at BBN, i.e. the number of relativistic particles present at BBN. Since the presence of those 
additional particles is equivalent to a change in expansion rate, we treat this parameter effectively 
as a change in the gravitational constant Gn- In SBBN, rj is the only parameter which enters, and 
Ni, is set to 3. 
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stays constant throughout the later evolution of the Universe, unless the number 
density of the photons is not significantly increased by later events, e.g. by the decay 
of particles. This is exactly what happens during e~^e~ annhililation as described in 
section [5.4.11 The effect on the photon temperature could be computed, giving that 
?7 at T > 10^° K was higher than the current r]o by a factor ^, ?7T>ioioi<' = x'^o- 
Including this correction factor, we can directly deduce the baryon density at BBN 
(for a given temperature T « 10^ K) from today's value of rj by applying Eq. (|6.5p . 

ssj^ 11 2C(3) /fcsTy 

--Vo-^[-f-) ■ (6.15) 

As this quantity determines the frequency of nuclear collisions, it is clear that a 
modified ub will change the "speed" of the BBN process. 

6.4.3 The element synthesis process 

During BBN, element production is dominantly driven by three types of processed 

• 2^1 fusion processes A + B ^ C , 

• 1^2 scatterings ^ + B ^ C + Z?, and 

• particle decays, A ^ B + C . 

The parameters describing these processes are the reaction cross sections ua+b—>c-, 
'^A+ B^c+D and the decay width \a^b+c = (inverse mean lifetime). Given these 
quantities, one derives simple differential equations that describe the time evolution 
of the number densities of the nuclei A, B, C, D, ua, ns, nc, 



For a decay A ^ B + C, 

duA 
"dT 

duB 
dt 

dnc 
"df 

For a 2 — + 1 process A + B ^ C 

duA 



-XuA (6.16) 
-XuA (6.17) 
-XuA (6.18) 



dt 
driB 

dt 
dnc 

dt 



-n^ns(o'u) (6.19) 
-nAnB(crv) (6.20) 
-nAUBicrv) (6.21) 



*Due to the comparably low temperature and density at BBN, three-particle interactions which 
are important in the stellar nucleosynthesis process do not play a role at BBN [AprahamianOS, and 
are hence not considered in BBN simulations. 
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• For a 2 2 process A + B ^ C + D 

duA 



dt 
duB 
IT 
dnc 
"dT 
dnu 



dt 

Here {av) is the Maxwell-Boltzmann averaged cross section. 



—nAnsicrv) (6.22) 

—nAnB{(yv) (6.23) 

+nAnB{ov) (6.24) 

+nAnB{(Tv) . (6.25) 



fi^/^{kBT)y- 



poo 

■ / a{E)Eexp(~E/kBT)dE 
Jq 



(6.26) 



where fi is the reduced mass of the reactants, and E is the reaction energy in the 
center of mass system. 

6.5 Nuclear reaction rates and the Q value 

For nuclear reactions, one of the most important quantities is the reaction Q value, 
which gives the amount of energy released (or absorbed) in a reaction. Working with 
positive binding energie^ Bi , the Q value for a reaction is defined as 

outgoing incoming 

Hence each reaction Q-value is determined by the masses of reactants and products. 
The Q-value of each reaction mainly affects the abundances via the reverse thermal 
reaction rate relative to the forward rate, and via phase space and radiative emission 
factors in the reaction cross-sections. 

The reverse reaction rate is simply related to the forward rate via statistical factors, 
due to time reversal invariance (see for example [NACRE99] ): the relevant dependence 
for a 12 ^ 34 reaction is 

« . (6.28) 

(av) 12^34: 

The Q-dependence of radiative capture reactions A+B C+j, assuming a dominant 
electric dipole, is 

aiE) <x E^ ^ (Q + E)\ (6.29) 
whereas for 2^2 inelastic scattering or transfer reactions the dependence is 

cr{E) (X (3 ^ {Q + E)^^\ (6.30) 

where f3 is the outgoing channel velocity. In the current treatment we assume E <^ Q 
at relevant temperatures, and simply scale rates by the appropriate power of P'^ 



^The positive binding energy is defined via B{A, Z) := Zirip + (A — Z)mn -— m{A, Z). 

® Clearly this breaks down when Q approaches zero, and we have not considered varying any 
binding energy to the point where this happens. A more accurate treatment would involve applying 
the phase space dependences directly to the cross-sections, for example in the S-factor description 
(see Sec. I6.6.3|l of charged particle reactions, which involves an expansion in E\ the dependence on 
(Q + E) can then be applied order by order. 
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Reaction 


Q value [MeV] 


D 


^He 


-iHe 




'Li 


n <-!■ p 


1.29 


-0.80 


-0.30 


-1.48 


-2.76 


-0.92 


p{n, 7)(i 


2.22 


-0.2 


0.1 





-0.2 


1.3 


7)^He 


5.49 


-0.3 


0.4 





-0.3 


0.6 


did, n)^Hc 


3.27 


-0.5 


0.2 





-0.5 


0.7 


d{d, p)t 


4.03 


-0.5 


-0.3 





-0.5 


0.1 


d{a,jf Li 


1.47 











1.0 





'^llc{n,p)t 


0.76 





-0.2 








-0.3 




18.35 





-0.8 








-0.7 


^He(a,7)^Be 


1.59 














1.0 


''Li(]5, a)^He 


4.02 











-1.0 





^Be(?i,]5)^Li 


1.64 














-0.7 


^Be(n,a)4Hc 


18.99 














-0.01 


T(p,g)^He 


19.81 














0.02 


'^Li(p,a)*He 


17.35 














-0.06 


T(a,g)'^Li 


2.47 














0.03 


T(d,n)4Hc 


17.59 





-0.01 








-0.02 


'^Be(d,pa)4Hc 


16.77 














-0.01 



Table 6.1: Leading dependence of abundances on thermal averaged cross-sections 
dhiYa/dhi{av)i for important reactions (1st part) and less important reactions (2nd 
part) 

A variation in binding energies can have two kinds of effect. The Q value can 
change the time when a reaction drops out of equilibrium, for instance at the n + p ^ 
d+j reaction. Or it can change the absolute rate of a reaction, and thus the production 
rate of a given species, for example at the ^Be-producing reaction whose cross-section 
varies with Q^. 

Whether the reaction matrix elements have a dependence on the binding energies 
and on Q is in general not clear because there is no systematic effective theory for 
multi-nucleon reactions. The exception is the np — s- d'j reaction, for which we have 
implemented a nuclear effective theory result (see Sec. I6.6.2|) . 

6.6 Nuclear reactions important for BBN 

Simulations of the BBN process usually track a large amount of elements and reac- 
tions. However, only a few of them turn out to be important for BBN. In the next 
section, we will be interested in the behavior of BBN under variations of parameters, 
hence we will be concentrating on those reactions where a variation of parameters 
will result in a variation in final abundances. In order to estimate which reactions are 
more or less important for the variation of the final abundances, we varied each ther- 
mal averaged cross-section (av) by a temperature-independent factor, preserving the 
relation between forward and reverse rates. The resulting values for dlnYa/d\n{av)i 
are given in Tab. 16.11 As usual in BBN simulations, the slow /3-decays of tritium 
and ^Be are accounted for by adding on the T and ^Be abundances to '^He and ^Li 
respectively at the end of the run, when other nuclear reactions have frozen out. 

For all reactions in Tab. l6.11 the comments on Q-value (and hence binding energy) 
dependence of the preceding sections do apply. In the following sections we will 
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give details on those reactions for which a more detailed dependence on fundamental 
parameters is known. 



6.6.1 The n p reaction rate 

The n ^ p weak interactions influence every abundance nontrivially. Opposed to 
most of the reaction rates which can only be determined experimentally, a closed 
formula is known for the n ^ p weak reaction jScherrer"83l [LopezQ?] |WeinbergGRT| . 

— >p ^e.7i — ^i^p ^ ^i^n — >ep ^ — >pei^ (6.31) 
— '71 — Agp — ^jyji -\- \i/p — 'en ^peu — *n (6.32) 



oo 



[H 


- exp(-ez)] [1 + exp((e - q)zi, - 






(e-g)2(e2_i)i/2e 




'[1- 


f exp(— ez)][l + exp((e — q)z^ 






(e + g)2(e2_i)i/2e 




[H 


- cxp(ez)][l + cxp(— (e + q)zi, - 


-Ce]) 



= K / de- ' ^' ^ ^ — (6.34) 



Xne^p.-Kj^ ^'[l+CXp(ez)][l + CXp(-(e + (7)z,-ee]) ' ^^'^^^ 

Here q := (m„ — mp)/me, z = me/T^, Zi, = mejTy and is the electron neutrino 
degeneracy parameter which we always set to zero. The constant K is obtained by 
demanding Xn^p(T — > 0) = t^^. The corresponding p n rates are obtained via 

Ap— = Xn^pi^ q^ ^e)- 

The equations given above derive from first-order electroweak theory and contain 
some approximations to allow an easy numerical evaluation. However, higher preci- 
sion of measurements demands a higher accuracy for theoretical simulations. Hence, 
modern simulations of BBN apply corrections to the equations given above. Those 
are radiative corrections in the form of zero-temperature and thermal Coulomb cor- 
rections [LopezQS] , as well as finite nucleon mass corrections given in |Lopez97] . For 
example, thermal Coulomb corrections are applied by multiplying the integrand of 
the rates for n pev and ep <-> nv with the Fermi factor 

^'^"/^ (6.36) 



1 — exp(— 27ra//3) 



where /3 = \J\ — e 2. Similar prescriptions exist for the other corrections [LopezQ?) 
|Lopez98| . 
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6.6.2 The n + p — > D + 7 reaction rate 

Chen et al. derived an effective theory for the strong 71 + p ^ Z? + 7 cross section 
which is now widely used [Chen99| . They use an effective theory which contains direct 
nucleon-nucleon interactions and photon exchange, but neglects pion exchange^. They 
estimate the error in the energy regime of BBN to be < 4%. Using nucleon-nucleon 
phase-shift data and the cross section for cold neutron capture as input data, they 
obtain the foUoing equations 



(T. 



np — 'd^ 



47ra(7^ -I- P 
j^m%P 



2\3 



X 



Ml 



X 



El 



(6.37) 



where P — \fmjqE^ 7 — ^/BjymN, is the deuteron binding energy and E is the 
cross section energy in the center of mass system. X^^^ and X^^ are given by 



~2 ^_P^m^ 

(72+P2)4 



1 + 7Pd + (7/Od)' + ilPd? + 



Gtt V 3 



Ci 



(6.38) 



«57'(i-7)- 

^Ml — 



1 +p2 (^2+p2)^ 



1 + 7Pd - T-o 



\ ai 



p2 p 



with the EFT fitting constants [ChenQSj IChen99| 



LnprriN 7^ + 

P2)(J._^) 2nn, J^-7_ 

(6.39) 



Pd 


1.764 fm 


(effective range parameter) 


Ci = 


-1.49 fm'' 


(P-wave interaction constant) 


Kl = 


2.352945 


(isovector nucleon magnetic moment) 


ai = 


-23.714 fm 


(scattering length in the "'S'o channel) 


ro = 


2.73 fm 


(effective range in the ^ 5*0 channel) 




-4.513 fm^ 


{Ml capture constant) . 



6.6.3 Charged particle reaction rates 

Theoretical models describing the rough shape of charged particle reaction cross sec- 
tions have long been available ^Fowler67j. In absence of resonances, the cross sections 
are a product of the Gamow factor and an "S-factor" , 

aiE) = S{Ef-^ , (6.40) 

where 

fl = aZiZ2^flIj2E , (6.41) 

Zi^2 the atomic numbers in the initial state and /i is the reduced mass. The S-factor 
may be expanded in a Maclaurin scries to quadratic order in energy, which is usually 
sufficient to account for any smoothly-varying dependence, 

S{E) = 5*0 + SiE + 52^2 . (6.42) 

'^In EFT treatments of processes where the external momenta are much smaller than the pion 
mass (as is the case at BBN with nucleon energies Ejv ^ 1 MeV, rrivr ~ 135 MeV), a pion- less EFT 
turns out to describe the process sufficiently well |Gegelia98| |ChenRupak99| . 
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To obtain the Maxwell-Boltzmann averaged cross section (Eq. (|6.26p ) one has to 
perform an integration over energy, resulting in a sum of terms for (av) with specific 
dependence on a and the reduced mass /i. A generic cross section has the form 
|Bergstrom99| !| 



[av 



a 



ao 



36 



Mo 



-K 3 



,1/3 

a 
ao 



Mo 



^ 89 1 

1 H K 3 

36 



35 
36' 



-K 3 



[resonance terms], (6.43) 



where 



2 „ 2 r^l ryl 



M 



2kBT 



(6.44) 



and the are fitting constants which are fit to the measured cross-sections. Some 
cross-sections are fit with an additional exponential term S'(0)e~^^ jFowler75] . In 
addition, non-resonant terms may be multiplied by a cutoff factor /cut = 6"^"^/"^'="*^ , 
where Tcut has been argued to be proportional to [Bergstrom99[ IFowler75j . 



Resonances 

Where the cross-section as a function of energy shows one or more resonances, they 
contribute to the thermal averaged rate as 

{c7v),,s = g{T)e-^/^ , (6.45) 

where g{T) and E are fitting parameters corresponding to the shape and position of 
the resonance. Usually a power-law is taken for g{T), thus g{T) = cT^. In principle 
one should consider the variation of the resonance parameters if this term is significant. 
But since the major contribution to the resonance energy E probably arises from 
^QCD which we take as our (non- varying) unit, it seems a reasonable first guess to 
keep the resonance parameters fixed. 

For the code, the NACRE formulae |NACRE99| fitted at the level of the thermal 
averaged cross-sections {av) are not suitable as they do not allow to incorporate 
the dependence on physical parameters. We use the functional forms of rates from 
|Smith92l |Bergstrom99| which have been described above in Eq. (|6.43p and fit the 
free parameters of these rates to reproduce the NETGEN rates |NETGEN] as closely 
as possible. We also checked that the resulting cross-sections are consistent with 
experimental data. In the case of c?(a, 7)®Li we found a set of parameters which 
seems to fit the experimental cross-section at low energies |Kiener6Li] better than the 
NACRE fit. But note that this cross-section is not measured directly, rather it is 
derived from experimental data under various assumptions. 



|Bergstrom99| only gives the dependence on a, we have recomputed the expansion including the 
dependence on the reduced mass fi. 
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Replacing the NETGEN rates in the code with our fitted rates, the obtained 
abundances change as follows: 

• Yu differs by -0.3%, 

• Yane by +0.9%, 

• Y4B_c by less than 0.1%, and 

• Yju by +3%. 

Hence we do not consider this refitting as significant, except in the case of the 
d(Q;,7)^Li reaction. Depending on whether this reaction is fit to NETGEN, or to 
the cross-section values of [Kiener6Li] . we found a ^Li abundance larger by a factor 
of 1.02, or 3.3, respectively. However, given the unclear observational status of ^Li 
this discrepancy is not currently worth pursuing (see Sec. I6.8."5|) . 



6.7 The simulation of the BBN process 

The BBN process has been described and simulated more than 40 years ago 'Wa goner66| . 
One starts with the initial conditions given in Sec. 16.4.21 and evolves the element syn- 
thesis processes for a set of nuclei as described by the differential equations in Sec. 
16.4.31 in the background of an expanding universe. As we are working at a time 
when the Universe was radiation-dominated, the expansion of the Universe as given 
by the Friedmann equation p.4p is dominated by the photon, neutrino and electron 
contribution^ Hence the energy density p entering in Eq. p.4p is 

p = + + (6.46) 

and the corresponding pressure 

1 2 1 2 . ^ 

P^-^P-yC + ^PyC +Pe- (6.47) 

Whilst photons and neutrinos are ultra relativistic throughout BBN, the electrons 
move from the relativistic to the non relativistic domain in the course of the BBN 
simulation (10 MeV }^ T > 0.001 MeV). Hence the electron density and pressure 
have to be described very carefully. 

For the electron energy and pressure density one introduces the electron chemical 
potential 

■K'^{hcfnBYp 
'^^^ 2{kBTzY 

2 

where z = . The electron and positron energy and pressure densities are then 
|FowlerHoyle64| 

p, ^ p,- + = -^r^ Y.^-T^' cosh(n</>e)M(nz) (6.49) 

^=$g?f:^ («») 



V„(-)»+ini(nz)_ 



(6.48) 



^The computer code also includes the evolution of baryonic and dark matter. However, this does 
not make any impact on the obtained abundances. 
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The functions L and M are related to the modified Bessel functions Ki via 

L{z):^K2{z)/z (6.51) 

(6.52) 



M(z) := - 

z 



^-K^(^z) + -K^{z) 



In our BBN simulation, we evolve the following set of parameters with time 

• temperatures T-y, T^, 

• electron chemical potential 0e 

• baryon density pB (or equivalently the scale factor a) 

• the abundances Yi 

according to the differential equations given in this chapter, and the Friedmann equa- 
tion (j3.4p is used to evolve the baryon density and temperatures according to the 
expansion of the Universe. 

6.7.1 Numerical aspects of the BBN simulation 

The main task of a BBN code is to numerically integrate the set of coupled differential 
equations which have been given in the preceding sections. In our code, the differential 
equations are solved using an adaptive second-order Runge-Kutta method, which 
is important to correctly account for the nuclear reactiona^^l. Given a differential 
equation 

V\t)^i(t,y) (6.53) 

and a starting value at time ioi 2/(^o) = 2/0: the task is to derive the value of y at time 
to + At. The second order Runge-Kutta method does this in a two-step approach. 
First one derives the value in linear approximation, i.e. 

y{to + At) = y{to) + At ■ y'{to) . (6.54) 

The "average" derivative between to and to + At is then defined as 

~, fito,yo) + f{to + At,yito + At)) 

y := (6.55) 

and the final value 

yito + At) =yo+y'-At. (6.56) 

The step width At is determined adaptively by demanding that the changes in abun- 
dances per time step do not exceed a certain value. 

We have adapted the code to the capabilities of modern computer technologies, 
which means that we have increased the internal numerical precision, implemented 
more precise integration routines and added several further numerical improvements. 
These corrections were necessary to remove numerical inaccuracies present in the 
available BBN codes. 



^"Some nuclear reactions can have both high forward and backward reaction rates. A naive solution 
of the differential equations for the nuclear reactions via y{to + At) = y{to) + At ■ J/' (to) would not 
take into account that produced nuclei might be immediately destroyed after creation via the back 
reaction process. 
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In general, the uncertainties of the input parameters for BBN resuh in much 
higher uncertainties of the resulting abundances than do the numerical inaccuracies 
of the code. However, in the next chapter we will derive linearized dependences of 
the final abundances Y on input parameters X, d\nY/d\nX. This will be achieved 
by slightly varying the input parameters X and observing the resulting changes in 
the abundances Y, where the variations in X are much below their uncertainty in 
order to get the correct linear order approximation. Hence, we have to ensure that 
the changes in the final results of the BBN simulations under small changes of the 
input parameters are not affected by numerical effects but solely by changes in the 
physical processes. 

6.8 Observational situation and uncertainties 

One of the biggest successes of standard BBN is the matching of theoretically pre- 
dicted and observed primordial abundances for major elements. For a review of the 
theoretical and observational status and obstacles see for instance |Steigman05| . 

6.8.1 ^He 

The highest precision observation is that of the ^He abundance (conventionally writ- 
ten Yp). The post-BBN evolution of ^He is simple. In stars, hydrogen is burned to 
^He which increases the abundance of ''He above its primordial value. Hence one 
expects the ''He abundance to decrease once one goes to stars with lower metallic- 
itjo, and a ''He "plateau" is expected for sufficiently low metallicity. Olive et al. 
|OliveSkillman04] analyze 8 systems with very low oxygen abundance, which are dis- 
played in Fig. 16.31 They argue that the observational data indicated a primordial 
abundance of 

Fp = 0.249 ± 0.009 , (6.57) 

which we take to be a la range. However, given the probable dominance of sys- 
tematic effects, instead of using 2-a bounds to later determine the range of al- 
lowed variationlil, we rather use the "conservative allowable range" of Yp given in 
|OliveSkillman04| as 

0.232 <Yp< 0.258. (6.58) 

6.8.2 Deuterium 

For deuterium, the post-BBN evolution is as clear as it is for "^He. Deuterium which 
was formed during primordial nucleosynthesis is burned in stars to ^He and higher 
elements. Also, any deuterium which is newly produced in stars is immediately burned 
into ■^He as deuterium is the most weakly bound light nucleus. Hence, the deuterium 
abundance should increase when going back in time (opposed to ^He, which decreases). 
However, the almost identical absorption spectra of D and "normal" hydrogen (^H) 
complicate spectroscopic determinations of deuterium significantly. 

Elements with higher mass number than Helium (which arc in astronomy called "metals") are 
only produced in stars. Hence a low metallicity of a system is a strong indication of element abun- 
dances which are assumed to be close to the primordial ones. 
12 See Sec. [911 
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10i'(0/H) 



Figure 6.3: ^He abundances versus oxygen abundance from SteigmanOS . The solid 
line is the weighted mean of for the 8 systems. Data from |OliveSkillman04j 



Due to this complication, the determination of the primordial deuterium abun- 
dance follows from a small number of observed systems. One basically observes ab- 
sorption spectra of interstellar deuterium from sightlines of Quasars which have to 
fulfill certain quality conditions (see ^Kirkman03j for details). As D/H is of the or- 
der 10^^, the hydrogen column density Nhi must be large enough in order to be 
able to observe deuterium with modern high-resolution spectrograph^. We use the 
determinations of [O'MearaOSi IKirkmanOS] which give a value of 



D/H = (2.8 ±0.4) X 10" 



(6.59) 



The analysis of [O'Meara06| is displayed in Fig. 16.41 where the large scatter between 
values determined from different systems should be noted. After our studies have 
been performed, a new analysis appeared [PettiniOS] which reduces the error of D/H 
by a factor of 2. 



6.8.3 ^He 

The post-BBN development of "^He is quite complex, as it is both produced and de- 
stroyed in stars. As quantitative analyses of these competing processes are quite model 
dependent, the primordial "^He abundance cannot be extracted easily from spectra of 
old stars. This complexity is revealed by determinations of the '^He abundance, which 
typically have a large scatter as can be seen in Fig. 16.51 Hence the common under- 
standing is that '^He abundance determinations cannot be considered as good tracers 
for the primordial abundance |Vangioni-Flam02] . 



^^The column density Nhi gives a two-dimensional measure of the density (here number density) 
of a cloud of material, measured in cm~^. 
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Figure 6.4: Deuterium versus neutral column density Nhi for a set of low-metallicity 
absorption spectra along QSO sightlines. From [O'Meara06| . {Nhi in cm~^, its value 
is of no relevance here). The horizontal line represents the weighted mean, the right 
axis gives the derived SBBN values for flth"^. 




10^(0/H) 



Figure 6.5: "^He abundances versus oxygen from Steigman 05]. Data from |Bama02j . 
The dashed lines show the la band given by |Bania02] . the blue spot indicates the 
■^He abundance for the pre-solar nebula. 
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Figure 6.6: ^Li abundances versus metallicity from [Steigman05| . loge(Li) 
log(Li/H) 
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6.8.4 ^Li 

The post-BBN ^Li abundance increases due to cosmic ray spallation/fusion reactions 
and ''Li production in stars. Even though ^Li is also easily destroyed inside stars, some 
stars observationally appear to be "super- lithium rich" Stei gman05| , supporting the 
assumption that (most) stars are in fact net producers of lithium. Hence, one expects 
the lithium abundance to decrease when going back in time, which is supported by 
observations of the ''Li abundance which show a plateau at old, metal-poor halo stars 
(see Fig. 16. 6|) . One assumes that the plateau value is closely related to the primordial 
one. 

The most recent determinations of the abundance are quite small: ^Li/H = (1.3 ± 
0.3) • 10~^° [Bonifacio06 (see also |Asplund05 ) . It has been suggested that there 
are unresolved systematic errors relating to the effective temperature of the stars 
|Charbonnel05| which may imply a value as large as ^Li/H = (1.64 ± 0.3) • 10~^°. To 
account for this possible systematic, we adopt a value 



^Li/H= (1.5 ±0.5) X 10" 



(6.60) 



It turns out that the observed ^Li abundance is about a factor of 3 smaller than the 
standard theoretical prediction, which is called the "lithium problem" . 
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Abundance 


Observational 


Theoretical 


D/H 


(2.8 ±0.4) X 10-^ 


(2.61 ±0.04) X 10-^ 


Yp 


0.249 ±0.009 


0.2478 ±0.0002 


'U/ H 


(1.5 ±0.5) X 10^1" 


(4.5 ±0.4) X 10-1" 



Table 6.2: Current observational and theoretical primordial abundances 



6.8.5 ^Li 

A possible detection of ^Li was discussed in |Asplund05| , who claim to have found 
signals of ^Li in nine stars at the > 2a significance level. Their observation suggests 
a ^Li plateau at the level of 

6Li/H« 6.2 X 10-12 (g gi) 

If this detection is correct, the ^Li abundance is about a factor 100 larger than the 
SBBN prediction. However, given the unclear observational status and post-BBN 
history of this isotope, we do not include ®Li in the final analysis. 



6.8.6 Theoretical predictions 



Theoretically predicted primordial abundances also come with an error, mainly due to 
cross-section uncertainties. Our numerical procedures do not provide error estimates, 
so we adopt the la ranges from SerpicoOl] , using a baryon density ^Ibh"^ = 0.0224 
|ManganoTalk07[ IHjns hawWMAPS] : 

D/H = (2.61 ± 0.04) X 10"^ 

3He/H = (1.03 ± 0.03) X 10"^ 

Yp = 0.2478 ± 0.0002 

^Li/H = (4.5 ± 0.4) X 10^1°. (6.62) 



A compilation of the recent observational and theoretical primordial abundances 
which we will use for BBN is given in Tab. 16.21 



Chapter 7 



BBN with varying constants 



7.1 Nuclear and fundamental parameters 

Big Bang Nucleosynthesis happens at energies in the keV regime, much below the 
energy scale of the Standard Model of particle physics where the parameters which 
are nowadays considered as fundamental (Tab. 12. ip are relevant. In particular, the 
theory of QCD cannot be applied at BBN as it happens well below the QCD scale 
^QCD ~ 200 MeV. Quarks have combined into nucleons and the strong coupling 
constant as is in the non-perturbative regime. Hence, rather than the parameters of 
the SM, a set of effective low-energy parameters (masses of nucleons, binding energies, 
cross-sections and lifetimes) enters in the BBN simulations. We call these parameters 
"nuclear parameters" which stresses that those are the actual relevant (effective) 
parameters which enter in the nuclear physics processes. 

7.2 Nuclear parameters relevant for BBN 

The following set of effective low-energy ( "nuclear" ) parameters enters the BBN sim- 
ulations: 

• Neutron and proton mass m„, mp, combined to 

• Average nucleon mass tojv := (f7?,„ -I- mp)/2 

• Neutron proton mass difference Qjq :— m„ — nip 

• Neutron lifetime t„ 

• Binding energies of D, T, ^He, ^He, ^Li, '^Li, ^Be. 

Additionally, several fundamental parameters are also used in the BBN simulation: 

• Electron mass 

• Gravitational constant Gn 

• Fine structure constant a. 
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The baryon to photon ratio 77 enters as a further important parameter which is of 
cosmological origin. If one had a closed picture for the origin of our Universe, this ratio 
would derive from "fundamental" parameters, but lacking any theory this value only 
enters as an observed cosmological quantity. We combine all parameters which enter 
in the BBN simulation into a set Xi which we will call "set of nuclear parameters" 

7.3 Nuclear parameter dependence and the response 
matrix 

We consider the set of primordial abundances Ya with a = (D, ^He, "^He, ^Li, ^Li) 
and study its dependence on the variation of our set of nuclear physics parameters 
Xi. Here the index i denotes the parameters which enter the calculation of nuclear 
abundances, listed in the preceding section. Our central quantity is the response 
matrix C with matrix elements |MSW04] 

It indicates the leading linear dependence for small deviations of the abundances about 
the values obtained given the nuclear parameters inferred from present laboratory 
experiments. The matrix C is extracted by varying the quantities Xi independently 
in the BBN code, a procedure which includes variation of the reaction cross-sections 
and rates that have a physical dependence on Xi. If all variations in parameters are 
taken to be small, all necessary information can indeed be extracted from the response 
matrix. 

Our results for the nuclear response matrix are shown in Tab. 17.11 The first 
thirteen rows constitute the transposed nuclear response matrix C"^. We also quote 
the dependence of the abundances on 77 in the last row. Values are quoted to 2 d. p. 
or to 2 sig. fig. when the magnitude exceeds 1. 



^Note that for instance variations with respect to the nuclear parameter a are different from 
variations with respect to the fundamental parameter a. This is due to the fact that variations 
w.r.t. the fundamental parameter incorporate all dependences, whereas variations w.r.t. the nuclear 
parameter only include some dependences. For instance the dependence of t„ on a will only be 
included in the fundamental parameter dependence. 
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dlnYa/dlnXi 


D 


^He 


^He 








0.41 


0.15 


0.73 


1.4 


0.43 


Qn 


0.83 


0.31 


1.55 


2.9 


1.00 


TUN 


3.5 


0.11 


-0.07 


2.0 


-12 


Bb 


-2.8 


-2.1 


0.68 


-6.8 


8.8 


B'Y 


-0.22 


-1.4 





-0.20 


-2.5 


^3He 


-2.1 


3.0 





-3.1 


-9.5 


-B4He 


-0.01 


-0.57 





-59 


-57 


-BeLi 











69 





BtLi 














-6.9 


BjBe 














81 


Gn 


0.94 


0.33 


0.36 


1.4 


-0.72 


a 


2.3 


0.79 


0.00 


4.6 


-8.1 


rUe 


-0.16 


-0.02 


-0.71 


-1.1 


-0.82 


V 


-1.6 


-0.57 


0.04 


-1.5 


2.1 



Table 7.1: Response matrix C, dependence of abundances on nuclear parameters. 



Chapter 8 

From nuclear to fundamental 
parameters 



8.1 From nuclear to fundamental parameters 

Looking at the Standard Model of particle physics, it is clear that the nuclear pa- 
rameters given in the preceding section (Tab. 17. ip are degenerate in a sense that the 
set of 10 non-fundamental parameters only depends on about 5 to 6 fundamental 
parameters. Hence, in a next step, we will derive relations between a set of Standard 
Model parameters Gk and the nuclear physics parameters Xi. This is encoded in a 
second response matrix F with entries 

The variation of abundances with respect to the fundamental parameters Gk is then 
expressed by the "fundamental response matrix" R with elements raki 

^Ya AGk ,^ 

-r^=r„fe-— . (8.2) 

The matrix R is obtained from C and F by simple matrix multiplication, 

R^CF. (8.3) 
We consider the following six fundamental parameters Gk- 

• Gravitational constant Gn 

• Fine structure constant a 

• Electron mass rrie 

• Light quark mass difference 6q = rud — m„ 



• Averaged light quark mass m = (m^ -t- r7i„)/2 oc m 

• Higgs v.e.v. (0). 



2 
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An additional parameter is the strange quark mass m^. We have omitted it from 
our treatment of BBN, as it enters the parameters (binding energies, neutron and 
proton mass) with high uncertainty. We have computed that the final dependence 
of abundances on the strange quark mass due to the known proton and neutron 
dependence is less than 3% of that of the light quark mass dependence (given in 
Tab. 18. 2p . hence the dependence on TOs is much lower than the model uncertaint}0 
{e.g. for nuclear binding energies). However, in our study of variations in Chap. [TT] we 
will include the strange quark mass contributions, as the parameters studied there are 
more directly influenced by strange effects. This is why we will also include the strange 
contribution in our following derivation of fundamental parameter dependences. 

In the next sections, we give details on how the nuclear parameters depend on 
our set of fundamental parameters. In effective theories of nuclear forces, the pion 
appears as an effective mediator of the strong force. Hence, we will also introduce the 
pion mass to^ as an intermediate parameter. 



8.1.1 Pion mass 

The pion as the lightest of all mesons is the dominant mediator of the strong inter- 
action in effective low energy theories of QCD. Its light mass is due to the fact that 
the pion is the pseudo Goldstone boson of the only slightly broken chiral symmetry of 
QCD. This fact implies that one can make several predictions about pion properties 
in chiral perturbation theory. In first order chiral perturbation theory, one obtains 
the famous Gell-Mann-Oakes-Renncr relation [Gell-Mann68j . which relates the pion 
mass to fundamental parameter j^ Gasser82] . 

~ -^(to„ + mrf)(0|uu + Jd|0) . (8.4) 

The nonvanishing v.e.v. of uu and dd is a measure of the chiral asymmetry of the 
vacuum, signaling the spontaneous breakdown of chiral symmetry |Gasser8"2) . As can 
be seen from Eq. (|8.4p . the values of {qq) are negative. is the decay constant of the 
pion, which does not depend on quark masses at first order and which has a value of 
|PD(;08j 

U ~ 130 MeV . (8.5) 

The formulas mentioned above are pure chiral perturbation theory equations 
neglecting any electromagnetic contributions. In fact, the pion mass difference of 
TO^± — TTiTrO = 4.6 MeV is dominantly due to electromagnetic effects |Gasser82] : up to 
a small uncertainty of 0.1 MeV the virtual photon cloud surrounding the Tr^accounts 
for the vr^-TT^mass difference, 

rrP± -rrPo =4.6 ±0.1 MeV 

with m^o « 0. We can safely neglect this contribution and only work with the first 
order dependence on the light quark mass, which can be read off from Eq. ()8.4p . 

(X rh (8.6) 

^However, it might be that binding energies get a substantial contribution from strange effects 
which we so far cannot quantify |FIambaum02| . 

^Note that there are in fact three pions, 77+, n~ and tt" with roughly the same mass {Tn^± = 
139.6 MeV, m^o = 135.0 MeV). We will only work with the first order dependences on fundamental 
parameters which are the same for all three pions. Further details on specific mass contributions can 
be found in IGasser8"2] . 
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or 

Alnm^ = -Alnm . (8.7) 
8.1.2 Neutron and proton mass 

The different contributions to tlie neutron and proton mass and the origin of the 
proton neutron mass difference have been studied by Gasser and Leutwyier more 
than 25 years ago |Gasser75l IGasser82) . Using "improved chiral perturbation theory", 
they estimate the contributions coming from eiectromagnetic and (effective) strong 
seif energy as weii as from different quark masses. Knowing the properties of eiectron 
proton scattering they could estimate the electromagnetic contributions to the masses 
at Born leveH from the electromagnetic form factors as 

ml, = 0.63 MeV 

(8.8) 

ml = -0.13 MeV 

to a very high accuracy. A further electromagnetic correction might come from pos- 
sible intermediate states, a contribution which could not be calculated but has been 
estimated to 

Amre, = ±0.2 MeV . (8.9) 

Thus, the total electromagnetic contributions to the nucleon masses are the values 
given in (|8.8p with an error of ±0.2 MeV. For the neutron-proton electromagnetic 
mass difference one obtains 

(m„ - mp)'< = -0.76 ± 0.30 MeV . (8.10) 

The bare QCD masses of the nucleons are hence 

m^^^ = 937.64 ± 0.20 MeV , 

(8 11) 

mf^^'^ = 939.70 ± 0.20 MeV . 

The difference in neutron and proton mass of 

(m„ - mp)^"^^ = 2.05 ± 0.30 MeV (8.12) 

can be explained in terms of the different quark content, p = (uud), n — (udd), which 
turns out to be in lowest order quark mass expansion |Gasser82] 

(mn—mp)'^'^^ = (md—mu)- < p\uu—dd\p >— 5„- < p\uu—dd\p > . (8.13) 

Equations (|8.10p . (|8.12p and (|8.13p can be combined into 

AQn ^ (-0.76Alna + 2.05A\nSq) MeV 

^ AIuQat ~ (-0.59Alna + 1.59Aln<5q) . (8.14) 

Subtracting the electromagnetic contributions, one is left with the pure QCD 
masses as given in Eqn (|8.1ip , which is far from only being the sum of the masses of 



^Born level contributions are proportional to a, so a rescaling of a by a certain factor scales the 
values given in Eq. 118.81 1 by the same factor. 
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the three valence quarks of the nucleons. A recent study using "heavy baryon chiral 
perturbation theory" |Borasoy96| yields that the baryon mass in the chiral limit is 

m= 767± 110 MeV. (8.15) 

This quantity is purely due to QCD effects which only depend on Aqcd, and since 

we are holding this quantity fixed, m remains constant. The contributions to the final 
nucleon mass come from the u, d, and s quarks, whose contributions can be estimated 
using the two quantities [Borasoy96] 

cr^^(O) := m < p\uu + dd\p 45 ± 10 MeV (8.16) 

y:=/^*^= 0.21 ±0.20. (8.17) 

< p\uu + aa\p > 

They are connected to the quark mass dependence of the nucleons via the Feynmann- 
Hellmann theorem [Feynman39[ IHellmann37| , 

a^^(0)=m^ (8.18) 
om 



and equivalently 

This can be translated into 

dlnniN 



9 In 



m 



<p\ss\p>^—^. (8.19) 
am.. 



= mjYV^Ar(O) = 0.048 ±0.011 (8.20) 



|^^ = ^f-]vV..(0) = 0.12 ±0.12 (8.21) 
a mm. m 2 



using ^ « 25 |PD(;08j . 



8.1.3 Neutron lifetime 



In the electroweak model the neutron lifetime can be evaluated analytically. One 
derives WeinbergGRT| (see also Eq. (|6.33p with zero temperature). 



r-' = A„_>p+,-+,^ = -^^Gf / {Qn - xf Jl - -^dx . (8.22) 

Using q :— ^f-, this integral evaluates to 
_ l + 3.9i , . 



y/q^ - 1 {2q^ - 9q^ - 8) + 15gln (g + y/q^-l^ 

.23) 

Here, Gp is the Fermi coupling constant, which is related to the Higgs v.e.v. via 
Eq. (I4.2ip and qa is the nucleon axial vector coupling constant, qa ~ 1.27 |PDG08| . 
Recent developments in lattice QCD, combined with chiral perturbation theory, allow 
to compute gA, yielding |Hemmert03| 

5^ = (1.2±0.1)-6.9(^^) ln^^-(10.4±4.8)m2 GeV-^ + 0(m3). (8.24) 
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Obviously the chiral value gA— 1.2 ± 0.1 dominates the physical value of qa, and for 
the derivative we obtain 

0.006 =^ ^= 0.003. (8.25) 
o m TO^ o In m 

Thus gA can be assumed constant in our treatment. Using the known dependences 
of Qjv from the preceding section, we arrive at a fundamental parameter dependence 
of the neutron lifetime of 

AlnT„ = 3.86Alna + 4Aln(0) + 1.52AlnTOe - 10.4Aln5g . (8.26) 
8.1.4 Binding energies 

The dependence of nuclear binding energies on the pion mass and a have been esti- 
mated in |Pudliner97] and [PieperOl] using quantum Monte Carlo calculations with 
realistic models of nuclear forces (similar values for the a dependence appear in 
|Nollett02] V We use the Pudliner and Pieper values which give for the dependence of 
binding energies on a: 

Aln(i?D, Bt, -Bshc, -B4HC, Bqu, Bju, Bjbc) — 

(-0.0081, -0.0047, -0.093, -0.030, -0.054, -0.046, -0.088) A Ina. (8.27) 

The pion mass determines the range of attractive nuclear forces, and the quantum 
Monte Carlo calculations of |Pudliner97| |PieperOl] which include pion exchange ac- 
curately reproduce many experimental properties. One-pion exchange and two-pion 
exchange are dominant contributions within the expectation values of the two- and 
three-nucleon potentials respectively. Currently such studies have not been extended 
to determine the functional dependence of binding energies on the pion mass in gen- 
eral |Flambaum07| IDamourOT] . This dependence would in any case have uncertainties 
due to subleading effects of pion mass (or equivalently light quark masses) on other 
terms in the nucleon-nucleon potential []3eaneSavage02j . 

However, the dependence of the deuteron binding energy on the pion mass has been 
extensively studied within low-energy effective theory |Epelbaum02| |BeaneSavage02| : 
the result may be expressed as 

AIuBd = ^-Alnm^ = ^Alnm (8.28) 

for small variations about the current value |YooScherrer02j . with —10 < r < —60 
We will also take this dependence as a guide for the likely pion mass dependence of 
other binding energies. Although the size of the deuteron binding appears due to 
an accidental cancellation between attractive and repulsive forces, its derivative with 
respect to (which is just Bd /m^ times r) is not expected to be subject to any 
cancellation. We also expect that the pion contribution to the total binding energy 
should increase with the number of nucleons; a proportionality to {A — 1) seems 
reasonable to obtain correct scaling at both small and large Hence to estimate 
the effect of pion mass on the binding energy of a nucleus Bi we set 

1^ = MA^ - 1)— r- ^ -0.13/.(A, - 1) , (8.29) 



^Our definition of r differs by a sign from I YooScherrer02 1 . 

^Flambaum et al. |Flambaum07) obtain rougfily a scaling of type (A—l) for nuclei with 3 < ^ ^ 7, 
but their results have a large uncertainty. 
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taking r ~ — 8. The numerical constants fi are expected to be of order unity, but 
will differ between light nuclei due to peculiarities of the shell structure, etc. Our 
normalization corresponds to /d — 1- Then the nontrivial dependences of nuclear 
parameters on rh are 

A \n{Bu , Bt , -BsHe , B4HC , Bqu , Bju , Bjbc , "Iat ) ~ 
(0.5r, 0.26/Tr, 0.29/3Hcr, 0.12/4Hcr, O.lT/gLir, 0.17 frur, O.lS/yBcr, 0.048)A Inm , 

(8.30) 

where the dependence of toat is taken from Eq. (|8.20p . For the m dependence of 
abundances due to the variation of binding energies we then have 



9 In J 



^ {A-l)BndlnY, 



where the dependence 1}^^° is obtained from the BBN code by varying the binding 
energies Bi. Taking account also of the small effect of rh on the nucleon mass m^v, 
the resulting dependence of abundances on rh is 



dlnrh 
d In YsHc 

9 In TO 
d In KtHc 

9 In TO 
dlnYeu 

dlnrh 
dlnYju 

dlnrh 



11 + 0.5/t + 5/3He 
8 + 3/t - 7/3He 

-2.7 

27 + 0.4/t + 7/3HC + 55/4HC - 96/6Li 

-36 + 5/t + 22/3HC + 54/4HC + 9/7Li - 115/7BC , (8.32) 



where we have neglected subleading terms. Even if we consider that some contribu- 
tions could cancel against one another due to the values of the fi, the magnitude of 
these variations is striking, particularly concerning the lithium abundances. To get 
an idea of the possible effect of cancellations, we may set all fi to unity and find the 
dependences 

A ln(yD, >3Hc, nnc, Yeu, >7Li) ^ (17, 5, -2.7, -6, -61) A Iuto. (8.33) 

One may also consider to what extent varying to or the pion mass may affect reaction 
cross-sections beyond the npd'j reaction. It seems very likely that matrix elements 
would acquire nontrivial dependence on to^; however, since the dependence of abun- 
dances on reaction cross-sections is relatively mild (see Table 16. 1|) , the dependence 
via reaction matrix elements is unlikely to compete with the very large effects arising 
through the variation of binding energies. 
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d\nXjd\nGk 


Gn 


a 


(</>) 






rh 


Gn 


1 

















a 





1 














T„ 





3.86 


4 


1.52 


-10.4 





me 











1 








Qn 





-0.59 








1.59 





rriN 

















0.048 


Bb 





-0.0081 











-4 







-0.0047 











-2.1,/t 







-0.093 











— 2.3/3H0 


S4HC 





-0.030 











-0.94/4HC 


-BeLi 





-0.054 











-1.4/6Li 







-0.046 











-1.4/7Li 


BjBc 





-0.088 











-1.4/7BC 



Table 8.1: Response matrix F, dependence of nuclear parameters Xi on fundamental 
parameters Gk 



dlnYa/dluGk 


D 


■^He 


^He 






Gn 


0.94 


0.33 


0.36 


1.4 


-0.72 


a 


3.6 


0.95 


1.9 


6.6 


-11 




1.6 


0.60 


2.9 


5.5 


1.7 


me 


0.46 


0.21 


0.40 


0.97 


-0.17 




-2.9 


-1.1 


-5.1 


-9.7 


-2.9 


rh 


17 


5.0 


-2.7 


-6 


-61 


V 


-1.6 


-0.57 


0.04 


-1.5 


2.1 



Table 8.2: Response matrix R, dependence of abundances Yi on fundamental param- 
eters Gfe. All fi are set to unity. 



8.2 The response matrices 

The results of the previous section can be summarized in the response matrix F 
defined in Eq. ()8.ip which relates variations in fundamental parameters Gk to vari- 
ations in nuclear parameters Xi. Its values are shown in Tab. 18.11 These results 
can be combined with the nuclear response matrix, Tab. 17. H to the "fundamental 
response matrix" R defined in Eq. (|8.2p according to Eq. (|8.3p . The matrix R relates 
variations in fundamental parameters to variations in final abundances and is shown 
in Tab. 18.21 This table is the central result of this part of the thesis. In treating the 
m-dependences, which arise from the nuclear binding energies with their uncertain 
values of fi, we have given the values which arise when setting all to unity. In 
our further treatment of BBN we will neglect the model uncertaintjO in the binding 
energies and always assume fi = 1. 



^Given the current status of low-energy QCD, it seems hard to quantify possible ranges of uncer- 
tainty for the fi. 
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dlnYa/dlnGk 


Gn 


a 


(0) 


nie 




m 


"He 


0.41 


1.94 


3.36 


0.389 


-5.358 


-1.59 



Table 8.3: Dependence of abundances Yi on fundamental parameters Gk found by 
|MSW04j 



ainy„/91nGfc 


D 


^Hc 


''He 


''Li 


'Li 


Gn 


0.4 


0.6 


0.3 


0.3 


-1.0 


a 


2.3 


1.0 


2.3 


7.4 


-9.5 


{4>) 


0.07 


0.1 


3.1 


4.1 


« 0.5 



Table 8.4: Dependence of abundances Yi on fundamental parameters from 
|Chamoun05| and [Landau04] 



8.3 Comparison to other studies 

The ^He dependence was previously calculated in |MSW04] by semi-analytic meth- 
ods [Esmailzadeh91| . Their findings are displayed in Tab. 18.31 our results for the 
dependence on fundamental parameters are similar. The dependence on Gn can 
be compared with the results of [Scherrer^ and [ChamounOSi ILandau04] , once one 
translates from units where Gn is constant to ours where Aqcd is constant. The 
latter ones also give values for the dependence on {(p) (variation of Gp c>c 
and on a. Their results are shown in Tab. 18.41 They roughly match with our esti- 
mates most times except for a few numbers, e.g. the much smaller variations of D 
and "^He under variations of the Higgs v.e.v. We assume that these deviations are 
due to the applied semi-analytical method of |Esmailzadeh91j which might not give 
appropriate results in the case of abundances other than "^He. Bergstrom and NoUett 
|Bergstrom99, Nollctt02j use the full BBN code [Kawano92] to constrain the varia- 
tion of a, incorporating a-dependent reaction cross-sections as described in Sec. 16.6.51 
From their graphical results we extract the dependences given in Tab. 18. 5) which again 
match with our findings. 



d\nYa/d\nGk 


D 


^Hc 


^He 




^Li 


a 


w 3.5 


w 0.9 


w 2.0 




w -7 



Table 8.5: Dependence of abundances Yi on a from |Bergstrom99| and [Nollett02] 



Chapter 9 



Constraints on variations of 
fundamental parameters 

9.1 Bounds on separate variations of fundamental 
couplings 

The first application of our results from the preceding chapter is in setting bounds 
on the variation of each fundamental parameter considered separately, under the 
assumption that only one parameter varies at once. We may consider three observa- 
tional determinations of primordial abundances (see Sec. 16. 8^ : deuterium, *He and 
""Li. However, the observed ''Li abundance deviates by a factor two to three from 
the value predicted by standard BBN theory (SBBN), and systematic uncertainties 
related to stellar evolution exist [Korn06| . Thus, we use the former two, D and *He, 
to constrain the allowed variations of the fundamental constants individually. For 
deuterium we take 2cr limits; for ^He we consider instead the "conservative allowable 
range " of [OliveSkillmaSoi] . The result ing constraints are given in Tab. 19.11 

9.2 Variations of abundances in unified models 

A major problem of using BBN to constrain variations of fundamental parameters is 
degeneracy. There are only three observational abundances (D, ^He and '^Li) but six 
fundamental parameters which can vary. Hence, the three observed values cannot be 
used to constrain a set of six parameters. In Chapter [3] we have introduced unified 



-19% 


< 


AIuGn 
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+10% 


-3.6% 


< 


A In a 


< 


+1.9% 


-2.3% 


< 


Aln(0) 


< 


+1.2% 


-17% 


< 


A In TOg 


< 


+9.0% 


-0.7% 


< 


A\nSq 


< 


+1.3% 


-1.3% 


< 


A In TO 


< 


+1.7% 



Table 9.1: Allowed individual variations (2(7 or "conservative allowable range", see 
Sec. 16. 8p of fundamental couplings 
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models where the variations of fundamental couplings satisfy relations that reduce 
the number of free parameters. Here we apply GUT models in order to be able to 
visualize our findings for BBN. In the simplest case every variation of a parameter Gk 
is determined by a single underlying degree of freedom, for instance a variation of the 
unified coupling A In ax ■ We may also eliminate A In ax in favor of some observable 
parameter, which we will choose to be the fine structure constant a. 

For grand unified theories, it makes sense to change from a system of units with 
constant Aqcd to a system where Mqut, the mass scale of the GUT, is set to 
constant. Furthermore, for the scenarios we will use in this section, we will take the 
Planck mass fixed relative to the unification scale, thus A{Mp/Mcut) = 0. We set 
all fi to unity and, for simplicity, we take the Yukawa couplings to be constan10, thus 
the electron and quark masses are proportional to {(f)), 

Aln-^ = Aln^ = Aln^=Aln-i^. (9.1) 

Mgut Mgut Mgut Mgut 

Finally we define an exponent 7 which relates the variation of (0) with respect to 
Mgut to the variation of Aqcd /Mgut as 



Mgut \ Mgut 

We study three non-supersymmetric GUT scenarios, 7 = 0, 7 = 1 and 7 = 1.5, 
and use the fine structure constant a to parametrize the variations. As has been 
shown in Sec. 14.61 the variation of a can be related to variations of particle masses, 
the GUT coupling ax and the QCD invariant scale Aqcd, see equation (|4.23p . 

Scenario 7 = 

In the first scenario, 7 = 0, the Higgs v.e.v. and elementary fermion masses are all 
proportional to the unification scale, A{Mgut/Mp, {(l))/Mp,me,q/A[p) — 0. Then 
the variations of fundamental couplings are related as 

A In = A In ^ = ^ A In a ^ 32.3 A In a . (9.3) 

Mgut {4) 40a ^ ' 

In a system with constant Aqcd we get 

A ln(GN, a, (<?!>) , m^, (5„ m) ~ (64.5, 1, -32.3, -32.3, -32.3, -32.3) A Ina. (9.4) 
We then obtain variations of abundances 

Aln(yD,1^3Hc,>4Hc,>6Li,l7Li) ^ (-450, -130, 170, 380, 1960)A ln« . (9.5) 

Scenario 7 = 1 

In the second scenario with 7 = 1 , all low-energy mass scales of particle physics are 
proportional to Aqcd- Setting Aqcd constant, we have 

A\n{GN,a,{<j)),me,Sq,m) ~ (78, 1, 0, 0, 0, 0)A Ina. (9.6) 

In this case the variations of abundances are not subject to the theoretical uncertainty 
of varying ttIs/ Aqcd- We obtain 

Aln(yD,l3Hc,i"4Hc,r6Li,l7Li) ^ (77, 27, 30, 120, -68) A Ina . (9.7) 



^See also the comments on variations of the Yukawa couplings in Sec. Ill.l1 
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Scenario 7 = 1.5 

In the third unified scenario with 7 = 1.5 we consider the case when the Higgs v.e.v. 
and fermion masses vary more rapidly (with respect to the unification scale) than the 
QCD scale Aqcd does. Converting back to a system where Aqcd is constant, we 
find that the variations of fundamental couphngs are related as 

Aln(GN,a, ((/)),me,(5q,TO) ~ (87, 1, 21.5, 21.5, 21.5, 21. 5)A Ina. (9.8) 

The variations of abundances are then 

Alii{Yu,Y3Hc,Y^nc,Yeu,Yju) ^ (430, 130, -65, -60, -1420) A Ina. (9.9) 

9.2.1 Linear results 

In Fig. 19.11 we show the abundance variations given by the three GUT models, as a 
function of the variation of a. First we plot only the linear dependence of abundances 
on a as given by the fundamental response matrix R (Tab. 18. 2p . We also show the 
Ict observational bounds as highlighted regions. Also included in the plot is the effect 
on the standard BBN predictions of varying the baryon-to-photon ratio 77 over the 2a 
range allowed by WMAP 3 year datfH, 5.7 < 10^°ri < 6.5. 

It can be seen that in the 7 = scenario a reduction of a by about 0.025% 
(i.e. a fractional variation of —2.5 x 10~^) would bring theory and observation into 
agreement within 2a bounds, while remaining in the linear regime. Conversely, in 
the 7 = 1.5 model an increase of a by about 0.04%, i.e. A In a = 4 x 10~*, brings 
theory and observation into agreement within la bounds. Considering the variations 
of fundamental parameters in the three scenarios Eqns. (|9.4p . (|9.6p and (|9.8p . the 
behavior of the weak scale {(jj) and fermion masses is decisive for the variation of 
abundances. However, if AlnYa becomes larger than 1 (as in the case of ^Li) the 
results are affected by higher order terms and the linear approximation can no longer 
be applied. 

9.2.2 Nonlinear results 

The linear analysis in the previous section suggests that it is possible, and may even 
be natural, to obtain a large negative variation in the ^Li abundance, and considerably 
smaller variations in other measurable abundances: positive in the case of deuterium 
and negative for "^He. Agreement between theory and data in all three abundances 
could then be possible for a narrow range of values in the variation of fundamental 
parameters. The required fractional variation in ^Li is so large (a factor two or more 
in IrLi) that a linear analysis using matrix multiplication is inaccurate. 

We improve the analysis by including the relations between nuclear and funda- 
mental parameters and the three GUT relations for the fundamental parameters in 
our BBN simulation code. Then we run the code for a set of parameters A In a and 
obtain the full nonlinear parameter dependence of abundances on variations of fun- 
damental parameters. Note that this method is impractical to investigating the full 
parameter space: fundamental parameters span a six-dimensional and nuclear param- 
eters a 13-dimensional parameter space which cannot be analyzed numerically due to 

^Updating to WMAP5 values does not lead to any significant change, there the 2a range is 
5.9 < IQiOr; < 6.5. 
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Figure 9.1: Variation of primordial abundances with a in three GUT scenarios. Red 
Hnes show the first 7 = scenario; green hnes the second 7 = 1 and blue lines the third 
7 = 1.5 scenario. Highlighted regions give the observational la limits (as explained 
in Sec. 16.8.31 no observational limits can be given for ^He). Error bars indicate the 
standard BBN abundances with theoretical la error [Serpico0"4| , for three different 
values of ry about the WMAP central value, as indicated on the upper horizontal axis. 
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time reasons. It is only practicable if the dimensionality of the parameter space is 
reduced by applying unification relations. 

For most nuclear parameters Xi the dependence on fundamental parameters Gk 
is only known to linear order or the nonlinear dependence involves a high uncertainty. 
However, we tested that for the unified models considered here, the fractional vari- 
ations in the nuclear parameters Xi remain small, well below 0.1. A linear approx- 
imation for the relation between nuclear and fundamental parameters is therefore 
appropriate. The main nonlinear effects enter at the level of nuclear reactions, in 
equations where nuclear parameters as binding energies enter in a known power law 
dependence. 

The nuclear parameters affecting most the large variation in ''Li abundance are 
mainly the deuterium and ''Be binding energies, with the ^He and ^He binding energies 
playing a smaller role. A decrease of Bd causes BBN to happen later, which means 
that the nucleon density is lower and reaction rates smaller. The abundances of 
A > A elements are rate- limited and thus decrease with decreasing Bd . This accounts 
for about two-thirds of the change in Yju. In addition, the cross-section of the 
^He(a,7)''Be reaction depends strongly on the Q-value, hence on the ''Be binding 
energy. Both these effects are computationally under control, therefore we believe 
that the specific nonlinear dependence in the scenarios we consider is well estimated 
within our code. We show in Fig. 19.21 the primordial abundances including nonlinear 
effects, i.e. without using a linear approximation for the relation between Ya and Xi. 
For our three GUT models we find a slightly different behavior of the ''Li abundance, 
which now has an approximately power-law dependence on variation of a. It is only 
slightly more difficult to bring the present observational abundances into agreement 
with standard BBN and the WMAP determination of rj ; still, if wc allow a variation 
of 0.00045 < A In a < 0.0005 in the 7 = 1.5 model, the predicted abundances are all 
very close to the la allowed regions. 
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Figure 9.2: Variation of primordial abundances with a in three GUT scenarios in- 
chiding nonhnear effects. Labels as in Fig. 19. f I 
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Chapter 10 

Experimental tests of 
variations 

In recent years possible variations in the constants of nature have been tested in vari- 
ous ways. Whilst direct laboratory measurements do not point towards any variation, 
some astrophysical tests yield slight variations. Still, those non-zero results neither 
constitute compelling evidence for variations since the applied fitting and analysis 
methods are still under debate, nor do the different non-zero results seem to be con- 
sistent with each othei0. Hence the question if there are and have been variations 
remains open from the experimental side. 

Here we review and discuss the observational data that we will consider in our 
effort to obtain a unified picture of time variation of couplings. We summarize the 
results that are most relevant for our analysis in Tab. 110. f I 

10.1 BBN 

The earliest processes for which Standard Model physics can be tested is BBN (z ~ 
10^°). The influence of varying constants on BBN has been laid out in chapters [7] to 
[S] In an extension of our previous treatment, we include in Sec. 110.2.11 the possible 
effect of varying constants at last scattering (formation of the CMB) on the input 
parameter rj of our BBN procedure. 

The uncertainty in the r] determination, rj — (6.20 ± 0.16) x 10^^" (WMAP5 plus 
BAO and SN, [WMAPS] ) yields a further correlated error for the abundances, which 
can be treated using the method of [Fiorentini98| . For any given set of fundamental 
variations we can define 

^ ^{Y, - Yr'nmAY, - Yf^), (10.1) 
with the inverse weight matrix 

(10.2) 



^For a thorough comparison of measurements of variations, a more profound theoretical back- 
ground is needed, as we will lay out in the rest of this thesis. 



2,r) 



bs. 
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where 

% - ^^^^ dluTj dlnrj [v ' ■ ^ ^ ^ 

We take, as in Tab. Ol 

aiii(D/H,yp,^Li/H) 



i91nr/ 



= (-1.6, 0.04, 2.1). (10.4) 



The l(2)cr error contour is given by x^/^ ^ 1(4) where v is the number of degrees of 
freedom. As the final abundances depend on variations of all fundamental constants, 
we have to evaluate the variations allowed by BBN for every model separately. 

In the light of complex astrophysics which may affect the extraction of the pri- 
mordial ^Li fraction, we also consider bounding the variations using deuterium and 
■^He alone (see also Sec. 19. ip . This yields a value consistent with zero for variations 
at BBN, since these abundances are consistent with standard BBN. 



10.2 CMB 

A further far-reaching test of varying constants is the cosmic microwave background 
(CMB) (z ^ lO'^). In principle, a and Gn are bounded by CMB observations. A 
variation of a affects the formation of the CMB through Thomson scattering and 
the recombination history. However, bounds on variations from CMB are typically 
very weak as there are significant degeneracies with other cosmological parameters 
[Martin sOSl IRocha03| . see also the discussion in Sec. 110. 2TT] Current bounds on a are 
[Rocha03] 

0.95 < ^2M^ < 1.02 {2a). (10.5) 

The CMB anisotropics may also be used to constrain the variation of Newton's con- 
stant Gn ■ The resulting bound depends on the form of the variation of Gn from the 
time of CMB to now. Using a step function one finds |Chan07| IZahn02| 

0.95 < < 1.05 (2cr), (10.6) 

where the instantaneous change in Gn may happen at any time between now and 
CMB decoupling. Using instead a linear function of the scale factor a, the bound is 

G 

0.89 < < 1.13 (2cr). (10.7) 

Note that here, as in most studies of time-dependent Gn, units are implicitly defined 
such that the elementary particle masses (and thus the mass of gravitating bodies, 
if gravitational self-energy is neglected) are constant. The relevant bound on dimen- 
sionless parameters concerns G-^^m^ = (TOAr/Mp)^(87r)~^. 
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10.2.1 Effect of "varying constants" at CMB and rj 

In our study on BBN we used the WMAP determination of the baryon number density 
parameter rj = ns/n-y directly to reduce by one the number of unknown parameters. 
However, we should also consider the effect of possible variations of fundamental 
parameters at the epoch of CMB decoupling, as CMB measurements are used to 
derive the value of ij. Hence, it seems important to study how possible variations 
at the epoch of CMB decoupling affect the determination of rj and hence also the 
outcomes of BBN simulations. 

Fundamental parameters affecting the CMB are the proton and electron masses, 
the gravitational constant and the fine structure constant, as well as the mass of any 
dark matter particle present. In Planck units, these reduce to the particle masses 
and a. The relevant cosmological parameters are the amplitude, spectral index (and 
possible running, etc.) of primordial perturbations; the baryon, dark matter, dark 
energy (cosmological constant, etc.) and curvature densities normalized to the critical 
density; the Hubble constant; and the reionization optical depth. Of these, the baryon 
density flbh^ will vary linearly with the proton mass in Planck units, for a fixed 
baryon-to-photon ratio ?/. Conversely, given a measurement of f^h/i^, the correct 
value of 77 varies inversely with the proton mass. The conversion factor between fl^h^ 
and ryio = 10^°?? is then (see Eq. ^J^) 

273.9imp^/G^)\o{mp^/G^)-^ ~ 273.9(1 - A ln(mw/Mp)|CMB) , (10.8) 

where we approximate the proton and neutron masses by their average mjy. 

If, therefore, we allow the proton mass (or the gravitational constant, in QCD 
units) to vary arbitrarily at the CMB epoch, t] is undetermined by WMAP and we 
must consider it as an extra free parameter or try to impose independent cosmological 
bounds. However, we impose that the size of variations away from the present value 
of rUp/Mp is a monotonically decreasing function of time: thus Aln(mAr/Afp)|CMB ^ 
Aln(mAr/Mp)|BBisi. Hence we would have a self-consistent treatment of this pa- 
rameter if the secondary discrepancies in primordial abundances due to an incor- 
rectly estimated 77 were smaller than the primary effect of varying mjsf/Mp at BBN. 
The maximum effect due to rescaling of rj would occur when Aln(mAr/Mp)|cMB — 
Aln(TOAr/Mp)|BBNi adding the nuclcon mass induced variation in 77, 

A In 77 = Aln(r77jv/Mp)|BBN • (10.9) 

In our treatment of BBN in the next chapter, we will consider this effect by studying 
two limiting cases. First, when A(77ijv/Mp)|c]viB ^ A(?77Ar/A/p)|BBN; then our pre- 
vious results hold. In the second case, with A(r7ijv/Mp)|cMB — A(r7ijv/Afp)|BBNi the 
value of 77 may be significantly rescaled. 

10.3 Quasar absorption spectra 

The observation of absorption spectra of distant interstellar clouds allows to probe 
atomic physics over large time scales. Comparing observed spectra with the spectra 
observed in the laboratory, together with the in general well known dependence of 
the spectra on fundamental constants, gives bounds on the possible variation of cou- 
plings. Different kinds of spectra (atomic, molecular, . . . ) are sensitive to different 
parameters, which we will list in the following paragraphs. 
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Tests of a 

Atomic spectra are primarily sensitive to a. Several groups using various methods of 
modeling and numerical analysis have published results; we quote here only the latest 
bounds. Murphy and collaborators |Murphy03.2| studied the spectra of 143 quasar 
absorption systems over the redshift range 0.2 < Zabs < 4.2. Their most robust 
estimate is a weighted mean 

— = (-0.57±0.11) X 10"^ (10.10) 

a 

Dividing the data into low {z < 1.8) and high (z > 1.8) redshift subsamples, they 
obtain 



z < 1.8, N,y, = 77, (zabs) = 1.07, = (-0.54 ± 0.12) x lO"'^ 

a 

z>l.8, Nsy, = 66, {zabs) = 2.55, ^ - (-0.74 ± 0.17) x 10"^, (10.11) 

where N^ys is the number of absorption systems in the sample and (zabs) is the 
averaged sample redshift. 

In discussing unified models in Sec. 111.31 we will define various "epochs" for the 
purpose of collating data and comparing them with models over certain ranges of 
redshift. The 143 data points are then assigned to different epochs: we choose to put 
boundaries at z = 0.81 and z = 2.4, thus we obtain three sub-samples 

z<0.81, iV,j„ = 18, (z)=0.65, ^ = (-0.29 ±0.31) X lO"'^ 

Aa 

0.8K z < 2.4 Nsys=S5, (z) = 1.47, = (-0.58 ± 0.13) x 10"^ 

a 

Aa 

z > 2.4, N^ys = 40, (z) = 2.84, — = (-0.87 ± 0.37) x 10"^ (10.12) 

Here we have used the "fiducial sample" of Murphy03.1| , the weighted average has 
been taken, and we have included |Murphyprivate| the 15 additional samples used in 
|Murphy03.2] . For convenience we will refer to these results as "Ma" . 

Further results have been obtained by Levshakov et al. [Levshakov07. 1| . and re- 
ported in [Fujii07| : 

— = (-0.01 ± 0.18) X 10"^ Zabs = 1.15 

a 

Aa 

= (0.57 ± 0.27) X 10"^ Zabs = 1.84. (10.13) 

a 

We note that the value for z = 1.84 has an opposite sign of variation to the Ma result, 
though the variation does not have high statistical significance. The observational 
situation is clearly unsatisfactory. 



Tests of 

Vibro-rotational transitions of molecular hydrogen H2 are sensitive to /i = mp/me. 
From H2 lines of two quasar absorption systems (at z = 2.59 and z = 3.02) a variation 
is found [Reinhold06] of 

— = (2.4 ± 0.6) X 10~^ (10.14) 
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taking a weighted average. We will refer to this result as "R/i" after Reinhold et al. 
The individual systems yield |Reinhold06| 

— = (2.78 ± 0.88) X 10"^ zabs = 2.59 

— = (2.06 ± 0.79) X 10"^ Zahs = 3.02. (10.15) 

Recently the z = 3.02 system has been reanalyzed ' WendtOS] . with the result that the 
claimed significance of Eq. (|10.15p was not reproduced, and the absolute magnitude 
of the variation is bounded by |A/i//i| < 4.9 x 10^^ at 2ct, or 

|A/i//i| < 2.5 X 10^^ z„fc, = 3.02 {la). (10.16) 

Very recently, a new determination of the variation of /i appeared |King08| report - 
ing a reanalysis of spectra from the same two H2 absorption systems as jReinhold06j , 
and adding one additional system at z ~ 2.8. The results of the new analysis are not 
consistent with the previous claim indicating a nonzero variation, either considering 
all three systems or the two previously considered. Instead, |King08| obtain a null 
bound, A/i//i = (2.6 ± 3.0) x 10~^. As these results were published after this study 
has been finished, they are not considered any further. 

The inversion spectrum of ammonia has been used to bound precisely at lower 
redshift [FlambaumNH3 . Recently the single known NH3 absorber system at cosmo- 
logical redshift has been analyzed [MurphyNH3| , yielding 

— = (0.74 ± 0.89) X 10"^ z = 0.68. (10.17) 

Tests of y 

The 21cm HI line and molecular rotation spectra are sensitive to y = a^gp, where gp 
is the proton g-factor. Bounds on this quantity from [MurphyOT] are 

— = (-0.20 ±0.44) X 10"^ z = 0.247 

y 

— = (-0.16 ±0.54) X 10"^ z = 0.685. (10.18) 

y 

Tests of X 

Further, the comparison of UV heavy element transitions with HI line probes for vari- 
ations of a; = a^gp^~^ |Tzanavaris06j : the weighted mean of nine analyzed systems 
yields 

— = (0.63 ± 0.99) X 10"^ 0.23 < Zabs < 2.35. (10.19) 

X 

However, we note that i) the systems lie in two widely-separated low-redshift (0.23 < 
z < 0.53) and high-redshift (1.7 < z < 2.35) ranges; and ii) these two sub-samples 
have completely different scatter, x^/i^ about the mean for the low- and high-redshift 
systems being 0.33, and 2.1, respectively. Hence we consider two samples, with average 
redshift z = 0.40 (5 systems) and z — 2.03 (4 systems). With expanded error bars in 
the high-redshift sample (after "method 3" of [Tzanavaris06| ) we find 

— = (1.02 ± 1.68) X 10^^ (z)=0.40 

X 

^ = (0.58 ± 1.94) X 10"^ (z) = 2.03. (10.20) 
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Tests of F 

The comparison of HI and OH lines is sensitive to changes in F = gp [a^/i] 
[KanekarOS] and yields 

AF 

— = (0.44 ±0.36'*°' ±1.0''^^) X 10"^ z = 0.765. (10.21) 
F 

Tests of F' 

A similar method comparing CH and CO lines has very recently been proposed at 
high redshift [Levshakov07.2| yielding the best bound at redshifts > 4.5. The following 
bounds on F' = a^//i are obtained for two systems: 

(0.1± 1.0) X IQ-^, z = 6.42 

(1.4 ± 1.5) X IQ-^, z = 4.69. (10.22) 



AF' 
AF' 



10.4 The Oklo natural reactor 

In Oklo/Gabon, a natural fission reactor formed by naturally enriched uranium in 
a rock formation with a water moderator was operating about 2 billion years ago 
{At ~ 1.8 X 10^ y, z ^ 0.14 with WMAP5 best fit cosmology). The resulting isotopic 
ratios in this rock nowadays differ radically from any other terrestrial material. By 
modeling the nuclear fission process, one can in principle bound the variation of a 
over this period. The determination of A In a at the time of the reactions results 
from considering the possible shift, due to variation of electromagnetic self-energy, in 
the position of a very low-lying neutron capture resonance of ^^^Sm. The analysis of 
|Petrov05] gives the bound (taken as la) 

- 5.6 X 10"® < Aa/a < 6.6 X 10"^ (10.23) 

For a linear time dependence this results in the bound 

\a/a\ < 3 X 10""y-^ (10.24) 

Note that these results concern varying a only. If other parameters affecting nu- 
clear forces, in particular light quark masses, are allowed to vary, the interpretation 
of this bound becomes unclear [Olive021 IFlambaum02] since it depends on a nuclear 
resonance of ^^°Sm whose properties are very difficult to investigate from first princi- 
ples. In the absence of a resolution to this problem we consider Oklo as applying only 
to the a variation in each model. In scenarios where several couplings vary simulta- 
neously we do not consider strong cancellations. Nevertheless, we allow for a certain 
degree of accidental cancellation and therefore multiply the error on the bound Eq. 
I|10.23p by a factor three. 
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10.5 Meteorite dating 

Meteorites which have formed at about the same time as the solar system, ^Met — 
4.6 X 10^ y ago (z ~ 0.44) contain long-hved a- or /3-decay isotopes. The decay rates 
of those isotopes may be sensitive probes of cosmological variation [Olive02|, lOliveOS) 
ISisternaQO] . Their (generaUy) smaU Q- values result from accidental cancellations 
between different contributions to nuclear binding energy, depending on fundamental 
couplings in different ways, thus the sensitivity of the decay rate may be enhanced by 
orders of magnitude. 

The best bound concerns the ^^^Re /3-decay to osmium with = 2.66 keV. The 
decay rate Aist is measured at present in the laboratory, and also deduced by isotopic 
analysis of meteorites formed about the same time as the solar system, 4.6 x 10^ years 
ago. More precisely, the ratio Aigy/Au, averaged over the time between formation and 
the present, is measurable [Olive03[ |Fujii03| , where Au is the rate of some other decay 
(for example uranium) used to calibrate meteorite ages. 

The experimental values of Aigy imply (setting Xu to a constant value) 

*Mct /° ^^i^d^ = 0.016 ±0.016. (10.25) 

Since the redshift back to ^Mct is relatively small, we obtain bounds on recent time 
variation by assuming a linear evolution up to the present, for which the left hand 
side is — (tMct/2)Ai87/Ai87 and the fractional rate of change is bounded by 

^ ~ (-7.2 ±6.9) X 10-12 y-i. (10.26) 

Al87 

Projected back to ^Mot this gives the bound 

A In Ai87 =i 0.033 ±0.032 (z ~ 0.44). (10.27) 

This is a conservative bound unless the time variation has recently accelerated, or 
there are significant oscillatory variations over time. 

Since the possible dependence of "control" decay rates Xtj/itin on nuclear or fun- 
damental parameters is much weaker than that of Xisr/rri]^, we use this result for the 
variation of Ai87 in units where Au is constant, i.e. Aln(Ai87/Au) — Aln(Ai87/?7iAr). 
We find the decay rate dependence to be [DSW08.I] 

Aln — ~ -2.2xl04Alna-1.9xl0''Aln-^+2300Aln— ^2 530 A In 

'7^^f ^QCD J^QCD ^QCD 

(10.28) 

10.6 Bounds on the variation of 

Variations of Newton's constant have been studied in the solar system and in astro- 
physical effects. Whilst all references give bounds exclusively on a potential variation 
of Gn, one should note that besides Gn also nuclear parameters (neutron / proton 
masses and parameters of nuclear forces) can vary, which would in general add de- 
generacies and make the results less stringent. It has generally been assumed that 
particle masses are constant, thus the resulting bounds actually constrain variation 
of G^m\i oc {rriN /M-pY . 
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In the solar system, changes of Gn induce changes in the orbits of planets. Range 
measurements to Mars from 1976 to 1982 can be used to obtain Hellings83] 



Gn/Gn = 2±4 X lO^^V^- (10.29) 

Lunar laser ranging from 1970 to 2004 yields [Williams04j 

Gn/Gn = (4 ± 9) X 10-"y-\ (10.30) 

The stability of the orbital period of the binary pulsar PSR 1913+16 [DamourSS] may 
be used to deduce 

Gn/Gn (1.0 ± 2.3) X 10~"y-^ (10.31) 

All these results apply at the present epoch z = 0. 

A bound on the behavior of Gn over the lifetime of the Sun (approximately 4.5 x 
lO^y, z = 0.43) was found by Guenther et al. |Guenther98| by considering the effect 
of the resulting discrepancy in the helium/hydrogen fraction on p-mode oscillation 
spectra. The claimed constraint is 

|Gn/Gn| < 1.6 X 10-12 y"' 

lAlnGNl < 7.2 X 10-3 z = 0.43, (10.32) 

where the assumed form of variation is a power law in time since the Big Bang, which 
may be approximated over the last few billion years as a linear dependence. For 
models with significantly nonlinear time dependence the bound may be reevaluated: 
since the bound arises from the accumulated effect of hydrogen burning since the 
birth of the Sun, it may be expressed as an integral of the variation over the Sun's 
lifetime analogous to Eq. (|10.25p . 

The mass of neutron stars is determined by the Chandrasekhar mass 

where m„ is the neutron mass. This may be reexpressed in terms of the baryon 
number of the star ub oc Mc\i/mn oc (Gn™^)"'^^^, which is expected to be constant 
up to small corrections from matter accreting onto it. Thus the relative masses of 
neutron stars measured at the same epoch probes the fractional variation of G^ra^ 
between their epochs of formation. From the comparison of masses of young and old 
neutron stars in binary systems (where the oldest neutron stars are up to 12 Gy old, 
z ^ 3.3), it is found |Thorsett96| that the variation of the average neutron star mass 
is /i„ = -1.2 ±4.0(8.5) X lO-i2M0y-i at 60% (95%) confidence level. In units 
where particle masses are constant, we have 

Gn/Gn = -0.6 ± 2.0 (4.2) x lO^^^y'S (10.34) 

where the averaging is performed over the last 12 x lO^y, and the bound should be 
reinterpreted for variations which are not linear in time. The absolute variation over 
this period is then bounded at la as 



AIuGn = (-0.7±2.4) X 10-^ z = 3.3. 



(10.35) 
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10.7 Atomic clocks 

Atomic transitions can be measured in the laboratory to very high precision over peri- 
ods of years. As each atomic transition depends differently on fundamental constants 
{e.g. a, fi), comparisons of different atomic transitions over long periods of time give 
very sensitive results. 

Recently, stringent bounds on the present time variation of the fine structure 
constant and the electron- proton mass ratio have been obtained by [BlattOS] , 



Fortier et al. |Fortier07| obtain stronger bounds, \a/a\ < 1.3 x 10^^^ y^^, if other 
relevant parameters are assumed not to vary. If other atomic physics parameters are 
allowed to vary, this bound becomes considerably weaker, depending on a possible 
relative variation of the Cs magnetic moment and the Bohr magneton. Direct com- 
parison of optical frequencies may yield bounds at the level of 10""'^^ per year; limits 
on variation of a from this method are reported with uncertainty 2.3 x 10~^^y~^ 
|Rosenband08| but designated as preliminary. If these bounds are used then our 
limits from atomic clocks via a variation should be tightened by about an order of 
magnitude. 

Extrapolating the results of [BlattOSj to the time of Oklo {z = 0.14, t = l.Sx lO^y) 



d\na/dt= (-0.31 ±0.3) x 10"^^ y~^ 
d\n^i/dt = (1.5 ± 1.7) X 10-^^ y-^. 



(10.36) 



gives 



Alna = (-0.56 ±0.54) x 10 
Aln^ = (-0.27 ±0.31) x 10 



,-6 



,-5 



(10.37) 



Method 


redshift 


A In a 

[10-6] 


Aln^ 

[10-5] 


A In Gt^m'j^ 
[10-2] 


A Inx 

[10-5] 


A\ny 

[10-5] 


Aln^^ 

[10-5] 


Ahii^' 

[10-4] 


AlnAi87 

[10-2] 


Oklo a 'Petrov05| 


0.14 


0.00 ±0.06 
















21cm 1 Murphy 01 
Sun (Guenther98] 
Heavy /HI. low-z |Tzanavaris06j 
Meteorite [OhveOSj 
Ma epoch 2 'Murphy03.2" 
Ammonia [FlambaumNHS, 
21cm (Murphy 01] 

ni / v^n ivdjiicAdji uoi 


0.247 
0.43 
0.40 
0.44 
0.65 
0.68 
0.685 
n 7fi^ 

U. ( UJ 


-2.9 ±3.1 


0.06 ±0.19 


0±0.72 


1.0± 1.7 


-0.20 ± 0.44 
-0.16 ± 0.54 


zL -1- 1 1 




3.3 ±3.2 


Absorption ,Fujii07, 
Ma epoch 3 Murphy03.2 
Absorption ILevshakovOT.lJ 
Heavy/Hi, high-z TzanavarisOS 


1.15 
1.47 
1.84 
2.03 


-0.1 ± 1.8 
-5.8 ± 1.3 
5.7±2.7 






0.6 ± 1.9 










H2 [ReinholdOe 
Ma epoch 4 |Murphy03.2| 

H2 ;Reinhold06j 
Neutron stars (Thorsett96j 
CII / CO 'Levshakov07.2; 
CII / CO J.cvshakov07.2 


2.59 
2.84 
3.02 
3.3 
4.69 
6.42 


-8.7±3.7 


2.78 ±0.88 
2.06 ±0.79 


-0.7 ±2.4 








1.4± 1.5 
0.1 ± 1.0 




CMB |Martins03|. |Chan07| 


10^ 


'-'-3xl0*> 




0+^ 













Table 10.1: Observational la bounds on variations. Observables are defined a,s fi = nip/me, x = a^gpfi^^, y = a^gp, F = gp[a'^fi]^'^'^, 
F' = o? I [i. The given redshift may denote a single measurement, or an averaged value over a certain range: see main text. The two 
CMB bounds are independent of each other. Our BBN bounds cannot be displayed in this form. 



Chapter 11 



Variations from BBN to 
today in unified scenarios 

As has been laid out in Sec. 14.51 we will use the concept of grand unification to 
reduce the number of potentially varying parameters. The main idea is that the GUT 
relations interrelate variations of the "fundamental" parameters Gk which we defined 
in Sec. O 

In this chapter, we consider the hypothesis that, for all redshifts, all fractional 
variations in the "fundamental" parameters Gk are proportional to one nontrivial 
variation with fixed constants of proportionality. If the variation of the unified gauge 
coupling A In ax is nonvanishing, we may write 

AlnGfc = 4Alnax (H-l) 

for some constants d^, assuming small variations. Different unification scenarios cor- 
respond to different sets of values for the "unification coefficients" dk- Considering 
the values of AlnG^ as coordinates in an A^fc-dimensional space, this assumption re- 
stricts variations to a single line passing through zero. The variation then constitutes 
exactly one degree of freedom. We will go beyond this hypothesis in the next chapter 
where we also consider models with growing neutrinos and oscillating variations (see 
Sec. I5.4p for which a fixed linear relation (jll.ip is not realized for all z. 

11.1 GUT relations 

GUT relations have the property that variations of the Standard Model gauge cou- 
plings and mass ratios can be determined in terms of a smaller set of parameters 
describing the unified theory and its symmetry breaking. Hence, if nonzero varia- 
tions in different observables are measured at similar redshifts, models of unification 
may be tested without referring to any specific hypothesis for the overall cosmological 
history of the variation. We need only assume that for a given range of z the time 
variation is slow and approximately homogeneous in space, hence A In ax depends 
only on redshift z to a good approximation. The relevant unified parameters are the 
unification mass Mgut (relative to the Planck mass), the GUT coupling ax defined 
at the scale MquT: the Higgs v.e.v. (0) and, for supersymmetric theories, the soft 
supersymmetry breaking masses to, which enter in the renormalization group (RG) 
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equations for the running couplings. Then, for the variations at any given z we can 
write 

A In "^'^^^ =dMl, Ahiax^dxl, A In Jf^ = dnl, A\n—^ — = dsl, 

Mp MguT i^'J^GUT 

(11;2) 

where l{z) is the "evolution factor" introduced for later convenience. If ax varies 
nontrivially we may normalize I via dx = 1. In supersymmetric theories we set 
ax = 1/24, in nonsupersymmetric theories we set ds = Q and ax = 1/40 (see Sec. 



We make the simplifying assumption that the masses of Standard Model fermions 
all vary as the Higgs v.e.v., i.e. Yukawa couplings are constant at the unification scale: 

Aln-^ = Aln-^=Aln-^ = Aln-^ = Aln-^. (11.3) 

Mgut J^Igut Mgut Mgut AIgut 

Like the gauge couplings, also the fermion masses vary under variations of the unified 
coupling ax due to the renormalization group running of fermion masses. However, 
we have explicitly calculated the effect of varying couplings and found that it is at the 
order of a 1% correctioiQ, which is already smaller than our uncertainties in hadronic 
and nuclear physicsH. Hence we can apply the assumption pi.3p . Using the relations 
([4231) and ([42^1) . one finds for the QCD scale 

AlnjAgc d/Mgut) 2n 2 4 

1 ^ g^'^^ + 9^^ + 9^^ ^^'-^^ 

and for the fine structure constant, 

A In a 80a , 43 a , 257 a , 

= dx -\ dn -\ ds. (11-5) 

I 27ax 27 27r 27 27r ^ ^ ^ 

For the nucleon mass we include possible strange quark contributions. In our 
treatment of BBN, we neglected strange quark contributions, as the final dependence 
on rris was much below the model uncertainty. Here we include the roughly known 
strange contribution to the nucleon mass. The uncertainty in the strangeness content 
is an indicator of the overall uncertainty that may arise due to rris variation. We 
found (Eqs. (jOHl . KT^i and jHIl) 

Aln-!!^ = 0.048Aln-^ + (0.12 ± 0.12) A In (11.6) 

^QCD AqCD AqGD 

Aln-^^ = -0.59Alna + 1.59Aln— (11.7) 

^QGD AqGD 



and thus 



^^^^ - (0.58 T 0.08)— + (0.37 T 0.05)^5+ (-0.65 ±0.09)dij, (11. 



I ax 
Aln(GNm^) , n .a^n^.^d 



2dM + (1.16 T 0.17)— + (0.74 t 0.11)^5 + (0.71 ± 0.19)dH, 
I ax 

(11.9) 



^For low-energy observablos such as rnq{Q^) / Kqcjj we consider an RG scale that is fixed 
relative to Aqqd, ~ ^qc d ■ Thus the variation of mq(Q^) /mq(AlQjjrp) is entirely due to 

the dependence on asiMQux), which is suppressed by a loop factor ax /ti" compared to the nonper- 
turbative dependence of f^QCD /^GUT on ax- We find A ln(rnq(Q^)/mq(Mg,j^j,)) = 2/7Alnax — 
(9ax/77r) A ln(AQc7£)/MG[/T) under variation of axi where rhq is the running quark mass. 

^Langacker et al. |Langacker01| arrive at the same conclusion. 
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where the upper or lower signs correspond to the positive or negative signs in Eq. 
(|11.6p respectively. 

The largest contribution to variations of the proton g-factor gp has been argued 
to arise from the pion loop |Murphy03.2| , yielding at first order a dependence on the 
light quark mass of 

Alngp ~ -0.087Alnm/AQCD, 

^}l}^^Qm—~0.Q7dH + QMds. (11.10) 
I ax 

Hence the variations of observables including gp are 

= (-0.48 ± 0.08)— + (0.59 T 0.09)^^ + (-0.31 ± 0.05)ds 

I ax 

= 0.10— - G.Oerfff + G.06rfs 

I ax 

= (1.04 =F 0.13)— + (-1.08 ± O.U)dH + (0.65 T 0.08)ds 
I ax 

; = (-0.54 ± 0.08)— + (0.65 T 0.09)^^ + (-0.35 ± 0.05)ds. (H-U) 

I ax 

We have now expressed the variations accessible to observation in terms of three 
(four) variables: I, dx, dn (and ds), where one parameter may be eliminated by 
normalization. Different unified scenarios will be characterized by different relations 
among these parameters. 



11.2 Variations in six different unified scenarios 

We will now investigate six different scenarios for the variation of the grand unified 
parameters ax, Mgjjt/Mp, [cf)) / Mqjjt and rh/MQUT- These will fix the unification 
coefficients dk ■ For each unified scenario we display the z-dependence of the fractional 
variation. Each figure shows the available information from observations of different 
couplings, interpreted as constraints on the variation of a single parameter. Since we 
have only one free variable we can plot all observations simultaneously as a function 
of redshift. Inspection "by eye" permits to judge if a smooth and monotonic evolution 
of the varying parameter is consistent or not. Most data points are upper bounds 
on a possible variation, and for the non-zero variations we can study two immediate 
questions in each scenario: 

First, whether claimed nonzero variations of a |Murphy03.2| and fjL |Reinhold06| 
at redshift 2-3 are compatible with one another, since the ratio of their fractional 
variations is predicted in each scenario. 

Second, we consider whether there is an indication of nonzero variation at BBN. 
For no variation at BBN we obtain = 17.9 for 3 measured abundances ("^He, 
D, ^Li). This discrepancy between theory and observation is exclusively due to ''Li. 
(Considering only "^He and D, the value of is 0.24.) If we wish to solve or ameliorate 
the "lithium problem" by a nonzero variation, we will require /v to be not much 
larger than unity, taking v = 2 as appropriate for one adjustable parameter. If there 
is no significant range where the three abundances have a 2a fit [x^/v < 4) then 
we give up the hypothesis that the ^Li problem is solved by coupling variations and 
instead assume that the observed depletion is due to some astrophysical effect. In 
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this case we consider only D and ^He abundances as observational bounds on the size 
of variations at BBN. 



11.2.1 Varying a alone 

Before describing the six different grand unified scenarios, we consider a variation of 
the fine structure constant a alone. Clearly here we are unable to account for any 
nonzero variation in /i or other quantities independent of a. The cosmological history 
is dominated by the nonzero variation of the Ma values at redshifts z ~ 1 to 4. We 
find that there is almost no 2a match of the BBN values jv > 3.9): the 2-sigma 
range is 

3.25% > AlnasBW > 4.06%. (11-12) 

Hence it seems unlikely that the "lithium problem" can be solved by a variation of a 
alone. If we regard the ^Li discrepancy as due to systematic or astrophysical effects 
we can set a conservative bound on a variation from ^He and D abundances 

-3.6%> AlnaBSJV > 1-9%, (11.13) 

where we imposed that neither the D nor *He abundance should deviate by more than 
2(7 from observational values. See Fig. Ill.ll for a summary of the bounds in this case. 



11.2.2 Scenario 1: Varying gravitational coupling 

In this scenario we have only d]\{ nonvanishing, 

dn = ds = dx = 0, (11.14) 

therefore 

Aln - 

Mp 2 



^GUT 1 , , ^ ,2 

We find that there is no value of AlnGNAnrri for which BBN is consistent with 



^ln^^^ = -AlnGNA^CD- (11-15) 



the three observed abundances within 2(7. The best fit values are x^/'^ ^ 7.7 for no 
variation of ttin/Mp at CMB and x^/i^ > 5.9 if the variation of rriN/Mp has the 
same size at BBN and CMB. Assuming that the discrepancy in the ^Li abundance is 
due to some other effect, we find the allowed region of variation of Gn at BBN under 
which primordial D and ^Hc abundance lie within the observed range at la {2a), 

- 5% (-13%) < AIuGnAqcd < 12% (22%) (11.16) 

If the variation of mjv/Afp has the same size at BBN and CMB one finds 

-4% (-11%) < AIuGnA^c^ < 10% (16%). (11.17) 

The bounds on time variation of GnAq^^^ are much weaker than for many other 
varying couplings. This scenario also predicts a vanishing value of rj in Eotvos exper- 
iments (see Sec. 112.31 for details). Thus, to any one of the following scenarios we may 
add an additional nonzero dM of similar size to dx, dn or ds without changing the 
results significantly. 
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Figure 11.1: Variations for varying a alone. Only observations constraining a vari- 
ation are shown; the BBN fit including ^Li is poor jv > 7.8/2) hence we also 
display a conservative bound from ''He and D abundances neglecting ^Li. 
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Figure 11.2: Variations for scenario 2; BBN bounds are 2a bounds. 
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11.2.3 Scenario 2: Varying unified coupling 

In the first GUT scenario without SUSY we consider the case when only dx is non- 
vanishing, 

dn ^ds = 0, ax = 1/40. (11.18) 

Within a supersymmetric theory the same relations will apply except that ax = 1/24 
and the variations of observables are scaled by a factor 24/40 relative to A In ax : we 
designate this as Scenario 2S. In both cases we find here 



A In a 



27. (11.19) 



It is then highly unlikely for the nonzero Ma result for variation of a to coexist 
with the determination of fi at redshift around 3 [Reinhold06| . even if the latter is 
interpreted as an upper bound on the absolute size of variation WendtOS] . 

For the BBN fit, we find without SUSY (excluding modifications of the baryon 
fraction rj due to varying toa?) no range of values fitting at la level jv > 2.3). At 
2(7 the abundances, including ^Li, become consistent for the range 

- 5.7 X 10""* < Alnax < -1.7 X 10~^ (2cr). (11.20) 

If one includes a variation of mx at the time of CMB with the same magnitude as at 
BBN the result remains unchanged {x^/j^ > 2.45), with the same 2a range. For this 
scenario we may consider a nonzero variation at BBN, but more recent probes must 
all be viewed as increasingly tight null bounds. 



11.2.4 Scenario 3: Varying Fermi scale 

In this scenario we consider the case when the variation arises solely from a change 
in the Higgs expectation value relative to the unified scale, thus only dn is nonzero: 

ds = dM = dx = 0, ax = 1/40. (11.21) 

This scenario implies 

^ = -325. (11.22) 

A In a ^ ' 

Whether we interpret the determination of fji |Reinhold06| as a detection or an upper 

bound, any variation in a at large redshift should be orders of magnitude smaller 

than current observational sensitivity. 

We find for BBN including ''Li {u = 2) no la range {y^ /v > 1.95) but 

6 X 10"^ < Aln(0)/AfGc/T < 22 X 10"^ {2a). (11.23) 

A variation of mx at the time of CMB with the same magnitude as at BBN does not 
change this result. 



11.2.5 Scenario 4: Varying Fermi scale and SUSY-breaking 
scale 

This scenario corresponds to scenario 3, but includes supersymmetry and assumes 
that the mass-generating mechanism for SM particles and their superpartners gives 
rise to the same variation: 

dM = dx=0, ds = dH, ax = 1/24. (11.24) 
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We find here 

^--21.5, (11.25) 
A In a 

such that again the claimed nonzero variations in a and fi cannot be compatible and 
the variation in a at redshift 3 must be below current sensitivities. We demonstrate 
this in Fig. 111.41 where we show for this scenario the bounds on the variable dnl = 
Ahi{{(p) /Mgut) that arise from various observations. 

We find for BBN including ^Li(j/ ^ 2) no Icr fit {x^/iy > 1.60), while at 2a 

1.25 X 10"^ < Aln((/))/AfG;7T < 5.4 X 10"^ {2a). (11.26) 

If one includes a variation oi ttin at the time of CMB with the same magnitude as at 
BBN the allowed range becomes slightly restricted jv > 1.72), 

1.20 X 10^2 < A\n{(j))/MGUT < 4.9 x 10"^ {2a). (11.27) 



11.2.6 Scenario 5: Varying unified coupling and Fermi scale 

In this scenario we study a combined variation of the unified coupling and the Higgs 
expectation value: 



M = ds = 0, dH = ldx, ax = 1/40. (11.28) 

A ln(0) /Mgut 
A In Aqcd /Mgut 



The parameter 7 can be related to the parameter 7 = in{<p)/MGUT vi'hich was 



/ 27r 2 , , , 



introduced in Sec. 

\y 9ax ' 9 

There we examined the cases 7 = (0, 1, 1.5) which correspond to 7 = (0, 36, 63). Here 
we find that the best BBN fit is reached for 7 « 50 with x^/'^ — 1-45. Note that we 
have the freedom to adjust 7 such that nonzero variations of a and /i at redshift ~ 3 
are consistent with each other. We have 

AluM ^ 23.2- 0.657 . . 

A\na 0.865 + 0.0027 ■ ^ ' ' 

We choose for illustration 7 = 42, for which 

Aln^ = -5.6Alna (11.31) 

and the 2a contour for BBN is 

7.5 X 10"'* < A In ax < 28 X 10"''. (11.32) 

For a variation of at the time of CMB with the same magnitude as at BBN the 
fit becomes worse {y^ jv > 1.68). However, a 2a fit to BBN is obtained over a wide 
range of < 7 < 26 (negative A In ax) and 40 < 7 < 00 (positive A In ax)- 

Assuming that the apparent ^Li mismatch at BBN is due to systematic astrophys- 
ical effects, we may bound ax with only D and ^He abundances. Here we find at 
Icr 

- 5.5 X 10~* < Alnax < 1-44 x 10"^ (11.33) 

In Fig. lll.5l we again plot simultaneously all observations for this scenario. This shows 
that the bound from BBN including ^Li is not consistent with the claimed nonzero 
variations of a and fi for a monotonic evolution over z. 



11.2. VARIATIONS IN SIX DIFFERENT UNIFIED SCENARIOS 



89 



100- 

80- 

60- 

40- 

20- 

0- 

-20- 

-40- 

-60- 

-80- 

-100- 
0,1 



• Met 

G 

X 

y 

F 
F' 

BBN 



CMB BBN 



Figure 11.3: Variations for scenario 3; BBN bounds are 2a. Note that due to the very 
large ratio Aln/z/Alna in this scenario, points indicating any nonzero variation of 
a fall well outside the range of the graph. 
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Figure 11.5: Variations for scenario 5, 7 = 42; BBN bounds are 2a 
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Figure 11.6: Variations for scenario 6, 7 ~ 70; BBN bounds are 2a 
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Figure 11.7: Variations for scenario 6, 7 = 25; BBN bounds are 2a 



11.2.7 Scenario 6: Varying unified coupling and Fermi scale 
with SUSY 

In this scenario we study a combined variation of the unified couphng and the Higgs 
v.e.v. including SUSY, where as in Scenario 4 we tie the variations of the superpartner 
masses and Fermi scale together: 

dM = 0, ds^dn^ldx, ax = 1/24. (11.34) 

Now the relation to 7 is modified as 



/ 27r 2 
7 - 



One may again adjust 7 to make nonzero variations in a and /i self-consistent, where 
now 

Aln^ 14-0.287 



A In a 0.52 + 0.0137 



(11.36) 



We find that a good fit to BBN is obtained over a large range of 7, ranging from 
7 = 100 to infinity with minimal x^/v = 1.45. This shows that the main effect in 
the SUSY model comes from the variation of the Higgs v.e.v. Including a variation 
of m^r at the time of CMB with the same magnitude as at BBN the fits gets worse 
(x^/iy > 1.8). A 2(7 fit can be obtained for < 7 < 28 (for negative Alnax at BBN) 
and for 58 < 7 < cx) (positive A In ax)- 

First, we study the case 7 = 70 for which 



A In /I = -3.9 A In a (7 70) 



(11.37) 
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and BBN is fit with a 2a range 

5.5 X 10"* < A In ax < 18 X 10"''. 
Neglecting ^Li, we obtain a la bound from BBN 

~ 3.5 X 10"* < Alnax < 9.3 x 10" 



(11.38) 



(11.39) 



Secondly, we study the case 7 = 25 where 



Aln^ = 8.3Alna 



(7 = 25), 



(11.40) 



and where the 2a contour for BBN is 



13 X 10"* < Alnax < -7 x 10 



,-4 



(11.41) 



In this second case the Murphy a measurement and BBN point into the same direc- 
tion. The difference between the two values of 7 can be seen from a comparison of 
Figs. nH and [III71 

11.3 Epochs and evolution factors 
11.3.1 Epochs 

As a next step, we group the information on experimental bounds on variations of 
couplings in the different unified scenarios (displayed in Figs. lll.l1 to lll.7]) into differ- 
ent cosmological epochs. This produces a first quantitative estimate of the possible 
time evolution for the various unified scenarios. The choice of epochs is somewhat 
arbitrary. Two epochs are singled out by events in early cosmology, namely the last 
scattering surface of CMB, and BBN. The very recent epoch comprises present day 
laboratory experiments and the Oklo natural reactor, for which a linear interpola- 
tion to the present rate of varying couplings seems reasonable. We further divide the 
observations at intermediate redshift into three epochs. 

• Epoch 1: Today until Oklo 

Contains Oklo and laboratory measurements. For the laboratory measurements, 
we extrapolate the rate of change of the couplings to finite changes at the redshift 
z = 0.14 {t = 1.8 x lO^y) of the Oklo event. 

• Epoch 2: 0.2 < z < 0.8 

Contains absorption spectra and isotopic abundance measurements in mete- 
orites. We chose a boundary z — 0.8 since the Murphy dataset [Murphy03.2| 
has relatively few systems around this redshift, making a natural division. 

• Epoch 3: 0.8 < z < 2.4 

Contains several absorption spectra measurements. The end of the Tzanavaris 
dataset |Tzanavaris06| sets the cut at z = 2.4. 

• Epoch 4: 2.4 < z < 10 

Contains absorption spectra measurements and bounds on Gn from neutron 
stars. 

• Epoch 5: CMB, z w 1100 



Epoch 6: BBN, z w 10 
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II. 3.2 Evolution factors 

We define "evolution factors" Z„ for epochs n = 1, . . . , 6 by 

A\nGk,n = dkln. (11.42) 

For each unification scenario we proceed to a quantitative estimate of shown in Tab. 

III. 11 The usefulness of considering the evolution factors In is that the unknown (and 
possibly not monotonic) behavior of the mechanism driving the coupling variations is 
rolled into a finite number of parameters. For a monotonic behavior they satisfy In < 
/p whenever z„ < Zp. The basic assumption remains the proportionality A In G'fe(z„) — 
dkl{zn) — dklm with constant unification coefficients dk independent of the epoch. 
The normalization of In is arbitrary, and we take for scenarios 2, 5 and 6 

ln = A\nax,n, (11.43) 
while for scenarios 3 and 4 we take 

ln^Alll{{(j))/MGUT).n- (11.44) 

For each epoch and scenario, we compute the evolution coefhcients Z„ as a weighted 
average over the measurements in the epoch. The representative redshift Zn is the 
average over the redshifts of observations inside the corresponding epoch. It is shown 
together with the resulting values for In in Tab. 111.11 This table summarizes our 
results under the assumption of proportionality. 



Rates of time variation in the present epoch 

For epoch 1 we incorporate the laboratory measurements for rates of varying couplings 
by linear extrapolation in time to the Oklo redshift zi ~ 0.14. The logarithmic time 
derivatives may be approximated by linear interpolation 

^ = d,lnGk^--^, (11.45) 

where ti = 1.8 x lO^y is the time corresponding to the redshift zi = 0.14. 



Method of averaging 

We evaluate the weighted average using all values listed in Tab. 110.11 This procedure 
may be quite problematic, since sometimes different observations are in manifest 
contradiction. We take the attitude that, given the possible presence of systematic 
effects both in spectroscopic determinations of nonzero coupling variations and in the 
primordial ^Li abundance, a viable model need not fit all data points. However, even if 
any given nonzero claimed variation is actually due to systematic error, we still expect 
the size of the error to be comparable to the size of the claimed variation. Thus, such 
claims are most conservatively interpreted as bounds on the absolute magnitude of 
variation. The surviving nonzero variation(s), in addition to the null bounds at other 
epochs, define a set of evolution factors which must be satisfied by any explicit model 
of evolution. 

For some scenarios we therefore also evaluate the evolution factors that are ob- 
tained by considering that some of the claimed observations of nonzero variation may 
instead be due to an underestimated systematic error. These alternative evolution 
factors are given in square brackets, corresponding to the following replacements: 
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Epoch 


1 


2 


3 


4 


5 


6 




0.14 


0.53 


1.6 


3.8 


10^ 


IQio 


Scenario 


h X 10« 


I2 X 10^ 


^3 X 10^ 


I4 X 10^ 


^5 X 10^ 


le X 10^ 


Q/ only 


— 01 ± 06 


— 1.1 ± 1.0 


— 26 ± 10 


—0 85 ± 37 


— 150 ± 350 


5 ± 34 


2 


-0.1 ± 0.1 


0.04 ±0.03 


-0.15 ± 0.08 


0.10 ± 0.03 


0.9 ± 14 


-0.37 ±0.20 


3 


4.1 ± 4.8 


— 1.5 ± 1.2 


42 ± 3 3 


—3.6 ± 0.9 


69 ± 920 


14 ± 8 


4 


3.9 ±8.5 


-3.4 ±2.7 


-8.4± 5.1 


-8.7±2.1 


31 ± 450 


33 ±21 


5, 


-0.02 ±0.18 


-0.24 ±0.18 


-0.25 ±0.10 


-0.61 ±0.13 


0.6 ±8.6 


1.7± 1.1 


(7 = 42) 












[0.4 ± 1.0] 


6, 


-0.02 ±0.12 


-0.10 ±0.07 


-0.17±0.07 


-0.44 ±0.10 


0.3 ±5.0 


1.2 ±0.6 


(7 = 70) 












[0.3 ±0.6] 


6, 


-0.12 ± 0.18 


0.04 ±0.12 


-0.30±0.11 


0.29 ±0.08 


0.7± 10 


-1±0.3 


(7 = 25) 








[-0.43 ±0.28] 







Table 11.1: Redshifts and evolution factors for each epoch, for the scenarios defined 
in Sec. 111.21 In the first row the values of In give the fractional variation of a; in 
Scenarios 2, 5 and 6 that of ax', and in 3 and 4 that of {(/>)/ Mgut- Values in brackets 
give, for BBN (1^) the evolution factors neglecting ^Li; or for I4, the evolution factor 
with the A/i//i value of |Reinhold06| substituted by that of [WendtOS] . 



Scenario 5, 7 = 42: Neglecting 'Li-abundance at BBN 
Scenario 6, 7 = 70: Neglecting Li- abundance at BBN 

Scenario 6, 7 = 25: Replacing the /i measurements of |Reinhold06j by the conserva- 
tive upper bound of [WendtOS] . 

In the case where a alone varies, since the fit including ^Li is poor we calculate a 
2(7 range using observational central values and errors of D and *He abundances as 
explained in Sec. 111.2.11 

11.3.3 Monotonic evolution with unification 

It seems natural to expect that variations of constants, if they occurred, evolve mono- 
tonicalljd. Looking on Figs. 111.11 to 111.71 and Tab. 111.11 we can ask the question 
whether the claimed variations are consistent with monotonic variation within the 
specific GUT scenarios. Here we briefly summarize whether the unified scenarios we 
consider can be consistent with a monotonic evolution of the single underlying varying 
parameter. 

Varying a only 

Although variation of a alone does not help to account for deviation of BBN abun- 
dances from standard theory, or for any nonzero variation of /i, the cosmic history 
is interesting due to the significant nonzero value in Epochs 3 and 4. The Oklo 
bound in Epoch 1 restricts the present time variation to 3.7 x 10^^^ y~^ (assuming 
no acceleration of dta). 



^In Sec. 112.21 we will also consider scenarios of quintessence wliere the implied variation is non- 
monotonic. 
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Scenario 2 

Scenario 2 favors a negative variation of ax at BBN, and a negative variation may 
also fit the Ma results. However, the Reinhold /z measurement indicates a positive, 
but much smaller, variation. The R/i results dominate the weighted average for Z4 
due to their small error on A In ax ■ The ratio A In /^/ A In a = 27 makes this scenario 
unlikely to fit the reported signal of nonzero Aa. 

Scenario 3 

In scenario 3 a positive variation of {(j))/McuT is favored by BBN. The high ratio 
A In ^/A In a ~ —325 makes the bounds obtained on a variation of /x strongly in- 
consistent with the claimed size of variation of a. The Reinhold et al. values again 
dominate the results for I4. 

Scenario 4 

In this scenario, the ratio A In/i/A Ina — —22 is again large and makes any observa- 
tion of significant nonzero A In a unlikely. Both the Ma and the R/i measurements 
point in opposite direction to BBN; however the two spectroscopic observations are 
also inconsistent with each other, within the scenario. Again, the R/i results dominate 
the determination of I4 due to the small error. 

Scenario 5, 7 = 42 

In this scenario the variation of ax favored by BBN is positive (Zg = (1.7 ± 1) x 
10'^), however both nonzero variations from spectroscopic data Ma and R/i require 
negative variations. With A In fx/A In a = —6 the spectroscopic measurements appear 
consistent with each other. Hence one would require some non-monotonic evolution 
to fit nonzero variations both at BBN and at moderate z. In Table fTOI we have also 
evaluated Iq using only the constraints given by D and "^He (in brackets). 

Scenario 6, 7 = 70 

As in the preceding scenario, BBN favors a positive variation in ax, but Ma and Kfi 
favor negative. Again, Fig. 111.61 may suggest a non-monotonic evolution. Fitting to 
BBN including '^Li we would obtain Iq = (1.2 ± 0.6) x 10"^; Tab. [TO also displays 
in brackets the value of Iq obtained from D and ^He bounds only. 

Scenario 6, 7 = 25 

In this scenario, both BBN and the Ma signal favor a negative variation of ax, 
whereas the R/i observations point towards a positive variation. Following the argu- 
ment of Wendt et al. (WendtOS] . we substitute the R/i value by the null constraint 
|A/x//z| < 2.5 X 10-^ [WendtOSj to obtain the bracketed value of U in Table fTTTl In 
this scenario the evolution factors show a crossover from negligible variation at low 
redshift, to strong and monotonically increasing negative variation at z ^ 2. 
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11.3.4 Tension between the ^Li problem and variation of /i 

Measurements of the primordial ^Li abundance show that the BBN abundance needs 
to decrease below the standard value to fit the observations, whereas the Reinhold fj, 
measurement indicates to increase at ^ ~ 3. We find that for all our unification 
scenarios the sign of the dependence on the fundamental parameter is the same for 
H and ""Li. Moreover, the coefficients of this dependence are nearly identical up to a 
common factor; hence the induced variations for fj, and ^Li point in the same direction, 
in contradiction to the tendency inferred from the observations. For example, for 
scenario 5 we find 

Aln/i= (23.2 - 0.657)Alnax, 
Aln^Li= (1692-497)Alnax. (11.46) 

These expressions change sign at 7 = 35.7 and 34.5, respectively. For a monotonic 
evolution, there is no possibility to have both a significant variation of /i and a vari- 
ation of opposite sign in the ^Li abundance. (In the regime 7 « 35 there is no 2a fit 
to BBN.) A similar result can be found for scenario 6 (including the SUSY partner 
mass dependence, which shows the same sort of degeneracy). Note that scenario 2 
and 3 are just limiting cases of scenarios 5 and 6. 

The main reason for this behavior is that variations of ^Li and /i are dominated 
by the variations of tti/Aqcd and mg/ Aqcd, respectively, with the same sign of 
prefactor. This degeneracy can be broken if me varies differently from the quark 
masses, a possibility that we do not consider in this thesis. For our scenarios with 
constant rh/rrie, the conflict between a monotonic time evolution and the fj,- and 
''Li-observations is reflected in the opposite signs of I4 and Iq. 

This observational tension for monotonic behavior is clearly depicted in Fig. I11.8[ 
where we plot simultaneously the averaged observational values of evolution factors 
^i/ln(l -|- Zi), normalized to ^4/ln(l -|- Z4). For Scenario 6, 7 = 25, we also display 
the result obtained by substituting the Wendt et al. value of /i variation for that 
of [Reinh old06] . The factor ln(l -I- Zi) is introduced as a convenient normalization 
to avoid compressing the scale of variations excessively in recent epochslf] For the 
purpose of a quick inspection we have omitted the error bars, which are of course 
necessary for a quantitative interpretation. 



^In quintessence-like theories, if the scalar field contributes a constant fraction of the total energy 
density of the Universe, as in so-called "tracker" models, the evolution of the field is typically also 
proportional to ln(l + z). This is an additional motivation for our normalization. 
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Chapter 12 

Probing quintessence models 



In this chapter we will describe how measurements of varying constants can be used 
to constrain models of quintessence under the assumption of grand unification. For 
that purpose we will use the measurements which we listed in Chap. \W\ in the six 
unified scenarios studied in the preceding chapter to constrain the free parameters of 
the quintessence models described in Chap.[5l Our aim is to see to what extent such 
models can be consistent with the behavior of time variations that we have outlined. 

12.1 Crossover quintessence 

As we have seen in Sec. 15.31 crossover quintessence models yield a monotonic evolution 
of the cosmon and hence also a monotonic variation of constants, assuming a constant 
coupling S to the fundamental varying parameter, 

A\naxiz)^S{ip{z)^ip{0)). (12.1) 

We discussed the viability of monotonic evolution in Sec. lll.3751 where a first judgment 
can be made by inspection of Fig. 111.11 to 111.71 for the various unification scenarios, 
or by inspection of Table 111.11 

To allow us to easily compare with the observational results, we observe that in 
Table [TTTT] the constraints for 1^ are considerably weaker than those for the remaining 
evolution factors. Furthermore, the one sigma range for is nonzero for all scenarios. 
Hence we shall compare observational and theoretical values for the ratios I1/I4, h/h, 
I3/I4 and le/U. We note the opposite sign of Z3 and I4 for scenario 2, 3 and scenario 
6 with 7 = 25, which disfavors any monotonic evolution for these scenarios. The 
averaged observational values in each epoch are given in Table 112.11 Note that the 
coupling S drops out of the ratios h/lj. Considering these ratios allows us to probe 
quintessence directly without knowing the absolute size of the coupling, only assum- 
ing its (approximate) constancy over the relevant range of evolution. In view of its 
monotonic evolution, crossover quintessence cannot give negative ratios h/lj. Hence 
it cannot be a good fit to the lithium abundance within the unification scenarios 2 to 
6 that we consider. This reflects the tension between the ^Li problem and a positive 
variation of fi discussed in Sec. 111.3.41 

For a simulation of the crossover quintessence model, we follow the procedure 
described in Sec. 15.31 By specifying the present densities of matter, radiation and 
dark energy and the model parameters w^q and 2:+, we can trace back the evolution 



98 



12.1. CROSSOVER QUINTESSENCE 



99 



Scenario 


^1 /u 

" ± / "4 


h/k 




h/l 


4 


/, K / l,A 

"O / "4: 





00 ± 01 


0.13 ±0.13 


0.31 ±0.18 


—592 ± 4031 


2 


—0 10 ± 12 


0.04 ±0.04 


-1.59 ±0.88 


—380 ± 228 


3 


—0 11 ± 13 


0.04 ±0.03 


-0.12 ±0.91 


—387 ± 241 


4 


-0.04 ±0.10 


0.04 ±0.03 


0.97 ±0.64 


-381 ± 259 


5, 7 = 42 


0.00 ±0.03 


0.04 ±0.03 


0.41 ±0.19 


-280 ± 191 


without BBN 












-66 ± 165 


6, 7 = 70 


0.00 ±0.03 


0.02 ±0.02 


0.39 ±0.17 


-275 ± 150 


without BBN 












-69 ± 139 


6, 7 = 25 


-0.04 ±0.06 


0.02 ±0.04 


-1.04 ±0.49 


-343 ± 142 


with Wendt 


0.03 ±0.05 


-0.01 ±0.03 


0.70 ±0.52 


231 ± 163 


Table 12.1: Ratios of the evolution factors from observations. 




Who z+ 


h/k 


h/k 


k/k 


k/k 






-0.95 3 


0.12 


0.38 


0.67 


26 






-0.99 3 


0.09 


0.28 


0.55 


38 






-0.9999 3 


0.02 


0.10 


0.37 


54 






-0.95 7 


0.15 


0.47 


0.80 


13 






-0.99 7 


0.15 


0.47 


0.80 


24 






-0.9999 7 


0.08 


0.27 


0.52 


125 






-0.999999 7 


0.02 


0.09 


0.34 


177 





Table 12.2: Ratios of evolution factors from crossover quintessence. The li are eval- 
uated by averaging over the variations evaluated at the same redshift as the data in 
each epochs (weighting by the number of absorption systems if appropriate). 



of quintessence and the other components of our Universe. In Tab. ll2.2l we display the 
ratios h/k expected from crossover quintessence for various parameters Who and z^., 
which can be compared with the observed ratios displayed in Tab. 112.11 Considering 
the lithium abundance to be affected by astrophysical systematics in scenarios 5 and 
6 (7 = 70), or using the null result for Aln/j, at intermediate redshift for scenario 
6 (7 = 25 "with Wendt"), we find that some crossover quintessence models indeed 
yield the observed order of magnitude for the ratios h/k- We conclude that crossover 
quintessence could, in principle, reconcile a coupling variation of the claimed size in 
epochs 3 and 4 with the bounds from late cosmology, i.e. epochs 1 and 2. This is due 
to the "slowing down" of the cosmon evolution, as noted in [Wetterich03| . Values of 
the present equation of state Who quite close to —1 would be required, however. In 
other words, an observation of coupling variations would put strong bounds on the 
dynamics of the cosmon fielcfl and provide for an independent source of information 
about the properties of dark energy. 

Note that observational probes of dark energy would not give results for Who that 
coincide with the values that we take in our model. Such probes do not actually 
measure the present-day equation of state, rather they extrapolate wq from past 
epochs under some parameterization. 



^The current WMAP5 bound on the equation of state parameter of quintessence is w = —0.972 it 
0.06 |WMAP5| . but the derivation assumes a constant w over a certain range of redshift, while our 
uiftO gives the equation of state at the present time. 
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12.2 Models with growing neutrinos and oscillating 
variation 

In this section we investigate the two models of growing neutrinos which were laid out 
in Sec. 15.41 These models do not obey the proportionality of all coupling variations 
for all redshifts and do not show a monotonic evolution of the cosmon field. A 
systematic analysis of all such models seems difficult, hence we concentrate on the 
specific examples given in jAmendolaOSj and |Wetterich08j . 

The combined variation is different for each of the unified scenarios. For all unified 
parameters except (0) /Mqjjt we still have a proportionality for the variations at all z, 
since the variations are proportional to Lp. The variations due to the direct coupling 
(|5.2ip can be described by our method of evolution factors Im even though the Z„ 
need not be strictly monotonic due to the oscillations in ip. However, for {(p) /Mqjjt 
we now have one variation linearly proportional to the variation of Lp, Eq. (|5.21|) and 
an additional one with a nonlinear dependence, Eq. (|5.30II5.3T|) . A simple treatment 
with common evolution factors for all variations will no longer be applicable. Due to 
the additional (/s-dependence of {(p) /A'Iqut we may have separate evolution factors for 
{(p) /Mqjjt, different from the (common) evolution factors for the other couplings. 

For example, the "linear contribution" (|5.2ip may dominate at BBN and induce a 
positive Iq common for all couplings. In the range 2 < 10 the "non-linear contribution" 
(I5.30|) could be more important, leading to effectively negative ^3,4 for {(p) /Mqjjt- 
(Such an effect could, in principle, relieve the tension between ^Li and a positive 
/i-variation at high z which was explained in Sec. 111.3.41 ) In practice, we calculate 
A In ax and A In j^ j^'^^^ at each epoch directly from the model, and extract the varying 
couplings and observables as explained in Chap. 1111 Then we may search for a set of 
parameters S, Rq which minimizes the for all measured variations. 

12.2.1 The stopping growing neutrino model 

The stopping growing neutrino model which was introduced in Sec. l5.4.T] has an oscil- 
lation in {(f)) that grows both in frequency and amplitude at late times as tp approaches 
its asymptotic value. Such oscillations must not be too strong as measurements be- 
tween z = 2 and today would measure a high rate of change. The oscillation may be 
made arbitrarily small by choosing small Rq. The restrictions from the low-z epochs 
are actually so strong that to a good approximation the non-linear contribution Rq 
can be neglected. However, the linear variation (|5.2ip is independent of Rq. It can 
be described by our method of evolution factors and yields for the set of parameters 
given in IWetterichOS] {ipt ~ 27.6, a = 10, e = -0.05) the ratios 



Comparing this with the numbers given in Table [12TT] shows that this model naturally 
yields evolution factors which are of the correct order of magnitude. We emphasize 
that no new parameter has been introduced for this purpose. 



h/h 
h/h 
h/h 
h/h 



0.008, 
0.09, 
0.44, 
175. 



(12.2) 
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12.2.2 Global fit to the scaling growing neutrino model 

Each growing neutrino model contains a few parameters that determine the cosmolog- 
ical evolution of the cosmon (p, and the Higgs v.e.v. {(j>). Two coupling parameters give 
respectively the relative strength of variation of ax with if, and the relative strength 
of the additional variation of the Higgs v.e.v. due to the varying triplet. For each 
example of cosmological evolution we may calculate the observables directly in terms 
of the two coupling parameters and make a global fit for their values. In performing 
the fit we take the 125 systems of the Murphy et al. a determination [Murphy03.2| 
within Epochs 3 and 4 (Eq. (|10.12p ) and further split them into 5 subsamples each 
with 25 absorption systems, since the data set extends over a wide range of redshift 
where there may be significant oscillations. 

For the global fits, we take the scaling growing neutrino model (Sec. I5.4.2P with 
(3 = —52, a = 10, m,y.o — 2.3eV. With zero variation at all times (no degrees of 
freedom), we find = 3.25 including ''Li at BBN, and = 2.40 neglecting "^Li. The 
results of the best fits with varying couplings are given in Table 112.31 



Scenario 


5 X 10-* 


i?o 






2 


-0.019 


0.045 


3.09 


0.16 


2 without Li7 


-0.040 





2.09 


0.31 


3 


1.64 





2.85 


0.40 


3 without Li7 


1.55 





2.03 


0.37 


4 


3.79 





2.85 


0.40 


4 without Li7 


3.80 





1.98 


0.42 


5.42 


0.30 





1.96 


1.29 


5.42 without Li7 


0.24 





1.87 


0.53 


6.70 


0.20 





1.93 


1.32 


6.70 without Li7 


0.16 





1.87 


0.53 


6.25 


0.18 


0.090 


2.73 


0.52 


6.25 without Li7 


-0.055 


0.061 


2.20 


0.20 



Table 12.3: Fitting parameters and minimal x^ values for the different unification 
scenarios for best fit to the scaling growing neutrino model [Amendola08j . The last 
column gives the increase in x^ produced when 5 and Rq are forced to vanish, i. e. for 
zero variation. 

It turns out that x^ cannot be reduced by more than 1.3 in the fit including ^Li 
and 0.53 in the fit neglecting ''Li, which we do not consider as convincing evidence 
for coupling variations within this model. We have investigated some other choices 
of parameters for the cosmological evolution and also the stopping growing neutrino 
model, without a substantial change in the overall situation. In view of the unsettled 
status of the observational data it seems premature to make a systematic scan in 
parameter space. Our investigation demonstrates, however, how a clear positive signal 
for a coupling variation could restrict the parameter space for quintessence models. 

Most of the additional variation of the Higgs v.e.v. occurs at later epochs, z < 2, 
thus recent observational bounds rule out any significant additional growth in (</)). 
We considered fitting the observational values excluding BBN, as a function of the 
model parameters 5 and i?o, and we find always that the value of i?o at the minimum 
of x^ is unobservably small. 
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12.3 Tests of the weak equivalence principle 



In Sec. 12.41 we have explained that variations of constants will influence the outcome 
of tests of the weak equivalence principle (WEP). Besides variations of couplings, the 
cosmon coupling to atoms also influences the outcome of tests of the weak equiva- 
lence principle, as it may also produce local gravitational effects which violate the 
WEP due to their interactions with matter. In Sec. 112.41 we use both aspects to 
set further constraints on a possible time variation of couplings. Test bodies with 
different composition have in general different couplings to the cosmon, and will 
hence experience different accelerations towards a common source. Usually, devia- 
tions from the universality of free fall are measured in terms of the Eotvos parameter 
r] ^Wetterich02.H,Dent06j 

rf-^ ^ 'P-^ , (12.3) 

\ab + ad 

where a^^c are the accelerations towards the source of the two test masses. The 
experiment setting the currently tightest limits on rj |Schlamminger07| has the result 

= (0.3± 1.8) X 10^^^ (12.4) 

for test bodies of Be {A = 9, Z = 4) and Ti {A = 48, Z = 22) composition, where the 
gravitational source is taken to be the Earth. 

In contrast to the direct observations of time varying couplings, tests of the uni- 
versality of free fall do not determine directly the values of fundamental 'constants' or 
their possible variations. However, given our basic assumption of a slow time variation, 
driven by a light scalar degree of freedom, the current limits on composition-dependent 
long range forces put bounds on the scalar couplings to different constituents of mat- 
ter. In our language, they measure or constrain the coefficients /3fe at z = 0, which 
relate the evolution factors and the changes in the cosmon held (see Eq. (|11.42p ). 

AlnGfc = 4/ = /3feA(^(z„) . (12.5) 

These constraints then imply bounds on the time variation of constants: Differenti- 
ating with respect to time, Eq. ()12.5p becomes 

^=f3k^- (12.6) 
Applying the conservative bound |Dent06| (p/Ho < 0.7, we derive 

< '/'max ^ 5 X 10-"y"' • (12.7) 

Pk can be derived as a function of r;, utilizing a model of nuclear masses and bind- 
ing energies and using GUT relations to reduce the number of free parameters (see 
|DSW08.^ for details on this treatment). It turns out that WEP violation places 
significant bounds on the present-day values of scalar couplings, as will be shown in 
the next section. 

The differential acceleration rj for two bodies with equal mass but different compo- 
sition and therefore different "cosmon charge" is related to the present time variation 
of couplings and cosmological parameters via [DSW08.^ 



12.4. BOUNDS ON PRESENT-DAY VARIATION 



103 



Scenario 


2 


3 


4 


5, 7 = 42 


6, 7 = 70 


6, 7 = 25 


F (Be-Ti) 


95 


-9000 


-165 


-25 


-26 


41 



Table 12.4: Values of F for a WEP experiment using Be-Ti masses for the unified 
scenarios studied in this thesis. 

Here a/ a is the present relative variation of the fine structure constant in units of 
year~^, Wh is the present dark energy equation of state and flh « 0.73 the present 
dark energy fraction. The "unification factor" F encodes the information on the 
specific relations between variations in fundamental parameters as implied by the 
different GUT scenarios, and the composition of the test bodies. The factor F is 
displayed in Tab. ll2.'ll for our different unified scenarios and for the Be-Ti test masses 
of |Schlamminger07| . For typical test mass compositions we find 1 < F < few x 10^. 
The very large value of F for Scenario 3 reflects an accidental cancellation of terms 
which we do not consider to be typical. 

Once F is fixed, the relation (|12.8p allows for a direct comparison between the 
sensitivity of measurements of ry versus the measurements of a/ a from laboratory ex- 
periments, or bounds from recent cosmological history, such as from the Oklo natural 
reactor or the isotopic composition of meteorites. 

12.4 Bounds on present-day variation 

Within our theoretical framework there exist three distinct ways to bound or measure 
the present-day rate of variation of fundamental parameters. The first is a direct 
measurement, of the type probed by atomic clock experiments (see Sec. 110. 7p . If 
one or more nonzero variations are found in this way, bounds on unified models may 
immediately be set. The second method is by combining information on the size 
of scalar field couplings from WEP tests (Section 112. 3p with a cosmological upper 
bound on the kinetic energy of scalar fields |Wetterich02.H IDent06| . Such bounds 
on scalar couplings will depend on the choice of unified model and in general will be 
independent of those derived from atomic clocks. Thirdly, under the assumption of 
a monotonic variation (that also does not significantly accelerate with time) , we may 
convert any "historic" bound on the net variation of a fundamental parameter since 
a given epoch into a bound on the present rate of variation: 

|Gfc|<(to-t„)"MGfe(to)-Gfc(t„)| = ^^, t,,<to. (12.9) 

Here to denotes the present, and any nonzero rate of variation should have the ap- 
propriate sign, i.e. Gk has the opposite sign to AGk referring to some past epoch. 

For any given unified model of time variations, the three bounds on present-day 
evolution will have different sensitivity. Therefore if one method gives a nonzero 
variation we would (in some cases) be able to distinguish between models due to the 
fact that the other bounds are still consistent with zero. To give a simple example, the 
direct detection of a nonzero time variation in atomic clocks near the present upper 
bound would immediately rule out a large class of models that cannot account for 
such a variation without leading to WEP violation above current bounds; and would 
also rule out models in which such nonzero variations extrapolated to past epochs 
tn < to would exceed observational bounds. 
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However, this chain of inference does not function equally well in all directions. A 
nonzero finding of differential acceleration violating the WEP would indicate nontriv- 
ial scalar couplings, but need not imply nonzero time variation since the rate of change 
of the scalar is not bounded from below. Also, a nonzero variation at some past epoch 
tn would not necessarily imply a lower bound to the present-day rate of variation or 
size of scalar couplings, since the variation could have slowed substantially since then 
(either due to nonlinear scalar evolution or a nonlinear coupling function). Only with 
the assumption of a reasonably smooth and monotonic variation of the scalar field 
and its coupling functions, one can find, for any given unified model, where the first 
signals of present-day or recent variation are expected to appear. 

At present these methods give null results up to redshifts of about 0.8, but if 
a nonzero time variation exists, we can determine for each unified scenario which 
observational method is most sensitive. Thus if a nonzero signal of late time variation 
arises it may be used to distinguish between models. We assume for this purpose an 
approximately linear variation over recent cosmological times, thus measurements of 
absolute variation at nonzero redshift z imply time derivatives 

d\nX AlnA:(z) 

— ^ 7-^- 12-10 

dt to - t{z) ^ ' 

Here X is the fundamental varying parameter: we consider first AT = a, if only the 
fine structure varies; in scenario 1 X = G^mj^, in scenarios 2, 5 and 6 A = ax and in 
scenarios 3 and 4 A = {<j>) / Mcut ■ Then Table [1231 gives the precision of bounds on 
time derivatives for the unified scenarios we consider, except scenario 1 (varying Gn) 
which is probed by a quite different set of measurements. As explained in Sec. 110.41 
we take the Oklo bound as applying directly to the variation of a, and increase its 
uncertainty by a factor 3 to account for possible cancellations when other parameters 
also vary. We present the recent Rosenband et al. [RosenbandOSj Al/Hg ion clock 
bound separately to illustrate to what extent it improves over previous atomic clock 
results. 

Extending this method beyond z ^ 0.5 becomes questionable. One could use 
linearity in ln(l + z) instead of t, but even this improvement may lead to unreliable 
extrapolations for models with a particular dynamics of the scalar field, as crossover 
quintessence. 
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Scenario 


X 


Clocks 


Al/Hg 


WEP 


Oklo a 


Meteorite 


Astro 


a only 


a 


0.13 (a) 


0.023 


6.2 


0.033 


0.32 


0.44 (y) 


2 


ax 


0.074 ill) 


0.027 


0.007 


0.12 


0.015 


0.006 (NH3) 


2S 


ax 


0.12 ip.) 


0.044 


0.012 


0.19 


0.026 


0.010 (NH3) 


3 


m/Mx 


2.6 i^i) 


12.4 


0.33 


54 


0.53 


0.22 (NH3) 


4 


m/Mx 


6.2 i^i) 


1.78 


0.35 


7.7 


1.2 


0.51 (NH3) 


5, 7 = 42 


ax 


0.32 (a) 


0.024 


0.013 


0.11 


0.069 


0.035 (NH3) 


6, 7 = 70 


ax 


0.21 (a) 


0.016 


0.008 


0.070 


0.049 


0.025 (NH3) 


6, 7 = 25 


ax 


0.25 ifi) 


0.027 


0.011 


0.12 


0.056 


0.021 (NH3) 



Table 12.5: Competing bounds on recent (z < 0.8) time variations in unified scenarios. 
For each scenario we give la uncertainties of bounds on d{lnX)/dt in units 10~^^y~^, 
where X is the appropriate fundamental parameter. For the Oklo bound we multiply 
the uncertainty by a factor 3 except when only a varies. The column "Clocks" indi- 
cates whether a or /x gives the stronger bound; the recent Al/Hg limit |Rosenband08| 
is given a separate column. The column "Astro" indicates which measurements of 
astrophysical spectra are currently most sensitive in each scenario. 



Chapter 13 

Conclusion and outlook 



We have developed a systematic method to relate variations of fundamental parame- 
ters of particle physics to the primordial isotope abundances produced by BBN. The 
main advantage of the method, which is laid out in part two of this thesis, is that 
we are able to vary every parameter independently, both at the level of fundamental 
Standard Model parameters and of nuclear physics parameters, thus we are not depen- 
dent on any particular theoretical model which enforces particular relations between 
the variations. 

We follow a two step approach, first extracting the nuclear parameter dependence 
(without major theoretical uncertainties) and in a second step relating this to vari- 
ations of fundamental Standard Model parameters. We define two linear response 
matrices, where the first, C, encodes the change in predicted abundances produced 
by small variations away from the current values of nuclear physics parameters which 
enter the BBN integration code. These parameters comprise the gravitational con- 
stant, fine structure constant, neutron lifetime, electron, proton and neutron masses, 
and binding energies of ^ < 7 nuclei. The dependences of nuclear reaction rates on 
these parameters are also implemented insofar as they are calculated within some 
efi'ective theory. One notable result is that the ^Li abundance depends heavily on the 
binding energies of ^He, ^He and ^Be. 

We also investigated possible further effects of variations in nuclear reaction rates 
on predicted abundances by varying each rate {i. e. thermal integrated cross-section 
(av)) separately by a temperature-independent factor. We find that the *He abun- 
dance is insensitive to nuclear rates, and only eight reactions could lead to significant 
variation of the D, "^He or ""Li abundances. 

The second response matrix, F, relates variations in nuclear parameters to the 
fundamental parameters of particle physics, comprising the gravitational constant, 
fine structure constant, Higgs vacuum expectation value, electron mass, and the light 
(up and down) quark masses. At this point theoretical uncertainty enters into the 
relation between quark masses and nuclear binding energies. We parameterize the 
dependence of binding energies on the pion mass (and hence on light quark masses) 
by the deuteron binding, which has been treated by a systematic expansion in effective 
field theory. 

The resulting fundamental response matrix R = CF allows us, first, to bound 
the variations of the six fundamental couplings individually, some bounds being at 
the percent level. Secondly, studying three exemplary unified scenarios, we can also 
bound correlated variations affecting many couplings at once. We find that one sce- 
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nario allows us to fit observed D, ^He and ^Li abundances within 2a bounds, given a 
variation Aa/a = —2 x 10~^ away from the present value; another fits these obser- 
vational abundances within almost la bounds, given a variation Aa/a = 4 x lO"**. 
For this analysis we have left the linear matrix approach and make use of our BBN 
code which allows us to derive the full nonlinear dependence. 

In a next step (part three of the thesis), we combine our findings for BBN with 
further measurements of the variation of fundamental constants from BBN to today. 
Within grand unified theories, this set of different varying parameters can be consis- 
tently reduced to a variation of a few "unification parameters" , namely the unification 
scale Mgut/Mp, gauge coupling ax, the Fermi scale {4))/Mgut and SUSY-breaking 
masses rh/MGUT- We define various GUT-scenarios for varying couplings by the 
assumption of proportionality of fractional variations of the unification parameters. 

Assuming that couplings really vary, this is a way of excluding such GUT scenarios 
by demanding consistency of the implied variations. The assumption of proportion- 
ality permits us to project all observations into constraints on a common evolution 
factor l{z) for each scenario. We show that different GUT scenarios yield different 
time evolutions of l{z) assuming that certain claimed measurements of varying con- 
stants are correct. We confirm that "simple" models which have only one fundamental 
parameter varying [ax or (0) /Mqjjt) result in inconsistent variations. However, com- 
bined variations of these two parameters, as described in scenarios 5 and 6, lead to 
results more consistent with the possible quintessence-induced time variations of fun- 
damental couplings. For instance, some models of crossover quintessence and models 
of growing neutrinos naturally yield ratios of evolution factors which are of the same 
size as those derived from measurements of varying couplings. 

Still, we have not found a scenario with a monotonic time evolution l[z) that 
makes the two main signals or hints of variation ( |Murphy03.2| , [Reinhold06j ) and 
BBN mutually consistent. A monotonic evolution requires either to discount one of 
the "signals" by substantially increasing its uncertainty, or to alter our assumptions 
by including additional time variation of some Yukawa couplings. 

In a last step we demonstrate how a clear observation of time variation of funda- 
mental couplings would not only strongly disfavor a constant dark energy, but also 
put important constraints on the time evolution of a dynamical dark energy or quint- 
essence. We have shown that for a given unified scenario, the bounds on the time 
variation of various couplings in different cosmological epochs can strongly restrict 
the possible time evolution of the cosmon field and put very strong bounds on late- 
time dynamics of the cosmon field. Hence, once at least one irrefutable observation 
of some coupling variation at some redshift becomes available, our method provides 
a powerful tool for testing extensions of the Standard Model or, vice versa, allows us 
to control consistency of claimed variations under the demand of unification. 

We have demonstrated this by an analysis that implicitly assumes a nonzero varia- 
tion, considering both general features and specific quintessence models. However, we 
are aware that the actual values for the evolution factors l{z) from this analysis may 
be premature, since the observational situation is unclear and on moving grounds. For 
example, taking the recent reanalysis of the variation of the proton to electron mass 
ratio /i in Ref. 'KingOSj instead of the results in Ref. [Reinhold06| used in this thesis, 
would strongly influence the values of the evolution factors. We have demonstrated 
this in a somewhat different way by investigating the change in the evolution factors 
if some claimed observations of varying couplings are omitted. 
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Outlook 

Progress in the field of BBN requires both observational and theoretical improve- 
ments. Both statistical and systematic errors in abundance measurements could be 
improved, for example observations to better determine the nature of systems where 
"^He is measured SteigmanOS, , or stellar modeling to test possible solutions of the 
^Li problem. On the theoretical side the relation between quark masses and nuclear 
physics remains unclear beyond the level of the two-nucleon system: the largest un- 
certainty in our BBN bounds arises from the poorly known dependence of the binding 
energies on the fundamental couplings. BBN is already the most powerful probe of 
fundamental "constants" in the early Universe, and precision bounds may well be 
obtained, given continued efforts in observation and theory, to rule out or confirm the 
presence of a cosmological variation. 

Our investigation has further shown how the variations of different couplings in 
the Standard Model may be compared. If the observational situation becomes clearer 
and at least one nonzero time variation is established, such methods may be used for 
new tests of the idea of grand unification and, even further, models of quintessence. 
The presented method could then easily be applied to constrain both theories beyond 
the Standard Model like grand unification and theories beyond the concordance model 
of cosmology like quintessence. Hence, it can establish relations between two a priori 
distinct theories, namely high-energy physics and general relativity, a capability which 
is so far dominated by theories with radically new concepts of physics like string 
theory. 
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Appendix A 

Conventions 



A.l Symbols and abbreviations 



Symbol 


Explanation 




fine structure constant 


as 


strong coupling constant 


ax 


grand unified coupling constant 


Bi 


binding energy of nucleus i 




light quark mass difference 


V 


baryon-to-photon ratio; Eotvos parameter 


Gn 


Newton's constant 




QCD scale 


m 


average light quark mass 


md 


d- quark mass 


me 


electron mass 


rUn 


neutron mass 




average pion mass 


rrip 


proton mass 


rus 


strange quark mass 


m„ 


u-quark mass 


niN 


average nucleon mass, rnjv = 5('Ti.„ + mp) 


Mqut 


GUT scale 


Mp 


Planck mass 


Mp 


reduced Planck mass 


Mx 


GUT scale (= Mgut) 




proton to electron mass ratio, := mp/rrie 




baryon fraction 


Qn 


neutron proton mass difference, Qn '■= mn — nip 



Table A.l: List of symbols 
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Abbreviation 


Explanation 


BAO 


baryon acoustic oscillations 


BBN 


Big Bang Nucleosynthesis 


EFT 


effective field theory 


CMB 


cosmic microwave background 


DOF 


degree of freedom 


MSSM 


minimal supersymmetric standard model 


QCD 


quantum chromodynamics 


QED 


quantum electrodynamics 


SBBN 


standard big bang nucleosynthesis 


SM 


Standard Model (of particle physics) 


SN 


supcrnovae 


v.e.v. 


vacuum expectation value 


WEP 


weak equivalence principle 
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